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ABSTRACT 

This  thesis  investigates  the  feasibility  of  adapting  Naval 
Architecture  and  Marine  Engineering  design  processes  to  digital  computer 
solution© 

Two  problems,  the  determination  of  critical  whirling  frequencies 
of  propeller  shafts  and  the  determination  ^f  stresses  in  a  transverse  web 
frame  are  programmed  for  a  digital  computer.  A  detailed  explanation  of 
the  programs  and  a  discussion  of  the  methods  of  approach  and  difficulties 
encountered  in  adapting  the  solutions  of  these  problems  to  a  digital 
computer  is  included.  The  most  important  results  obtained  from  programming 
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these  problems  are  experience  with  and  knowledge  of  digital  computers. 
Applying  this  knowledge  and  experience  to  the  general  field  of  naval 
architecture  the  following  conclusions  and  recommendations  are  presented! 

Conclusions: 

a.  Digital  computers  are  applicable  to  naval  architecture  problems. 

b.  The  average  naval  architect  has  sufficient  background  to  pro- 
gram a  digital  computer. 

c.  A  programmer  should  have  a  complete  engineering  understanding 
of  the  problem  and  some  familiarity  with  numerical  methods  of 
mathematical  processes. 

d.  In  general  the  limitations  imposed  by  the  computer  on  a 
problem  are  no  more  restrictive  than  those  imposed  by  manual 
methods. 

e.  The  digital  computer  used  should  be  able  to  handle  numbers 
usually  encountered  in  engineering  problems. 

Recommendations : 

a.  Digital  computer  should  be  used  in  naval  architecture, 

b.  Digital  computers  are  suited  for  theses  and  research  problems. 

c.  Any  future  thesis  should  only  consider  the  computer  as  an  aid 
in  the  solution  of  a  specific  problem  and  not  general  appli- 
cability of  the  computer. 

d.  For  engineering  applications  the  programmer  should  have;  (a)  a 
good  engineering  understanding  of  his  field;  (b)  knowledge  of 
numerical  methods  of  mathematical  processes;  (c)  organizational 
abilities. 


Thesis  Supervisor:  Robert  J.  Hansen 

Title:  Associate  Professor  of  Structural  Engineering 
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Introduction 

The  field  of  naval  architecture  contains  some  of  the  most 
complex  problems  facing  the  engineering  profession  in  this  modern  day 
and  age,  A  ship  is  a  structure  containing  a  network  of  mechanical 
equipage  that  is  unequaled  by  any  other  movable  structure  in  the  world0 
The  complexity  of  such  a  structure  almost  defies  analytic  approach  and 
as  a  result  much  of  the  design  work  is  based  on  previously  built  vessels. 
A  parent  design  is  altered  by  repeating  various  numerical  procedures 
until  the  desired  result  is  obtained.  The  practice  of  naval  architecture 
can  best  be  described  as  laborious  and  repetitive  with  a  necessity  for 
a  background  of  experience. 

In  an  effort  to  eliminate  the  necessity  of  utilizing  personnel 
to  perform  laborious  calculations,  following  World  War  II  a  committee  was 
appointed  by  the  Bureau  of  Ships  to  investigate  and  report  on  the  appli- 
cability of  utilizing  machine  aided  calculations  for  problems  arising  in 
the  preliminary  design  stage.  Due  to  the  stage  of  development  of  elec- 
tronic computers  at  that  time  an  analog  computer  was  developed  for 
stability  calculations  in  the  preliminary  design.  Although  this  Investi- 
gation did  consider  the  application  of  digital  computers  to  problems  in 
the  preliminary  design  stage,  the  application  to  the  general  field  of 
naval  architecture  was  not  fully  considered,  (ll) 

The  purpose  of  this  thesis  is  to  investigate  the  applicability 
of  digital  computers  to  the  field  of  naval  architecture,  with  primary 
emphasis  on  structural  design  problems,  to  acquaint  the  authors  and  other 


1.  Numbers  in  parentheses  refer  to  the  Bibliography  at  the  end  of  the 
thesis. 
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naval  architects  with  this  new  method  of  approach,  to  determine  the  type 
problems  best  suited  for  computer  solution  and  publish  some  of  the  diffi- 
culties encountered,  to  determine  the  type  of  person  required  to  utilize 
most  effectively  the  machine  aids  and  to  incite  interest  in  applying 
digital  computers  in  the  investigation  of  problems  heretofore  considered 
too  laborious  and  time  consuming  to  analyze  by  hand  calculations.  To 
do  this  the  authors  had  to  first  acquaint  themselves  with  an  existing 
computer  and  secondly,  to  pick  sample  engineering  problems  of  generic 
nature  for  specific  programming. 

The  authors  are  graduate  students  in  naval  architecture  whose 
studies  have  been  primarily  channeled  in  the  direction  of  hull  structures. 
The  background  in  mathematics  included  calculus  through  differential 
equations  and  a  basic  course  in  vector  and  complex  mathematics.  Their 
background  in  mathematics  is  no  better  than  can  be  expected  from  the 
practicing  naval  architect,. 

Two  problems  were  picked  as  samples,  one,  calculation  of  the 
whirling  frequencies  of  a  propeller  shafting  arrangement  and  two,  an 
analysis  of  an  indeterminate  structure  of  five  degrees  of  indeterminancy. 
These  problems  were  picked  because  they  are  examples  of  types  that  require 
reiterative  procedures  in  their  solutions.  Further,  both  riroblems  had 
solutions  that  had  been  previously  determined  by  hand. 

The  digital  computer  investigated  and  utilized  was  the  Whirlwind  I 
computer  at  MIT.  This  computer  started  operations  in  November  195>0  and 
has  been  used  as  a  guide  for  the  development  of  many  commercially  available 
machines  today.  For  this  reason  the  general  conclusions  drawn  from  an 
analysis  of  the  applicability  of  Whirlwind  I  should  apply  to  any  machine 
considered  for  particular  problems,  although  other  machines  may  have  different 
capabilities. 
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Basically,  all  digital  computers  utilize  the  same  theory  of 
operation.  They  all  contain  the  following  elements t     storage  element, 
control  element,  arithmetic  element  and  input  and  output  devices,  as 
depicted  in  Figure  I. 

Calculations  are  performed  in  the  arithmetic  element  utilizing 
the  numbers  and  instructions  contained  in  the  storage  element.   Input 
and  output  devices  are  used  to  get  information  in  and  out  of  the  computer, 
The  control  element  makes  the  components  function,. 


Arithmetic 
Element 

^ 

Control 
Element 

Storage 
Element 

V 

r 

> 

i 

r 

, 

> 

4 

\ 

> 

t 

a 

i 

1     DIGIT 

TRAN5M  ISSIO 

N       BUS 

1 

r 

T  ■■■'-«  .    .  -~ 

r\  i  .  -i-  i^  .  .  -i- 

4 

r 

1 

.  in  r 

"i 

J   1 

Ol 

>  i 

ru 

i 

^s 

Schematic  Diagram  of  Digital  Computer 
Figure  I 
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The  discussion  of  the  individual  elements  that  follows  pertains 
to  Whirlwind  I  only,  although  other  computers  have  similar  equipment  that 
performs  the  same  functions. 

The  storage  element  consists  of  20I4.8  registers  or  "cells." 
Instructions  or  numbers  are  stored  in  these  registers  by  impressing 
electrical  pulses  on  magnetic  cores.  Each  register  contains  16  such  cores 
and,  therefore,  the  basic  word-length  is  16  digits.  The  pulses  may  be 
either  positive  or  negative  and  the  resulting  direction  of  magnetization 
determines  whether  the  digit  represented  is  0  or  1.  This  means  that 
all  the  machine  claculations  are  carried  out  in  the  binary  number  system. 
In  every  register  the  first  binary  digit  represents  the  sign  of  the  number 
and  the  remaining  IS   binary  digits  can  represent  approximately  5>  decimal 
digits. 

Instructions  are  stored  In  the  registers  in  a  slightly  different 
manner.   Here  the  first  five  binary  digits  represent  the  mathematical 
operation  that  is  to  be  performed  when  the  instruction  is  obeyed.  The 
remaining  eleven  binary  digits  represent  the  location  of  the  register 
containing  the  number  on  which  the  operation  is  to  be  performed. 

The  permanent  storage  of  Whirlwind  I  may  be  increased  to 
approximately  22,000  registers  by  the  use  of  a  magnetic  drum.  This  device 
is  somewhat  sloxirer  in  the  transfer  of  information  than  is  the  permanent 
storage,  but  it  is  very  convenient  for  the  solution  of  long  problems. 

The  arithmetic  element  consists  of  three  registers.   In  one  of 
these,  called  the  accumulator*  or  "AC"  register,  the  basic  arithmetic 
operations  of  addition  and  subtraction  are  performed.  Multiplication  and 
division  are  accomplished  by  a  series  of  these  two  basic  operations. 
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The  input  devices  are  a  photoelectric  tape  reader  and  a  mechanical 
tape  reader.  These  devices  read  into  the  storage  element  coded  information 
from  a  punched  paper  tape.  They  convert  holes  in  the  tape  into  electrical 
signals  to  properly  magnetize  the  registers  of  the  storage  element.  A 
conversion  routine  is  available  to  transform  decimal  information  on  the 
tape  into  binary  information  for  storage  in  the  computer.  This  feature 
is  a  great  aid  in  programming  and  enhances  the  flexibility  of  the  computer. 

The  output  devices  available  are  flexowriters,  five  magnetic 
tape  units  and  an  oscilloscope  and  camera.  The  desired  output  information 
can  be  typed  out  directly  by  flexowriters,  or  it  can  be  stored  on  magnetic 
tape  and  typed  out  by  the  flexowriter  at  a  later  time.  The  latter  procedure 
saves  computer  time  since  it  is  much  faster  to  read  information  onto 
magnetic  tape  than  to  type  it.  Ano-ther  method  of  getting  the  information 
out  is  to  display  it  on  the  oscilloscope  and  photograph  it.  This  saves 
typing  time,  but  the  film  must  be  processed.   However,  it  is  convenient 
for  reproducing  the  results.  Standard  output  programs  are  available  or 
special  ones  may  be  written  to  get  the  information  out  in  the  form  desired. 
For  instance,  numbers  may  be  in  octal  or  decimal  system;  written  in  columns, 
lines,  blocks,  or  tabulated;  formed  raised  to  a  fixed  power  of  ten  or  with 
a  fixed  decimal  point.  This  flexibility  of  output  is  an  important  feature 
of  digital  computers. 

Control  Element 

Each  of  the  various  components  of  the  computer  consists  of 
combinations  of  electrical  circuitry  and  transfer  buses.   Only  through 
intelligent  control  of  these  devices  can  electrical  signals  be  made  to 
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perform  arithmetic  operations.  The  control  element  performs  this  function 
by  causing  the  flow  of  magnetic  signals  to  reach  the  proper  element  at  the 
appropriate  time.  Basically  then,  the  control  element  is  a  switching 
device  that  channels  the  information  through  transfer  buses  to  itself 
and  the  other  elements.  To  do  this,  the  control  element  must  interpret 
the  program  instructions,  which  are  themselves  magnetic  signals,  into 
other  magnetic  signals  which  control  the  operation  of  the  computer. 

The  speed  of  operation  of  the  computer  depends  upon  the  memory 
access  time  and  the  time  required  to  integrate  a  number  of  switching  com- 
binations into  a  useful  operation.  Switching  is  performed  by  utilizing 
electrical  impulses  from  a  high  frequency  oscillator.  A  new  set  of 
switching  combinations  can  be  established  for  each  time  pulse  received 
by  the  control  element.  The  time  required  to  perform  an  operation  depends 
on  the  number  of  impulses  necessary  to  perform  the  switching  and  the 
repetition  rate  of  the  impulses.  In  Whirlwind  I  the  pulse  repetition  rate 
is  1  micro-second.  At  this  rate  the  addition  of  two  numbers  requires  It 8 
micro-seconds  and  multiplication  requires  approximately  6£  micro-seconds. 
Using  the  conventional  rating  system  which  is  the  number  of  additions  a 
computer  can  perform  per  second,  the  Whirlwind  I  is  rated  at  20,000  opera- 
tions per  second. 

The  instructions  necessary  to  perform  the  desired  mathematical 
calculations  are  called  a  program  and  are  fully  explained  in  Appendix  A. 
For  each  computer  there  exists  a  standard  instruction  code.  By  ingenious 
combinations  of  these  standard  instructions  the  programmer  can  cause  a 
computer  to  perform  the  requisite  numerical  calculations. 
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One  of  the  more  important  instructions  is  the  "conditional 
response."  This  instruction  allows  the  computer  to  make  an  either-or 
decision  based  on  the  sim  of  the  number  in  the  accumulator.  As  a 
result  of  this  capability  the  computer  can  be  made  to  cycle  back  through 
a  series  of  instructions  any  desired  number  of  times.  The  cycling  ability 
is  ideally  suited  to  repetitive  processes  which  are  useful  in  design  calcu- 
lations. 

Digital  computers  have  been  widely  used  in  other  fields  of 
engineering  and  new  adaptations  are  being  developed  every  day.  Some  of 
the  more  important  uses  developed  are: 

(1)  Aircraft  weight,  balance  and  thrust  calculations. 

(2)  Steam  turbine  design 

(3)  Catapault  launching  analysis 

(k)     General  structural  application  in  the  aircraft 
industry. 
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PROCEDURE 
General 

In  order  to  accomplish  the  objectives  discussed  in  the 
introduction  it  was  necessary  that  some  knowledge  of  computers  be  obtained, 
A  twenty- five  hour  course  in  programming  for  Whirlwind  I  was  completed 
by  the  authors  to  partially  fulfill  this  requirement. 

Once  the  capabilities  of  the  computer  were  realized,  it  was 
necessary  to  decide  what  problem  would  best  prove  the  basic  contentions. 
Probably  the  biggest  single  factor  in  this  decision  was  time.  A  very 
complicated  problem  would  require  a  great  deal  of  time  to  program  while 
a  simpler  problem  would  not  test  the  authors'  ability. 

A  compromise  was  reached  by  deciding  to  program  two  problems. 
The  first  should  be  a  relatively  short  and  uncomplicated  problem  aimed  at 
giving  the  authors1  experience  with  the  computer  and  also  making  some 
contribution  to  naval  architecture.   The  second  problem  should  be  of  such 
size  and  complexity  that  its  solution  would  more  fully  describe  the 
abilities  of  the  computer  and  the  benefits  of  computer  applications  to 
naval  architecture. 

The  first  problem  programmed  was  the  determination  of  the 
whirling  frequencies  of  propeller  shafts.  The  method  of  solution  was 
developed  by  Professor  F.  M.  Lewis.   This  problem  was  selected  because 
the  method  of  solution  is  an  iteration  process  ideally  suited  to  computer 
application.  Furthermore,  the  problem  was  short  and  each  step  had  been 
delineated  by  hand  calculations  performed  by  Professor  Lewis.  This  last 
reason  was  very  important  because  it  enabled  the  authors  to  check  their 
work  and  to  ascertain  the  accuracy  of  the  program. 
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The  second  problem  programmed  was  the  determination  of  stresses 
in  a  transverse  web  frame.  A  procedure  similar  to  the  one  used  in  the 
Bureau  of  Ships  was  adapted  to  the  computer.  The  numerous  computations 
encountered  in  this  problem  made  a  computer  solution  very  desirable.  The 
calculations  were  ms.de  for  a  new  destroyer  type  for  which  a  partial  check 
of  the  numerical  results  was  available.  Another  important  reason  this 
problem  was  selected  was  to  further  the  authors1  knowledge  in  structural 
analysis  of  a  ship.  The  course  of  instruction  in  naval,  architecture  had 
not  included  the  calculation  of  transverse  strength. 

PROPELLER  SHAFT  WHIRLING  FREQUENCIES 

Description  of  the  Problem 

Basically,  the  problem  is  to  calculate  the  whirling  or  lateral 
frequency  of  a  propeller  shaft  system  in  which  the  excitation  is  induced 
by  the  propeller.  The  method  of  calculation  involves  a  numerical  solution 
of  a  fourth  degree  differential  equation  of  the  form 

k  2 

d  y   mw 

ji  -  w?  -  ° 

dx 

of  which  the  general  solution  is  of  the  form 


y  =  C  sin  Xx  +  C  cos  Xx  +  C  sinh  Ax  +  C.  cosh  Xx 


where 
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let  jx  s.  9  and  y"  «*  Y.     Applying  boundary  conditions  to  determine  the 
constants  of  integration:     at  x  =  0 
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O.cos  p  cosh  p-1,    k  9Q(coS  P  sinh  P-sin  R  cosh  p) 
9L  =  T^    sinh  p-sin^     }  +  sinh  p-sta  p 
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From  which 


o      cos  B  sinh  B-sin  B  cosh  ft 

1  '        sinh  B-sin  B 

o      cos  B  cosh  B-l 

2  "  "  B(sinh  B-sin  8) 

2B  sin  B  sinh  B 

S. 


3     sinh  B-sin  B 

The  above  equations  are  applicable  to  one  span  of  the  shaft. 
In  the  solution  of  this  problem  these  equations  were  used  to  determine  the 
moment  and  angle  of  rotation  at  the  end  of  one  span  and  then  carry  these 
results  over  to  the  next  span.  The  number  of  spans  to  be  used  is  unspeci- 
fied, however,  it  is  the  opinion  of  the  authors  that  four  spans  give 
sufficient  accuracy. 

Since 

various  values  of  B  are  assumed  and  the  residual  moment  at  the  free  end 
plotted  versus  the  frequency.  The  value  at  which  the  end  moment  becomes 
zero  is  a  resonant  point  or  natural  whirling  frequency  of  the  shafting 
system. 

In  applying  these  equations  it  is  useful  to  determine  certain 
parameters  in  terms  of  the  physical  properties  of  the  shaft  system.   These 
are  given  below  and  the  notation  is  explained  in  Figure  I. 
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Schematic  Diagram  of  Shaft 
Figure  I 

It  can  be  shown  from  the  equations  previously  developed  that 


where 
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and 


J  =  polar  moment  of  inertia  of  propeller 
and  entrained  water* 


q  =  the  order  of  vibration  being  investigated. 
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The  number  subscripts  refer  to  spans  as  shown  in  Figure  I, 
From  the  above  it  can  be  seen  that 


where 


912  "  "  (S11901  +  m01S21/cal) 


cal  -  Wl 


Here  the  second  number  subscripts  refer  to  the  span.  It  should 
be  noted  that  0   is  assumed  as  unity. 

For  determining  s,  ,  s  ,  s  for  different  spans  it  is  convenient 
to  use  the  following: 


*  ■  '•  V* 


where 

E  I 
_  1  1 

■  =  W  L3- 


c 


1  1 


The  procedure  described  is  repeated  for  each  end  moment  and 
angle  of  rotation  until  the  residual  end  moment,  m_j  is  determined.  Then 
a  new  value  of  B  is  assumed  and  the  process  repeated  until  a  series  of 
end  moments  versus  frequencies  results  are  generated.  By  plotting  these 
values,  or  by  interpolation,  the  critical  whirling  frequency  can  be  found. 
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General  Discussion  of  Program 

The  method  developed  for  solution  of  the  propeller  shaft 
whirling  frequencies  is  ideally  suited  to  computer  solution  because  it 
is  essentially  one  of  numerical  iteration.  Various  values  of  |3_  are 
assumed  and  the  residual  moment  calculated  for  any  number  of  spans.  Since 
the  labor  involved  in  programming  is  the  same  regardless  of  the  number 
of  times  a  process  is  repeated,  good  accuracy  can  be  obtained  with  no  extra 
effort  by  the  programmer  and  small  increase  in  cost  of  operation  of  the 
computer. 

The  underlying  guide  in  programming  this  problem  was  to  make 
it  as  general  as  possible.  This  was  done  so  that  anyone  wishing  to  use 
this  method  in  the  future  would  merely  have  to  change  a  few  numbers  to 
adapt  this  solution  to  his  particular  problem.  The  quantities  to  be 
changed  were  grouped  separately  at  the  beginning  of  the  program  and 
involve  such  factors  as  the  physical  characteristics  of  the  propeller 
and  shafting,  the  number  of  spans  to  be  considered  and  the  order  of  the 
vibration  being  investigated. 

The  first  step  in  programming  was  to  assign  floating  addresses 
or  "flads"  to  registers  which  were  to  hold  the  numerical  quantities  and 
key  instructions.  These  flads  contain  the  numbers  or  instructions  upon 
which  mathematical  operations  are  to  be  performed.  The  purpose  of 
floating  addresses  is  to  simplify  programming  by  enabling  the  programmer 
to  refer  to  a  flad  rather  than  an  absolute  address  each  time  an  instruc- 
tion is  written.  The  advantage  of  this  system  is  clearly  demonstrated 
in  the  programs  contained  in  the  appendix.  The  flads  used  in  pro- 
gramming and  the  numbers  they  represent  are  given  in  Appendix  B, 
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In  addition  to  storing  the  data  and  constants  of  any  problem, 
registers  must  be  assigned  to  store  intermediate  results  and  for  temporary 
storage. 

The  addresses  of  key  instructions  were  assigned  flads  to  aid  in 
writing  the  program.  For  example,  the  flad,  hi  is  the  address  of  the  first 
instruction  in  the  program.  After  all  the  information  has  been  fed  into 
the  computer  the  following  instruction  is  used,  "START  AT  hi"  and  the 
program  begins* 

After  the  flad  assignments  had  been  made,  the  next  step  was  to 
draw  a  flow  diagram  such  as  shown  in  Figure  III.  The  purpose  of  this 
diagram  was  to  aid  in  the  organization  of  the  program  and  to  provide  an  orderly 
procedure  for  the  programmer  to  follow.  It  was  a  picture  of  the  problem 
showing  the  sequence  of  operations  to  be  performed.  The  programmer  then 
wrote  a  short  routine  that  would  do^  the  calculations  of  each  section  shown 
in  the  flow  chart.  The  last  step,  then,  was  to  combine  these  separate 
routines  into  one  program  that  would  solve  the  entire  problem.  Of  course 
considerable  care  must  be  exercised  in  writing  and  combining  these  routines 
since  the  results  of  one  section  will  influence  the  calculations  in  the 
following  sections. 

When  the  values  of  frequency  and  residual  end  moment  had  been 
generated  in  the  arithmetic  element,  they  were  transferred  to  a  magnetic 
tape  via  an  output  procedure.  When  the  program  had  been  completed  these 
results  were  then  transferred  to  a  flexowriter  and  typed.  The  results  were 
presented  in  two  columns  of  frequency  vs.  end  moment.  Thus,  it  was  possible 
to  determine  the  frequency  at  which  the  moment  becomes  zero  by  either 
linear  interpolation  or  a  plot  of  the  data. 


Figure    HL         FLOW     CHART 
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STRESS  CALCULATIONS  FOR  A  NON-SYMMETRICALLY  LOADED  WEB  FRAME 
Description  of  the  Problem 

A  transverse  web  frame  is  a  part  of  a  system  of  frames  and 
bulkheads  that  maintains  the  shape  of  the  hull.  A  sketch  of  a  web 
frame  is  given  in  Figure  IVa.  To  do  this  it  must  provide  support  for 
longitudinal  strength  members  and  carry  vertical  shear  loads  from  the 
decks  and  hull  bottom  into  the  sides  which  are  the  main  vertical  shear 
carrying  members  of  the  structure.  The  web  frame  is  a  continuous  beam 
ring  that  encircles  the  hull  girth  in  a  transverse  plane.  All  loads 
applied  to  the  frame  are  assumed  to  lie  within  this  plane,  but  may  have 
any  orientation  within  that  plane.  Under  such  loading  a  web  frame  must 
resist  axial  loads,  shear  loads  and  bending  loads.  Since  web  frames  are 
generally  very  long  with  respect  to  depth,  bending  moments  are  by  far  the 
most  important  aspect  to  be  considered  in  the  solution. 

Due  to  the  complex  interaction  between  the  webs,  longitudinals, 
and  bulkheads,  the  exact  solution  of  the  problem  of  web  frames  and  the 
system  of  members  that  maintain  the  hull  shape  has  defied  rational  engin- 
eering theory  up  to  this  time.  Time  limitations  and  the  intent  of  this 
thesis  prohibited  an  investigation  of  possible  improvements  in  the  method 

of  the  transverse  strength  calculations.  It  was  decided  to  adapt  to  a 
computer  solution  a  proven  method.  Although  transverse  strength  calcu- 
lations are  not  always  performed  for  modern  ships,  the  analysis  normally 

f? ) (6) 
utilized,  if  performed,  is  a  long,  laborious  method  developed  by  J.  Bruhn      , 
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Figure  IVa  Transverse  Web  Frame,  Starboard  Side 


Figure  IVb  Web  Frame  Showing  Redundant  Forces 
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Bruhn's  method  is  an  adaptation  of  Castigliano's  theorem  of 
least  work  which  states  that  "The  redundant  reaction  components  of  a 
statically  inoeterminant  structure  are  such  as  to  make  the  internal 
work  a  minimum."      Determination  of  each  of  the  redundant  components 
requires  that  the  differential  of  the  total  internal  energy  with  respect 
to  that  redundant  be  equal  to  zero  which  is  a  minimum  energy  condition. 
Bruhn  neglected  the  axial  and  shear  components  of  energy  and  considered 
only  the  larger  bending  effects.  Bruhn's  basic  energy  equations  for  a 
itfeb  frame  with  two  stanchions,  as  shown  in  Figure  IVb,  then  become 
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Castigliano's  first  theorem  states  that,  "the  deflection  in  a 
certain  direction  at  a  certain  point  of  a  statically  determinate  structure 
is  equal  to  the  partial  derivative  of  the  total  external  work  or  the  total 
internal  energy  with  respect  to  a  load  applied  at  this  point  in  the 
direction  of  the  deflection,"      Conversion  of  a  statically  indeterminate 
structure  into  a  statically  determinate  structure  by  means  of  cutting 
the  structure  and  applying  redundant  forces  as  loads  allows  the  application 
of  Castigliano's  first  theorem.  The  above  equations  then  represent  the 
deflection  between  the  limits  of  the  integration.  Since  the  web  frame  is 
continuous,  the  deflections  and  rotations  at  the  point  of  action  of  the 
redundants  nu,  pn  and  q_  are  equal  to  zero.  For  the  redundant  column 
forces  R  and  S,  the  columns  are  considered  imcompressible  and  the  vertical 
deflection  between  the  column  ends  is  equal  to  zero.  This  idea  of 
deflections  and  rotations  is  used  in  the  description  and  analysis  of  the 
problem  in  Appendix  B# 

To  determine  the  redundants,  five  simultaneous  equations  must 
be  formed  and  solved.  In  a  problem  that  contains  more  or  fewer  redundant- 
supports  such  as  stanchions  or  decks,  the  number  of  redundant  equations 
required  to  solve  the  web  frame  would  of  necessity  be  changed  to  form 
one  equation  for  each  redundant  force. 

After  the  equations  had  been  developed  and  formed  into  a  matrix, 
Grout's  method  of  matrix  reduction  was  used  to  solve  for  the  value  of  the 
redundantSo  These  redundants  were  then  combined  with  the  applied  loads 
and  the  moments  recomputed.  Utilizing  the  computed  moments,  the  stresses 
due  to  bending  were  calculated.  The  direct  forces  were  converted  into  axial 
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and  shear  stresses.  These  stresses  and  the  bending  stresses  were  then 
added  algebraically  to  give  the  total  stresses  in  the  flange  of  the  web 
frame  and  shell  of  the  ship. 

General  Discussion  of  Program 

The  flow  diagram  contained  in  Figure  V  describes  the  sequence 
of  calculations  used  in  calculating  the  frame  stresses.  A  detailed 
description  of  the  instructions  utilised  in  each  sub-routine  and  the 
sequence  control  program  are  contained  in  Appendix  B. 

Due  to  the  normal  manner  of  construction  of  ships,  the  web  frames 
are  constructed  as  a  symmetrical  ring  frame  of  variable  cross-sections. 
The  loads,  however,  are  very  complex  and  for  a  general  solution  must  be 
considered  as  non- symmetrical.  Note  that  for  this  problem  only  the  simple 
static  loads  are  considered. 

The  contours  of  the  ship's  hull  must  be  represented  within  the 
computer  and,  in  order  to  find  the  neutral  axis  and  to  form  the  components 
of  the  loads,  the  slope  of  the  contour  must  be  computed  at  each  station. 
The  most  obvious  solution  for  the  novice  programmer  is  to  form  the  equation 
of  the  contour  and  differentiate  that  equation  to  determine  the  slope.  As 
novices  this  solution  was  attempted,  but  was  soon  abandoned  due  to  the 
complexity  of  the  program  necessary  to  form  the  equation  and  the  inherent 
error  involved  in  first  using  approximate  methods  to  form  the  equation 
and  more  approximate  methods  required  to  utilize  the  equation.  The  method 
finally  selected  was  to  choose  a  series  of  stations  around  the  girth  of 
the  web  and  assume  the  web  is  represented  by  a  series  of  short  straight 
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Read  in  Grout's  method  sub-routine  on  drum 
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Figure  Va  Sequence  of  Operations  for  Web  Frame  Program 
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Figure  Vb  Comparison  and  Alteration  Routine 
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beam  segments.  Stations  were  selected  in  such  a  manner  that  the  angle 
change  between  segments  was  both  small  and  nearly  equal  which,  due  to 
the  variable  curvature,  makes  each  segment  unequal  in  length.  A  station 
was  necessary  at  the  point  of  action  of  each  redundant.  Such  a  proced- 
ure follows  a  common  philosophy  of  digital  computers  in  that  where 
complexity  must  be  weighed  against  a  greater  number  of  calculations 
and  simplicity,  the  choice  must  be  made  in  favor  of  simplicity.  The 
number  cf  stations  required  to  properly  approximate  the  DD  hull  was 
1|9  which  should  be  sufficient  to  approximate  the  frame  of  any  vessel. 

Sin  9  and  cos  9  were  determined  directly  from  the  coordinates 
by  using  the  Pythagorean  Theorem  and  basic  trigonometric  relationships. 

As  =l/dx  +  dy    sin  9  =  J       cos  9  =  .  . 
V      ^  As  As 

This  was  done  for  two  reasons,  first,  only  the  sin  and  cos  are  used  in 
the  solution  and  secondly,  the  tan  9  =  -r — becomes  infinity  when  Ax  =  0 
which  is  the  case  on  ships  with  wall  sides  or  tumble  home.  The  computer 
will  not  handle  division  by  zero  and  the  operations  will  be  stopped  when 
such  a  caSv  occurs.   To  approximate  the  functions  of  the  tangent  angle  at 
each  station,  the  functions  of  the  chord  angles  were  averaged.  While 
sin  9n  +  sin  9/2  does  not  equal  sin(9  +  9  )/2,  the  error  will  be  small 
as  long  as  the  change  in  9  is  small  which  is  another  example  of  using  a 
simple  idea  many  times  to  obtain  accuracy  and  generality. 

Using  the  total  water  head  which  in  the  present  problem  was 
the  height  of  the  main  deck  and  the  average  of  the  y  coordinates  of  the 
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ends  of  the  segments,  the  water  pressure  loads  were  detennined  and  vertical 
and  horizontal  load  components  formed. 

pw  =  H  "  ™2~ 

f  =  p  ds 

w    w 

V  =  f  sin  0 
w    w 

H  =  f  cos  9. 
w    w 

Using  the  inputs  of  the  web  dimensions,  the  distance  of  the 
neutral  axis  (cLT  .  )  from  the  hull  shell  was  calculated  and  the  moment  of 
inertia  about  the  neutral  axis  determined  (!L  .  ) 

mt 

dN.A.  =  At 

"Sj.A.  =  ^hell  "  AtdN.A. 

Knowing  cL.  .   and  the  (sin  ft)   and  (cos  ft)   the  offsets  of  the 
neutral  axis  from  the  shell  could  be  determined, 

XN.A.  =  VdN.A.(sijl  & 

yN.A.  =7s"dN.A.(c0S  ^h 

Using  the  same  procedure  developed  for  the  shell,  the  segment  lengths 
and  functions  of  the  angle  were  formed. 
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One  difficult  engineering  decision  was  the  consideration  of  the 
sharp  corner  at  the  junction  of  the  side  shell  and  the  deck  edge.  In 
normal  practice  the  corner  is  bracketed  and  represents  a  point  of  complex 
discontinuity  x^rhich  nearly  defies  analytical  analysis.  After  a  number 
of  methods  were  considered,  it  was  decided  to  make  the  corner  a  "hard 
spot"  by  including  the  full  bracket  depth.  The  station  spacing  near  the 
corner  was  reduced  to  minimize  the  effect  of  the  corner.  In  all  cases 
the  tangent  at  the  corner  was  taken  as  -1  which  is  equivalent  to  a  slope 
angle  of  135  degrees  on  the  starboard  side.  This  was  done  to  form  a 
smooth  transition  of  the  neutral  axis  around  the  corner. 

The  moment  due  to  applied  loads  and  water  loads  was  determined 

by  m  =  H  dy  +  Vdx,   The  moments  due  to  the  redundant  forces  were: 
n 

V<jT  +  %x   +  R(x"x  )  +  s(x-x,  ). 

The  integrals  of  the  moments  were  formed  by  the  trapezoidal 
rule,  While  the  trapezoidal  rule  implies  linear  segments  between  moment 
ordinates,  this  method  was  necessary  to  accomodate  the  random  station 
spacing.  The  error  resulting  from  the  trapezoidal  rule  will  be  small  as 
long  as  the  station  spacing  is  close  and  the  moment  curve  is  smooth.  Both 
conditions  are  present  in  this  problem. 

After  solving  for  the  redundants  and  re  correal  ting  the  moments  and 
axial  loads  to  include  all  forces,  the  direct  stresses  and  bending  stresses 
were  combined.  The  program  as  written  rives  the  alternative  of  performing 
a  check  calculation  on  the  set  of  frame  scantlings  or  of  altering  the  frame 
scantlings  to  approach  a  more  efficient  design.  The  stresses  were  compared 
with  allowable  stresses  and  the  section  altered  by  changing  the  web  flange 
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Figure  VI  Arrangement  of  Matrix  Formed  by  the  Integration  Sub-routine 
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width  by  predetermined  incremental  amounts.  If  (0"nn   -  (*„)   >  0,  then  the 
J    -  allow   i     * 

flange  was  reduced.  If  (<TT_  -&D   <  0  then  the  flange  was  increased. 

allow  f 

The  outpu.ts  of  each  cycle  change  may  be  any  computed  value  the 
programmer  desires,  but  in  the  present  problem  the  outputs  were  flange 
width,  stress  at  the  flange  and  at  the  shell,  direct  stress  and  shear 
stress  at  each  point.  Other  methods  of  scantling  change  could  have  been 
utilized  or  the  scantlings  could  have  been  altered  by  a  complete  change 
of  scantlings  around  the  girth.  In  the  latter  case  each  cycle  would  have 
required  a  decision  by  the  programmer,  a  procedure  that  has  considerable 
merit  where  design  experience  in  the  mind  of  the  programmer  can  be  better 
utilized.  In  the  present  problem  it  was  decided  to  use  the  conditional 
response  as  another  example  of  the  digital  computer's  capabilities. 

In  the  web  frame  problem  the  quantity  of  generated  numbers 
necessary  to  produce  a  solution  was  so  large  that  the  auxiliary  drum  storage 
had  to  be  utilized.  Blocks  of  numbers  were  moved  from  the  core  memory  onto 
the  auxiliary  drum  to  make  room  for  calculations.  Further,  the  size  of 
the  program  was  so  large  that  optimum  data  output  procedures  could  not  be 
utilized  without  unwarranted  complication  of  the  program.  The  output  data 
could  have  been  typed  or  presented  on  films  in  columns  of  easily  discernable 
numbers.  In  the  present  case  the  size  of  the  program  required  output  via 
an  "fp"  tape  which  gives  five  columns  of  numbers  running  consecutively 
across  the  lines  and  the  address  of  the  registers  containing  the  numbers 
listed  in  the  first  column.  Judiciously  placed  zeros  provide  separation 
between  numbers  of  different  meaning. 
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To  make  the  program  as  general  as  possible  sub-routines  were 
utilized  to  perform  the  calculations,  A  full  description  of  the  use  of  a 
sub-routine  is  contained  in  Appendix  A,  The  sub-routines  used  will  perform 
the  calculation  of  stresses  for  a  web  frame  of  any  dimension  and  shape 
desired.  However,  the  number  of  redundants  of  the  frame  cannot  exceed 
five  and  the  water  pressure  head  must  equal  the  height  of  the  deck  edge. 
Due  to  limitations  of  storage  space,  the  s  equence  instructions  used  will 
permit  the  use  of  only  Ij.9  stations  around  the  girth.  All  of  the  above 
limitations  can  be  exceeded  by  minor  changes  in  the  program  should  the 
necessity  arise. 
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RESULTS 

The  numerical  results  obtained  from  the  computer  solution  of 
problem  one  are  contained  in  Table  I.  Numerical  results  were  not  obtained 
for  problem  two. 

More  important  than  the  numerical  results  are  the  programs 
themselves.  Problem  one  has  been  run  and  is  known  to  be  correct.  The 
sub-routines  of  problem  two,  with  the  exception  of  the  sub-routine  for 
forming  the  stresses,  were  run  and  found  to  check.  When  all  of  the  sub- 
routines were  run  in  conjunction  with  the  sequence  control  program,  several 
minor  programming  errors  were  found. 

The  first  program  will  determine  the  critical  whirling  frequency 
of  any  propeller  shaft  where  no  more  than  four  spans  need  be  considered. 
It  can  be  adapted  to  another  shaft  merely  by  changing  the  physical  con- 
stants of  the  shaft  involved. 

The  second  program  will,  if  the  above  mentioned  errors  are 
removed,  determine  bending,  axial  and  shear  stresses  in  a  transverse 
web  frame  and,  if  desired,  will  adjust  the  size  of  the  web  frame  until 
the  computed  stresses  approach  the  allowable  stresses.  The  program  is 
general  within  the  limitations  stated  below;  the  required  inputs  being 
the  dimensions  of  the  web  frame,  the  location  of  the  stanchions  and 
coordinates  or  points  on  the  transverse  section,  the  weight  and  shear 
loads  at  these  points  and  the  hydrostatic  design  head  above  the  keel. 
The  following  are  the  limitations  of  the  program: 


(1)  the  frame  scantlings  are  symmetrical  about  the  ship  centerline; 

(2)  the  points  defining  the  transverse  section  must  be  symmetrical 
about  the  ship  centerline  and  only  25>  may  be  taken  on  one 
side  ; 

(3)  one  of  these  points  must  be  on  the  centerline  keel,  one 
on  the  deck  centerline  and  one  at  each  end  of  all  the 
stanchions; 

(k)  the  hydrostatic  pressure  head  must  be  equal  to  or  greater 
than  the  height  of  the  deck  edge  and  hydrostatic  pressure 
loads  are  not  generated  for  the  deck; 

(5)  ship  must  be  in  upright  position; 

(6)  only  one  centerline  or  two  off-center  vertical  stanchions 
may  be  used. 


Computer  time  for  calculation  of  problem  one  was  10  minutes, 
but  the  typing  time  of  the  results  was  20  minutes.  Computer  time  for 
the  solution  of  problem  two  was  10  minutes. 
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During  the  period  that  the  thesis  was  being  printed,  the  pro- 
gramming errors  in  the  second  program  were  remedied  and  the  results  listed 
in  Table  II  were  obtained.  Table  II  contains  the  values  of  the  redundants 
£.nd  the  stresses  at  the  keel  and  the  deck  at  centerlinee  Due  to  the 
difference  in  station  spacing  between  hand  and  machine  calculations,  the 
centerline  points  are  the  only  coincident  stations  for  the  two  methods 
of  calculation,, 
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Results  of  Web  Frame  Stress  Calculations 
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As  shown  in  Table  II  there  is  little  correlation  between 
hand  and  machine  calculated  values.  However,  this  lack  of  comparison  is 
caused  by  two  conditions.   The  most  important  condition  is  that  the 
hand  calculations  have  been  found  to  be  of  doubtful  validity.   Several  hand 
calculation  errors  were  found  and  the  scantlings  used  are  not  defined. 
Recent  information  indicates  that  the  scantling  values  used  in  the 
Bureau  of  Ships  hand  calculations  do  not  correspond  with  the  scantling 
values  presented  in  the  contract  plans  used  as  a  basis  for  machine 
calculations o 

The  other  condition  that  causes  a  small  error  is  the  fact 
that  summation  of  vertical  forces  does  not  equal  zero.   The  difference 
between  weight  forces  and  buoyant  forces  must  be  balanced  by  vertical 
shear  forces  applied  along  the  vertical  sides  of  the  hull*  This  shear 
force  was  taken  from  the  hand  calculations  and  was  included  in  the 
input  values  for  the  program.  The  machine  computed  buoyant  forces  were 
different  from  hand  computed  values  due  largely  to  the  difference  in 
station  spacing  and  the  method  of  averaging  the  water  pressure.  Therefore, 
the  input  values  of  shear  did  not  balance  vertical  forces  computed  by 
the  machine.  It  is  now  obvious  that  the  shear  force  values  must  be 
computed  by  a  sub-routine  in  the  computer  program,, 

This  lack  of  the  means  of  proper  comparison  points  out  one 
other  general  rule  of  programming.  The  programmer  should  perform  hand 
calculations  for  the  problem  being  programmed  using  the  method  adapted 
to  the  computer o   The  authors  decided  that  the  use  of  existing  hand 
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computations  would  be  a  time  saving  device j  however,  as  the  results 
indicate,  this  was  a  poor  decision  and  it  is  hoped  that  future  prop-rammers 
will  benefit  by  the  author's  experience. 

With  the  exception  of  the  shear  error  discussed  above,  the 
program  does  work  and  produces  good  results.  With  slight  correction  the 
web  frame  program  is  suited  for  research  and  engineering  calculations. 

The  time  required  for  machine  calculations  was  25  seconds; 
however,  program  conversion  required  five  minutes  and  typing  time  for 
results  of  calculations  required  twenty  minutes., 

Time  limitations  did  not  permit  the  running  of  the  comparison 
and  alteration  program  for  the  web  frame. 
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DISCUSSION  OF  RESULTS 

The  subject  of  computer  application  to  naval  architectural 
problems  is  as  general  as  the  field  of  naval  architecture  itself.  Each 
proposed  application  will  require  a  separate  analysis  of  feasibility. 
However,  certain  general  conclusions  and  impressions  can  be  gained  from 
the  two  problem  examples  of  this  thesis.  To  each  general  conclusion  there 
exists  certain  possible  exceptions;  however,  these  exceptions  orHj   tend 
to  prove  the  contention  that  digital  computers  are  not  universal  panaceas, 
but  only  local  aids. 

The  programming  errors  found  in  problem  two  are  normal  errors 
that  occur  in  every  prograjn  and  do  not  reflect  upon  the  method  of  approach. 
The  nature  of  these  errors  are  incorrect  addresses,  typographical  mistakes 
and  human  errors.  The  correction  of  these  difficulties  is  only  a  matter 
of  time  and  patience. 

The  results  of  the  two  programs  indicate  that  digital  computers 
are  applicable  to  long  and  complicated  naval  architecture  problems  and 
that  the  numerical  results  are  more  than  comparable  with  those  of  hand 
calculations.  Of  greater  engineering  importance  than  accuracy,  was  the 
great  reduction  jn  calculation  time  obtained  by  computer  usage.  For  the 
first  problem  the  computer  required  10  minutes  as  opposed  to  an  estimated 
hand  calculation  time  of  thirty  (30)  hours  and  in  the  second  problem  a 
more  defined  superiority  is  estimated  in  that  the  comouter  solution  will 
require  approximately  ten  (10)  minutes  versus  an  estimated  sixty  (60) 
hours  by  hand.  However,  this  comparison  does  not  depict  the  true  time 
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picture  and  the  long  tedious  hours  of  program  writing  and  ever  present 
"debugging"  or  program  correction  must  be  considered.  Although  the  pro- 
gramming time  varies  greatly  with  the  programmer's  familiarity  with  the 
problem,  computer  and  also,  with  the  nature  of  the  problem  itself •   in 
general  the  programming  time  and  the  time  required  for  one  hand  calculation 
are  roughly  commensurate.  The  authors  estimate  that  a  practicing  naval 
architect  very  familiar  with  the  proposed  problem  and  experienced  in 
programming  could  program  problem  one  in  twenty  (20)  hours  and  problem 
two  in  forty  (1|0)  hours. 

The  "debugging"  time  for  each  problem  would  require  another 
ten  (lO)  hours.   The  fact  that  the  authors,  unfamiliar  with  both  the 
problem  and  computer,  required  six  times  that  amount  of  time  indicates  that 
experienced  personnel  are  a  very  important  application  factor.  In  the 
present  case  problem  understanding  and  organization  were  responsible  for 
the  greatest  loss  of  time.  As  a  general  statement  it  can  be  said  that  for 
an  engineer  to  properly  program  a  particular  problem,  his  knowledge  of 
that  problem  must  be  more  thorough  than  that  required  for  hand  calculations. 
Not  only  must  he  know  the  individual  steps  in  the  problem,  but  also  must 
predict  the  effect  of  each  step  upon  the  following  step  and  the  general 
value  of  the  numbers  generated. 

Also,  the  authors  found  that  one  programmer  working  alone  is 
just  as  effective  as  several  working  together  on  the  same  problem.  The 
reason  for  this  is  the  necessity  for  a  program  to  follow  a  single  trend 
of  thought  and  the  exact  accounting  required. 
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To  provide  a  more  complete  time  description  of  orosrraimiing, 
the  nature  of  the  proposed  problem  should  be  analyzed.  Problem  solutions 
that  involve  repetition  of  the  same  equations  will  tend  to  have  a  much 
shorter  prograjmriing  time.  For  instance  the  time  required  to  program 
problem  one  for  i|0  spans  of  shafting  would  be  only  slightly  longer  than 
that  required  to  program  h   spans.  Also,  the  program  required  to  design 
the  web  frame  by  cycling  and  altering  the  scantlings  is  little  longer 
than  the  program  for  check  calculations. 

After  the  solution  of  one  particular  problem  by  a  computer  a 
general  program  will  find  further  use  in  solving  other  similar  problems. 
The  programs  written  for  this  thesis  are  reasonably  general  and  the 
limitations  stated  in  the  results  are  not  of  serious  damage  to  the  program, 
1    The  limitation  of  problem  one  to  four  shaft  spans  in  a  reasonable  engin- 
eering division  in  that  propeller  excited  whirling  vibration  amplitudes 
. beyond  span  four  are  so  small  that  their  effect  on  the  end  fixity  of 
span  foui  is  insignificant.  In  problem  two  the  limitations  were  imposed 
to  save  storage  space  and  reduce  the  complexity  of  the  program.  This 
coincides  with  the  authors'  opinions  that  for  engineering  applications 
several  simple  programs  are  easier  and  faster  to  handle  and  program  than 
one  complex  all  purpose  program.  For  instance,  web  frames  are  normally 
constructed  symmetrically,  25  points  should  properly  define  most  normal 
web  shapes  and  to  have  a  station  at  each  major  load  or  redundant  is 
reasonable.  The  design  hydrostatic  head  is  normally  above  the  maindeck 
and  waterloads  are  normally  considered  a  constant  value  along  the  deck. 
While  the  listed  position  of  the  ship  is  an  important  consideration, 
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modern  naval  architects  consider  the  upright  position  for  practical 
calculations  and  should  the  listed  condition  become  necessary  for  re- 
search calculations,  rotation  formulas  and  non- symmetrical  water  loadings 
must  be  incorporated  into  the  program.  To  change  stanchions  and  add 
decks  as  redundant  supports  major  alteration  of  integration  procedures 
would  be  required. 

A  general  program  is  useful  in  design  and  research  to  compare 
and  generate  great  quantities  of  data  and  to  test  the  results  of  slight 
alterations  upon  a  major  problem.  For  instance,  program  one  could  analyze 
and  compare  the  whirling  frequencies  of  hundreds  of  ships  in  a  relatively 
short  time  or  the  effects  of  propeller  weight,  fixity  and  blade  number 
may  be  determined  for  numerous  combinations  of  values.   Problem  two 
would  find  use  in  generating  design  schantlings  for  various  load  and  size 
variations  such  as  those  found  in  classification  society  publications. 
As  a  sub-routine  in  conjunction  with  a  more  complex  program  the  stiffness 
of  web  frames  relative  to  transverse  bulkheads  could  be  determined  by 
application  of  the  theory  of  beams  on  elastic  foundations.  By  reiteration 
of  calculations  under  varying  loads  the  variant  or  impactive  loadings 
may  be  analyzed.  Since  the  longest  time  consumed  in  a  program  is  in  the 
program  writing,  problem  repetition  increases  the  efficiency  of  the  com- 
puter application. 

To  be  feasible  for  engineering  calculations  the  computer  must 
save  time,  money  and  increase  the  effectiveness  of  the  engineers  it  aids. 
Therefore,  the  computer  must  accept  numbers  of  engineering  magnitude  and 
quantity  in  a  convenient  form.  The  coded  instructions  must  be  straightforward 
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and  have  names  descriptive  of  their  use.  Also,  the  generated  results  must 
be  presented  in  such  a  manner  that  wide  spread  publication  is  possible. 
Whirlwind  I  and  some  commercially  available  machines  will  accept  decimal 
digits  of  sufficient  length  and  quantity  to  perform  most  practical 
engineering  and  research  problems.  Of  extreme  importance  is  the  necessity 
of  maintaining  a  large  library  of  sub-routines  capable  of  performing  such 
calculations  as  sin  9,  l/x7  etc.  Frequently  it  is  desirable  to  purchase 
the  services  of  an  experienced  programmer  to  act  as  a  consultant  when 
first  using  a  computer. 

In  computer  application  there  exists  two  basic  schools  of 
thought.  One  group  thinks  that  problems  should  be  programmed  by  experts 
who  perform  no  other  function  than  programming.  In  this  case  the  engineer 
presents  the  problem  to  the  programmer  and  at  some  later  time  is  riven  a 
solution.  The  other  group  states  that  a  computer  is  best  utilized  by 
having  engineers  familiar  with  the  problem  effect  the  programming.  While 
there  exists  a  middle  road  somewhere  between  the  two  extremes,  the  author- 
are  of  the  opinion  that  engineers  should  program  the  computer.  This 
performs  two  useful  functions  for  the  profession  as  a  whole.  One,  the 
programmer  can  make  the  necessary  engineering  decisions  to  adapt  the 
problem  more  successfully  and  secondly,  the  engineer  becomes  more  aware 
of  possible  computer  applications  in  his  field. 

The  qualities  required  of  a  programmer  are  no  more  exacting 
than  those  required  of  a  good  naval  architect.  He  must  have  almost 
infinite  patience,  imagination,  ingenuity,  good  ability  to  organize 
problems  and  a  good  background  in  both  cheoretical  and  numerical  methods 
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of  mathematics.  The  work  is  tedious  and  all  instructions  must  be  perfect 
since  the  machine  will  generate  wrong  answers  as  fast  as  correct  ones. 
Digital  computers  are  often  referred  to  as  "brains,"  but  nothing  could  be 
more  incorrect  in  that  thought  is  impossible.  Actually  a  computer  more 
nearly  resembles  a  desk  calculator  or  an  abacus  which  operates  at  a  high 
rate  of  speed.  Imagination  and  ingenuity  are  so  necessary  to  produce  a 
workable,  general  program  that  many  commercial  computing  companies  obtain 
their  programmers  from  chess  clubs,  cross-word  puzzle  enthusiasts  and 
other  hobbies  where  deep  concentration  and  original  thought  are  required. 

The  future  of  digital  computers  in  naval  architecture  appears 
very  promising  from  both  the  research  and  practical  viewpoints.  In 
research  they  can  be  used  for  the  rapid  assiinilation  and  processing  of 
large  quantities  of  technical  data  and  for  the  solution  of  problems  that 
require  an  inordinate  amount  of  time  to  calculate  by  hand.  Typical  examples 
of  the  d  ata  processing  type  of  problem  would  be  a  comparison  of  stability 
among  several  different  vessels  and  other  similar  problems  involving 
statistical  analysis  to  establish  design  criteria.  Since  the  same  pro- 
gram may  be  used  to  solve  a  number  of  similar  problems,  design  data  in 
the  form  of  charts  and  tables  can  be  generated  for  the  grillage  beam 
analysis  of  decks  and  other  complex  problems.  Such  a  procedure  permits 
the  determination  and  testing  of  important  parameters  for  later  applica- 
tion to  specific  designs. 

For  specific  design  applications  computers  may  be  used  to  effect 
large  savings  in  time  and  labor.  This  is  possible  because  the  computer 
can  approach  more  closely  an  optimum  design  with  little  increase  in 
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computing  time  and  expenditure  of  effort.  In  this  category  of  appli- 
cation are  midship  section  design,  transverse  framing  design  and  flight 
deck  bents.   In  addition,  refinements  intheory  and  in  the  extent  of  calcu- 
lations can  also  be  utilized  to  a  greater  degree  than  practical  by  hand. 
Examples  of  these  refinements  are  the  use  of  grillage  beam  analysis  for 
decks  and  bulkheads,  the  performance  of  stress  calculations  for  a  ship 
in  drydock,  design  of  propellers  by  circulation  theory,  and  the  investi- 
gation of  submarine  dynamical  stability.  In  addition  to  desifm  calcu- 
lations a  digital  computer  can  be  applied  very  readily  to  scheduling 
and  accounting  problems.   In  shipbuilding  this  system  can  be  used  in 
determining  the  weight  and  center  of  gravity  of  a  ship  as  it  is  being 
built.  Another  possibility  is  merchant  ship  cargo  arrangement. 

Theses  present  an  opportunity  for  wide  and  varied  use  of 
digital  computers.  In  this  application  the  computer  should  be  used  as 
an  aid  and  not  as  the  general  subject.  A  student  attempting  a  problem 
in  which  the  magnitude  of  the  calculations  is  the  dominant  factor  should 
consider  computer  aid.  A  desireable  arrangement  is  to  have  a  team  with 
one  member  learning  the  problem  and  another  member  learning  the  computer. 
In  this  manner  both  can  derive  maximum  benefit  from  the  thesis  and  a 
larger  contribution  to  the  profession  can  be  made.  Since  computers  are 
becoming  more  widely  used,  familiarity  with  their  capabilities  may  be 
of  considerable  future  value  to  the  student. 
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CONCLUSIONS 


1.  Digital  computers  are  applicable  to  naval  architecture  problems. 

2.  The  average  naval  architect  has  sufficient  background  to  program 
a  digital  computer. 

3.  Some  general  requirements  of  a  programmer  are: 

a.  They  should  have  a.  complete  engineering  understanding 
of  the  problem  to  be  programmed. 

b.  Familiarity  with  numerical  methods  of  mathematical 
processes  is  necessary. 

c.  They  should  be  amenable  to  tedious  work. 

ho     The  problem  to  be  programmed  must  be  suitable  to  comcuters  since 
it  may  take  longer  to  program  than  to  determine  one  solution  by 
hand. 

5>.  One  programmer  can  work  faster  than  two  or  more. 

6.  In  general  the  limitations  imposed  on  a  problem  a  re  not  serious  and. 
are  no  more  restrictive  than  those  in  hand  calculations, 

7.  Complete  generality  of  a  program  is  not  always  feasible. 

8.  For  engineering  applications  the  following  features  of  computers 
are  necessary: 

a.  Computer  word-length  should  be  able  to  handle  engineering 
numbers. 

b.  Library  of  sub-routines  should  be  available. 

c.  Computer  should  be  easy  to  program. 
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RECOMMENDATIONS 

1.  Digital  computers  should  be  used  to  solve  naval  architectural 
problems • 

2.  Digital  computers  should  be  used  in  theses  and  research  problems. 

3.  Any  future  thesis  should  consider  the  computer  as  an  aid  in  the 
solution  of  a  specific  problem  and  not  applicability  of  computers 
to  a  general  field. 

h»     For  engineering  applications  the  programmer  should  have; 

a.  good  engineering  understanding  of  his  field; 

b.  knowledge  of  numerical  methods; 

c.  organizational  qualities  and  a  great  deal  of  ingenuity; 

d.  a  working  knowledge  of  the  computer  utilized. 
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Supplementary  Introduction 

A  computer  program  is  a  series  of  logical  instructions 
telling  the  computer  what  mathematical  operations  are  to  be  performed 
on  the  data.   Thus  there  are  two  groups  of  inputs,  the  data  and  the 
instructions.  The  data  fed  into  the  computer  is  stored  in  the  core 
memory.  Each  number  or  piece  of  data  is  stored  in  a  register  or  "cell" 
the  location  of  which  is  designated  by  a  number  called  the  address  of 
that  register.  After  all  the  data  is  stored  the  instructions  are 
stored  in  the  unused  portion  of  the  core  memory.  The  instructions  are 
stored  in  sequence  and  will  be  performed  in  the  order  in  which  they 
are  stored.  No  unused  registers  are  permitted  to  interrupt  the  storage 
as  these  will  stop  the  computer.  The  programmer  may  wish  to  know  the 
address  of  key  instructions  such  as  the  one  that  starts  the  computation 
process o  He  may  assign  an  address  to  these  instructions  or  allow  the 
computer  to  do  so,  as  explained  later. 

The  instructions  given  to  the  computer  consist  of  two  parts, 
the  operation  and  address  sections.  The  operation  section  is  first  and 
uses  letters  to  tell  the  computer  what  mathematical  operation  is  to  be 
performed.  The  address  section  gives  the  address  of  the  register  contain- 
ing the  data  on  which  the  operation  is  to  be  performed.  The  instructions 
recognized  by  Whirlwind  I  are  listed  in  the  glossary. 

When  writing  a  program  for  a  digital  computer,  there  are  several, 
useful  devices  that  simplify  the  procedure.  The  programming  devices  described 
below  apply  to  Whirlwind  I  and  may  or  may  not  apply  to  other  computers. 

Probably  the  most  useful  device  is  a  floating  address.  A 
floating  address  or  flad  is  used  to  eliminate  the  necessity  of  referring 
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to  the  absolute  register  address  in  the  core  memory  when  writing  the 
program.  At  the  beginning  of  a  program  the  inputs  or  key  instructions 
on  which  mathematical  operations  are  to  be  performed  are  given  floating 
addresses.  These  addresses  consist  of  a  letter  and  a  number,  for 
example  al,  a2  etc.  Any  letter  except  0  and  L  may  be  used  and  the  sum 
of  all  the  numbers  behind  the  letters  may  not  exceed  255.  In  writing 
the  program,  then,  the  flad  is  used  in  the  address  section  of  the 
instruction  to  represent  the  absolute  address.  The  computer  will  assign 
absolute  addresses  or  register  numbers,  to  these  flads  unless  this  is 
done  by  the  programmer.  For  example,  if  al  is  the  first  flad  in  the  prog- 
ram, the  computer  will  assign  it  to  register  32.  If,  however,  the  prog- 
rammer wants  al  in  the  register  130,  he  writes  130/al  and  the  computer 
will  place  al  in  register  130.  A  table  is  usually  produced  by  the 
computer  showing  absolute  addresses  assigned.  It  should  be  remembered 
that  a  flad  can  have  only  one  value  in  any  one  program. 

Another  useful  device  is  preset  parameters.  These  are  numbers 
other  than  data  that  are  used  in  the  program.  Preset  parameters  are 
frequently  used  to  designate  the  number  of  times  a  computation  is  to  be 
repeated.  They  differ  from  flads  in  that  they  may  have  any  absolute 
value  desired  whereas  flads  must  have  the  absolute  address  of  a  register. 
Preset  parameters  must  be  assigned  a  value  before  they  are  used  and  they 
will  retain  this  value  until  it  is  specifically  changed.  Preset  parameters 
are  designated  by  two  letters  and  a  number.  The  first  letter  must  be 
p,  u  or  z  and  the  second  letter  and  the  number  may  be  anything,  except  the 
letters  o  and  1  cannot  be  used. 
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In  writing  a  program  for  a  digital  computer  subroutines  are 
very  helpful.  Basically  there  are  two  distinct  types  of  subroutines  which, 
for  lack  of  a  better  name,  shall  be  called  internal  and  external.  Internal 
subroutines  are  those  performed  automatically  by  the  computer  upon  receipt 
of  a  single  instruction.  An  example  of  an  internal  subroutine  is  the 
group  of  output  instructions  listed  in  the  glossary.  Each  of  these 
instructions  will  cause  the  computer  to  go  through  a  certain  predeter- 
mined procedure  without  any  further  information  from  the  programmer. 
External  subroutines  are  further  divided  into  two  categories ;  one  category- 
is  designated  library  subroutines  and  are  perfectly  general  for  any  prog- 
ram. The  other  external  subroutines  are  groups  of  instructions  that 
perform  a  repetitive  section  of  a  particular  program.  Library  subroutines 
are  groups  of  instructions  usually  written  by  an  expert  programmer  for 
the  purpose  of  performing  a  particular  mathematical  operation.  They 
may  be  an  integral  part  of  the  main  program  or  they  may  merely  serve  as 
an  auxiliary  program.  Examples  of  these  external  library  subroutines 
are  groups  of  instructions  used  to  determine  the  square  root  of  a  number, 
trigonometric  and  hyperbolic  functions  of  angles  and  the  solutions  to 
simultaneous  equations.  Most  computer  installations  have  libraries  of 
subroutines  available  for  programmers.  These  are  written  quite  arbitrar- 
ily by  using  relative  addresses  and  can  be  used  with  any  main  program. 

There  are  some  rules  that  must  be  followed  when  using  subroutines. 
The  main  program  must  be  interrupted  to  call  the  subroutine  into  use. 
After  the  subroutine  has  performed  its  function,  it  must  direct  the  con- 
trol element  back  to  the  proper  instruction  in  the  main  program.  An  example 
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of  a  subroutine  follows: 

Major  program 

ica  a2 
f6,        isp  g3  —  This  instruction  directs  control  to  follow  the 
instruction  in  the  register  designated  by  g3. 

SUB-ROUTINE 

g3,  ita  g3  plus  20  This  instruction 

places  the  address  contained  in  the 
control  element  in  g3  plus  20  and  adds 
one  to  that  address.  In  this  example 
the  address  contained  in  the  control 
element  is  f6. 

iad  t6 

more  instructions 


g3  plus  20,  isp  f6  plus  1  —  The  address  of 
the  isp  instruction  was  supplied  by  the 
ita  instruction  at  the  beginning  of  the 
subroutine • 
f6  plus  1,  iad  a3 

The  authors  do  not  wish  to  imply  that  a  program  should  consist 
of  a  series  of  subroutines,  but  they  are  very  useful  in  programs  that 
perform  similar  operation?  many  times  throughout  the  problem. 
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The  numbers  usually  encountered  in  engineering  are  too  large 
to  be  stored  in  one  register  which,  it  will  be  recalled,  is  l6-binary 
or  about  k   decimal  digits  long.  To  accommodate  these  larger  numbers 
a  technique  called  floating  point  arithmatic  is  used.  Two  consecutive 
registers,  giving  32  binary  digits  available  for  each  number,  are  used 
to  replace  the  single  register.  Some  of  these  32  digits  are  used  to 
record  the  actual  number  and  some  to  locate  the  decimal  point.  In  this 
manner  the  range  of  significant  digits  can  be  increased  by  restricting 
the  number  of  digits  used  to  locate  the  decimal  point.  In  the  binary 
system  the  decimal  point  is  written  as  some  power  of  2  and  usually  6 
binary  digits  are  used  to  express  this  power.  Since  the  first  digit 
of  each  register  is  used  for  the  sign  of  the  quantity  in  that  register, 
in  this  case  the  sign  of  the  second  register  is  the  sign  of  the  exponent 
of  2.  Twenty-four  digits  are  left  to  express  the  significant  number.  This 
particular  distribution  of  binary  digits  is  referred  to  as  the  (2I4.6)  sys- 
tem and  can  handle  numbers  equivalent  to  eight  decimal  digits  with  expon- 
ents  up  to  10  .  Other  systems  are  available  and  the  size  of  the  numbers 
in  the  problem  determines  which  one  should  be  used. 

In  order  to  manipulate  these  two  register  numbers,  an  inter- 
pretive mode  has  been  internally  programmed  to  function  as  a  part  of 
Whirlwind  I.  All  instructions  that  refer  to  double  length  registers 
are  preceded  by  an  i.   If  calculations  or  operations  are  to  be  performed 
on  single  length  registers,  the  same  instructions  may  be  used  without  the 
i  preceding.  Hence  there  are  two  modes  of  calculation]  interpretive 
for  double  register  length  numbers  and  whirlwind  for  single  length  register 
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numbers.  It  is  possible  to  pass  from  the  whirlwind  mode  to  the  inter- 
pretive mode  by  use  of  the  instruction  IN  and  to  reverse  the  direction 
by  use  of  the  instruction  OUT. 
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GLOSSARY  OF  INSTRUCTIONS 


TERM 


DEFINITION 


ica  x 


ics  x 


iad  x 


isu  x 


unr  x 


idv  x 


its  x 


lex  x 


isp  x 


lcp  x 


clear  MRA  and  add  contents  of 
register  x. 

clear  MRA  and  subtract  contents 
of  register  x. 

add  contents  of  register  x  to 
what  is  in  MRA. 

subtract  contents  of  register  x 
from  what  is  in  MRA. 

multiply  what  is  in  MRA  by  the 
contents  of  register  x  and  round 
off  product  to  fifteen  digits. 

divide  contents  of  MRA  by  the 
contents  of  register  x. 

transfer  the  contents  of  the  MRA 
to  register  x  losing  what  was 
previously  in  register  x.  The 
MRA  remains  unchanged. 

exchange  the  contents  of  the  MRA 
with  the  contents  of  register  x. 

take  the  next  instruction  from 
register  x.  Does  not  affect  the 
contents  of  the  MRA©  An  uncon- 
ditional instruction  used  to  break 
the  sequence  of  operations. 

if  the  contents  of  the  MRA  are 
negative,  take  the  next  instruc- 
tion from  register  x  and  continue 
from  there.  If  the  contents  of 
the  MRA  are  positive,  take  the 
next  instruction  in  sequence.  A 
conditional  instruction  used  to 
break  the  sequence  of  operations. 
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Glossary  of  Instructions  (Continued) 

TERM  DEFINITION 

isc  x  select  counter  number  x.  Without 

"this  instruction  counter  zero 
will  be  used  whenever  +C  appears. 

icr  x  cycle  reset.  Sets  index  register 

of  counter  to  zero  and  criterion 
register  to  x. 

ict  x  cycle  transfer.  Increases  con- 

tents of  index  register  by  one. 
If  contents  of  index  register 
greater  or  equal  to  contents  of 
criterion  register,  set  index 
register  to  zero  and  do  next 
instruction  in  sequence.  If  con- 
tents of  index  register  less  than 
contents  of  criterion  register, 
take  next  instruction  from  register  x. 

ici  x  increase  the  contents  of  index 

register  by  number  x. 

icd  x  decrease  the  contents  of  index 

register  by  number  x. 

iti  x  transfer  the  right  eleven  digits 

of  the  index  register  into  the 
right  eleven  digits  of  register  x. 

iat  x  add  contents  of  index  register  to 

the  contents  of  register  x  and 
store  the  result  in  the  index 
register  and  register  x. 

ita  x  replace  the  address  section  of  the 

instruction  in  register  x  with  the 
address  that  is  one  more  than  the 
address  of  the  register  containing 
the  last  isp  or  icp  if  the  contents 
of  the  MRA  are  negative. 
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Glossary  of  Instructions  (Continued) 

TERM  DEFINITION 

iTOA  record  the  contents  of  MRA  on 

direct  printer  (typewriter). 

iMOA  record  the  contents  of  MRA  on 

magnetic  tape  for  delayed  printing, 

iSOA  record  the  contents  of  MRA  on 

oscilloscope  for  photographing. 

iFOR  this  instruction  provides  an 

automatic  device  for  obtaining  a 
suitable  layout  of  output  data 
in  columns,  lines  or  blocks. 

Note:  The  above  output  instructions  are  usually  followed  by  a  series 
of  letters  and  numbers  that  indicate  the  desired  form  and 
arrangement  of  the  output. 
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Details  of  Program  for  Determining  the  Whirling 
Frequencies  of  Propeller  Shafts 

The  first  step  in  programming  this  problem  is  to  put  the 
numerical  values  of  various  quantities  into  the  computer*  This  is  done 
by  the  use  of  floating  addresses,  or  "flads."  The  proper  use  and  func- 
tion of  flads  is  explained  in  Appendix  A. 

The  following  table  lists  the  flads,  the  quantity  and  the 
values  used.  Figure H  defines  the  symbols  and  nomenclature. 

al,       I»n  =  77"  -  length  of  tail  section. 
a2,       1^  -   5l3«6,r  -  length  of  inboard  sections. 
L2  =  576.0" 

L3  =  51Uu6'« 

10      2 
bl,     ^r^-r\   *  17  «U  *  10   lb-in  -  E  times  moment  of  inertia  of 

tail  section. 

10      2 
b2,     EL  T  =  13.7  X  10   lb-in  -  E  times  moment  of  inertia  of 

inboard  sections. 
E2I2  m  12.1).  X  10xu  lb-in 

10  9 

EI     «  12.lt  x  10       lb-in 
cl,  W     =  ii7,8l5  lb  -  weight  of  propeller 

IT 

c2,  </..    ■  .281|  lb/in     -  specific  weight  of  material  used  in 

inboard  shaft  sections. 


y2  =  .281;  lb/in3 


</    «  .28ii  lb/in3 


2 
dl,      A_^  ■  193.5  in  -  average  area  of  shaft  cross  section 

for  each  span. 

A2  m  17U«5  in2 

A3  «  17li.5  in2 
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d2,  3  -  number  of  shaft  spans  used. 

d3,  135  -  number  of  increments  of  p  desired  to 

be  used.  $   is  increased  by  increments 
of  2  degrees. 

•3,  (|*i)   q  =  6 

q  is  the  order  of  frequency  being 
/l  __  L      /-   investigated.  Plus  sign  is  used 
2   q  for  backward  precession  and  the  minus 

sign  for  forward  precession. 

(2-q}  q=1 

eh9  J  ■  2I1J4..O  *  10  -  polar  moment  of  inertia  of 

propeller  and  entrained  water. 

e6,  P_  =  90  -  initial  starting  value  of  $   ,  in  degrees. 

The  above  listed  flads  are  the  only  ones  that  must  be  changed 
in  order  to  adapt  this  program  to  another  problem.  In  order  to  accomplish 
this  change  it  is  merely  necessary  to  add.  a  tape  with  the  correct  values 
at  the  end  of  the  program. 

The  following  table  lists  flads  that  contain  constants  which 
will  remain  unchanged  for  every  problem: 

el,  tt  =  3.1I1I6 

e2,  180 

e5>,  g  ■  386  in-sec 

e7,  P-,  -  current  value  of  p.      (p     +  increment) 

TT 

e9, 


(180)2 


elO,  2 

ell,  1 

el2,  3 

el3,  h 


$k. 


A  block  of  registers  consisting  of  f  and  g  flads  has  been 
set  aside  as  intermediate  storage ,  This  is  necessary  to  enable  the 
computer  to  calculate  quantities,  store  them  temporarily  and  use  them 
later  in  the  program.  Another  group  of  registers  consisting  of  t  and  u 
flads  are  set  aside  for  temporary  storage.  These  quantities  are  used 
once  and  then  discarded. 

The  last  group  of  flads  used  is  the  h  and  m  group.  These  are 
used  to  control  the  program  by  assigning  key  instructions  to  them. 

The  outputs  from  this  program  will  consist  of  a  series  of 
frequencies  and  corresponding  end  moments.  These  values  will  be  arranged 
in  two  columns  with  the  frequency  corresponding  to  the  moment  that 
appears  on  the  same  line.  It  is  very  important  to  include  the  sign  of 
the  moment,  for  we  are  more  interested  in  the  point  where  the  moment 
passes  through  zero  than  the  actual  values  obtained. 

In  performing  certain  routine  arithmetical  operations,  such 
as  finding  the  square  root,  use  was  made  of  the  sub-routine  library  in 
the  computer  laboratory.  These  sub-routines  are  very  reliable  and  save 
considerable  time  when  programming. 

The  actual  orders  used  for  programming  this  problem  follow. 
The  short  paragraph  at  the  beginning  of  each  section  explains  briefly 
the  function  of  the  orders  in  that  section. 
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These  three  orders  place  the  sub-routines  for  square  root, 
sine-cosine  and  sinh-cosh  at  the  beginning  of  the  program,  This  is 
very  desirable  as  it  eliminates  any  confusion  with  the  main  program,. 

Order 

51,  LSR,  FU8 

52,  LSR,  FUiia 

53,  LSR,   FU5 

Put  in  flad  assignments  al  through  glO  as  previously  noted. 

These  orders  are  assigning  values  of  3,  2,  3  in  the  flads 
HL,  M2,  M3,  respectively,  the  function  of  these  flads  is  to  set  the 
cycle  counter  criterion  registers  to  a  predetermined  value.  The  opera- 
tion is  done  in  WW1  since  the  numbers  are  single  register  numbers. 
The  last  order  is  used  to  enter  the  CS  mode. 

Order       Input  to  MRA      In  MRA      Transfer 

from  MRA 

hi,  ca,  d2  3            3 

td,  ml  3           3 

td,  m3  3           3 

su,  ell  -1           2 

td,  m2  2           2 

IN 

This  group  of  orders  will  calculate  the  weights  of  each  shaft 
span  and  store  the  results  in  fl.  The  first  order  says  that  the  opera- 
tion is  to  be  repeated  three  times.  The  +C  at  the  end  of  the  next 
four  orders  is  to  permit  the  address  section  of  the  instruction  to  be 
increased  by  two  registers  on  each  cycle  thus  allowing  three  different 
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weights  to  be  calculated  by  the  same  set  of  orders.  The  last  instruc- 
tion is  necessary  to  return  control  to  the  first  instruction  at  the  end 
of  each  calculation,  to  increase  the  address  of  each  instruction  oj   two 
registers  after  each  operation,  and  to  pass  control  to  the  next  order 
when  the  operation  has  been  performed  the  number  of  times  specified  by 
the  first  order. 


Order 

Input  to  MRA 

In  MRA 

Transfer 
from  MRA 

ml, 

icr 

h2, 

ica,  c2+c 

-/i       • 

ii 

imr,  a2+c 

*  h 

Vyi 

imr,  dl+c 

X   A 

V  V 

A 
1 

its,  fl+c 

wl 

wi 

This  group  of  orders  will  calculate  C  and  store  the  results 
in  f2.  The  cycle  will  be  performed  three  times  since  the  order  ml,  icr 
is  still  effective. 

Order      Input  to  MRA      In  MRA.      Transfer 


from  MRA 

ica,  b2+c 

\h 

hh 

idv,  a2+c 

+  Li 

Li 

its,  f2+c 

Yi 

h 

C 
a 

ict,  h3 
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This  group  of  orders  will  calculate  C  for  each  span  and 
store  the  results  in  f3.  The  cycle  will  be  repeated  three  times  since 
the  order  ml,  icr  is  still  effective. 


Order 

Input  to  MRA 

In  MRA. 

Transfer 
from  MRA 

hii, 

ica,  f2+c 

idv,  a2+c 

+  h 

idv,  a2+c 

*h 

"2 

idv,  fl+c 

^wi 

its,  f3+c 

cd 

1  1 


ict,  hk 

This  group  of  orders  will  calculate  the  constant  tt/(i8o) 
and  store  the  result  in  e°.  This  operation  will  only  be  performed  once 
as  no  cycle  is  indicated. 
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Order      Input  to  MRA     In  MRA     Transfer 

from  MRA 

ica,  el  tr  ir 

idv,  e2  4-  180  tt/i80 

idv,  e2  -T-  180  tt/(i80)2 

its,  e°  tt/(180)2  tt/(180)2 

This  group  of  orders  will  calculate  the  constant  (tt/180)  and 

store  the  result  in  e8.  The  device  used  here  is  that  when  a  number  is 

transferred  from  the  MRA  to  a  register  the  MRA.  retains  the  number.  The 

fourth  power  is  formed  by  multiplying  the  number  in  the  MRA  by  the  number 

that  was  transferred  which  is  the  same  as  squaring  the  number  in  the 

>2 


MRA,  in  this  case  (tt/180)' 


Order 

Input  to  MRA 

In  MRA 

Transfer 
from  MRA 

imr,  el 

X    TT 

(tt/180)2 

its,  e8 

(tt/180)2 

(tt/180)2 

imr,  e8 

X  (tt/180)2 

(rt/lQ0)h 

its,  e8 

(rr/180)1* 

(it/ISO)11 

This  group  of  orders  will  calculate  G  and  store  the  result 

w 


in  flu 
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Order      Input  to  MRA      In  MRA     Transfer 

from  MBA 

ica,  f  3  r  -> 


W  E. 


.^ 


imr,   cl  *   W  -E — =- 

1  1 

njnr,   al  *   L  _£ — 

W  L2E  I 

imr,   al  XL  -^ ^ — 

•     ,       * l        vM 

imr,   al  X    L  -*- =j — 

idv,  bl  -  EI  -B°AA 

O  ^       f-JL_^  P    Oil     /    TT     >ll 

*» e8       *  w    rrvr  w 

wiWo 

..       ri  Vo^l  /  tr  ^  C 

wiWo 


This  group  of  orders  will  calculate  C  and  store  the  result  in 
f6.  In  this  operation  the  square  root  sub-routine  is  used.  The  third 
order  is  used  to  enter  this  sub-routine,  which  is  in  the  Whirlwind  mode, 
requiring  the  IN  order  to  return  to  the  GS  mode  once  the  square  root  has 
been  calculated* 
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Order 

Input  to  MRA 

In  MRA 

ica,  e$ 

g 

Z 

imr,  f3 

W  I? 
11 

"A 

sp,  si 

wiLi 

Transfer 
from  MRA. 


IN 


wA 


imr,   e9 


X     (180)2 


idv,   elO 


-r  2 


2(180)' 


N 


\\g 


11 


its,   f6 


This  group  of  orders  will  calculate  C  and  store  the  result 
in  f?.  The  operation  is  performed  twice,  since  m2  is  2,  because  this 
constant  is  used  to  determine  j3  and  fL0  If  more  spans  were  used  this 
operation  would  be  performed  a  number  of  times  equal  to  two  less  than 
the  number  of  spans.  The  square  root  sub-routine  is  used  twice  in 
succession  here  to  form  the  fourth  root  of  a  number. 
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Order       Input  to  MRA.       In  MRA       Transfer 

from  MRA 


m2,   icr 

hS,   ica,   f3  k  ^ 


IN 


¥l 

Wxl3 

E  I 

W  L^ 

idv,   f>2+c  -^   -^-  \ 

2  2  112  2 


sp,   si 
sp,   si 


kklh^ 


2^2 
its,   f7+c  \| o 

WlW2 


ict,   h5> 


This  group  of  orders  will  calculate  C  ,  and  store  the  result 
in  f5.  This  operation  is  to  be  performed  three  times  as  indicated  by 
the  first  ordere  The  flad  ml  was  not  used  here  because  the  number  of 
times  this  operation  is  performed  depends  on  the  number  of  vibrational 
orders  being  investigated  and  is  independent  of  the  number  of  spans 0 
Usually  three  orders  of  vibration  will  be  investigated, (synchronous, 
forward  and  backward  precession).  The  flad  h6  is  the  beginning  of 
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the  first  loop  of  the  program.  This  means  that  when  the  program  has 
been  run  through  once  it  will  automatically  come  back  to  h6  and  do  the 
entire  program  again  using  the  second  value  of  e3.  This  loop  will  be 
performed  three  times  as  indicated  by  the  first  instruction. 


Order       Input  to  MRA         In  MRA         Transfer 

from  MRA 


icr,    3 


h6,   ica,   f3  , 

T.7    lO 


hh  ^h 


imr,   eh  *  J 


W  L3 
1  1 


E.LJL 


imr,   al  X  L 


W  L3 
W1L1 


idv,  bl  tEJ 


¥lJL0 

ETIL3 
0  0  11 


EI  JL       , 

^                          v    f1         1\               110    /l         In 
imr,    e3+c  X  («■♦—;        ~ ?  ("o  +  TV 

2        ^  EI  WL3    2        q 


0  0  11 


_     i         E_  I_  JL_       _        ..  . 

Q  V      f     n      ^  110        /I  lw     IT      sll 

^  e8  W     77773  (7 +  q)  W 

0  0  1  1 

E.  I  JL         _  . 

its,  fs  -^-t'I*^^       °- 
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This  group  of  orders  assigns  a  value  to  the  flad  ml;.  The  value 
used  is  that  contained  in  d3  which  for  this  problem  was  135>«  This  number 
was  selected  because  we  wanted  j3  to  vary  from  90  to  36O  in  as  small  an 
increment  as  possible  to  increase  the  accuracy.  Due  to  limitations  on 
the  magnetic  tape  output  the  smallest  increment  feasible  was  2  .  The 
first  instruction  is  used  to  select  the  Whirlwind  mode  since  the  trans- 
fer of  digits  as  indicated  by  the  third  instruction  cannot  be  performed 
in  the  cs  mode.  The  fourth  order  is  used  to  return  to  the  cs  mode 
before  proceeding  with  the  remainder  of  the  program. 

The  last  instruction  is  used  to  select  a  new  cycle  counter. 
This  instruction  permits  a  loop  within  a  loop  to  be  performed.  In  the 
preceeding  group  of  instructions  counter  zero  was  used  and  is  stepped 
once  each  time  the  sequence  is  performed,  but  a  different  counter  is 
used  to  allow  the  following  sequence  to  be  performed,  and  counted,  a 
specified  number  of  tirnes  without  interfering  with  the  previous  counter. 

Order      Input  to  MM      In  HRA      Transfer 

from  MRA 

OUT 
ca,  d3        135  135 

td,  mU  135         13< 

IN 
is,  cl 

This  group  of  orders  will  calculate  the  frequency  and  print 
the  result  on  the  magnetic  tape  output  device.  This  operation  will  be 
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performed  135  times,  indicated  by  the  first  order,  for  each  cycle  of 

2 
the  main  program.  The  second  through  the  fifth  orders  calculate  P 

and  p  and  store  them  in  gl  and  g2  respectively*  The  frequency  is  then 

2 
found  by  multiplying  C  by  p  .  The  last  instruction  is  an  order  to  the 

computer  selecting  the  output  device  and  giving  the  information  necessary 

to  print  the  frequencies  in  a  column,  the  number  of  digits  desired,  etc. 


Order 

Input  to  MRA 

In  MRA 

Transfer 
from  MRA 

mil,  icr 

h7,  ica,  e6 

90° 

90° 

imr,  e6 

X  90° 

(90°)2 

its,  gl 

(900)2 

h 

imr,  gl 

*  (90°)2 

(90°)^ 

its,  g2 

(90°)^ 

\ 

-ina   -P6 

tt    IW 

tt    |^8 

2(180)2|  W^  2(180)2/  W1I^ 


to,  gl  x  (90O)2         JH2221 W 

2(i8or/w1i^ 

imoap  123l|..56ssss 

This  group  of  orders  will  calculate  the  moment  at  the  inboard 
end  of  the  tail  shaft,  mn-. ,  and  store  the  results  in  g5«  This  operation 
will  also  be  performed  135  times  for  each  cycle  of  the  main  program. 
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&     and  <h T,  both  functions  of  p  are  comouted  first  and  stored  in  g3  and  gl). 
respectively*  Using  the  formula  in  the  flow-chart  they  are  combined  to 
calculate  m m •  It  should  be  noted  that  this  and  the  preceeding  calcula- 
tions do  not  require  more  than  two  double  length  registers  to  store  the 
results.   In  the  frequency  calculations  the  results  are  printed  out  as 
soon  as  they  are  obtained.  Here  the  results  stored  in  g£  are  used  immedi- 
ately in  the  subsequent  calculations  after  which  they  are  no  longer  desired. 
Therefore,  the  m„  calculated  is  lost  as  each  succeeding  moment  is  determined. 


Order       Input  to  HRA        In  MRA  Transfer 

from  MRA 

W  L?K  L      ,   W  L?E,  I. 


ioa  ih       iOiif'f    ^aii(«)U 
wA3Vo  18°      ViVo  18° 


h  \Ll\h  ,jta_>U 


w  ^    180 


L?EL 


•+  n  pill    /JLPxU  / 


1  10  0 


ETll3  ^        q     180  EIWI3  2        q     180 


0  0  11  0  0  1 


*i 


unr     ec 


p  T  y  T3  2        q     160 


EIJL       ,  ,  . 

..  ,  11       0     /l  lw     TT     vh  . 

E  I  W  L?  2  q    18°  J 

0  0  1"! 


Order 


Input  to  MRA 


In  MRA 


(continued) 


>o. 


Transfer 
iron*  MRA 


ica     g3 


w 


w 


idv     el2 


-5-3 


w 


its     tl 


A 


w 


3 


idv     el3 


±k 


_J1 
12 


imr     git 


*  1>. 


iad 

Lh9 

+  1 

isu 

tl 

3 

isu 

gh 

-i 

its 

t2 

ica 

tl 

K 

w 


Vj 

12 

1 

4A 

12 

1 

12 

1 

12 

1 

Vj 

w 


rf 


--^--f 


4 


12 


-  —  -  p 


i  +  -j^- 

J  12 


A 


H 


unr 


sU 


*  /■ 


isu     gU 


-/■ 


J  -  jzT 
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Order 


isu  g3 


Input  to  MRA         In  MRA 
(continued) 


Transfer 
from  MRA 


-A 


j  - 


w 


J         w 


idv     t2 


imr     bl 


idv     al 


its     g£ 


■M  + 


A 


V 


12 


h  - 


t,  -  A 


"Vo 


0 


1  + 


J7j   K 


w 


12 


0  0 


3      "'j 

w 

L+     12     " 

w 
3 

* 

Vo 

3        "   'j   ~  K 

Lo 

,                VJ                  W                 / 

li  +     12     "    3    "  'd 

Vo 

4A 

V-4,-4, 

Lo 

t        12     "   3    "  'jj 

m, 


01 


This  group  of  orders  will  first  convert  the  value  of  p  currently 
stored  in  e6  from  degrees  into  radians  and  store  the  result  in  e7«  The 
last  two  orders  will  place  the  initial  value  of  9  _ ,  which  is  assumed  as 


unity,  in  g9. 
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Order 

ica  e6 

imr  el 

idv  e2 

its  e? 

ica  lh9 

its  g? 


Input  to  MRA 


In  MRA. 


180 


180 

180 


Transfer 
from  MRA 


Jt£_ 
180 


9 


01 


This  group  of  orders  will  calculate  various  parameters,  such 
as  cos  (3,  sin  {3  and  store  them  temporarily  for  future  use  in  the  calcu- 
lation of  end  moments.  The  first  instruction  selects  a  third  counter 
to  be  used  for  this  and  subsequent  calculations 0  This  is  done  for  the 
same  reason  as  before,  to  permit  this  sequence  to  cycle  three  times, 
once  for  each  inboard  span,  without  affecting  the  cycle  count  of  the 
previous  loop.  The  second  instruction  sets  the  number  of  times  the  loop 
is  to  cycle.  This  number  equals  the  number  of  spans  since  m3  takes  its 
valud  from  d3  which  is  the  flad  assigned  for  the  number  of  spans e  The 
manner  in  which  the  subroutines  are  utilized  is  shown  here  and  their 
advantage  becomes  evident. 
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Order       Input  to  MRA       In  MRA       Transfer 

from  MRA 


isc  2 

m3,  icr 


h8,  ica 

el 

Pi 

h 

isp 

s2 

cos  P 

its 

ul 

cos  P- 

cos  p.. 

ica 

el 

h 

h 

isp 

ls2 

sin  p. 

its 

u2 

sin  px 

sinP1 

ica 

el 

h 

h 

isp 

s3 

sinh  p_ 

its 

u3 

sinh  p_ 

sinh  p 

ica 

e7 

h 

Pi 

isp 

£s3 

cosh  p_ 

its 

uU 

cosh  p 

cosh  p 

ica 

u3 

sinh  p 

sinh  p 

isu 

u2 

-sin  p.. 

sinh  p  -sin 

h 

its 

x6 

sinh  p  -sin 

h 

sinh  p  -i 

ica 

ul 

cos  p 

cos  p 

imr 

u3 

X  sinh  pn 

cos  p_  sinh  | 

\ 

its 

u6 

cos  8  sinh  1 

\ 

cos  P„  s: 

Input  to  MRA 

In  MM 

Transfer 
from  MRA 

sin  pi 

sin  p1 
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This  group  of  orders  will  calculate  the  parameters  si,  s2,  s3, 
and  store  the  results  in  g6,  g7,  and  g8  respectively.  This  operation  will 
be  repeated  three  times,  once  for  each  span,  since  the  first  two  orders 
of  the  previous  section  are  still  in  effect. 

Order 

ica  u2 

imr  uU        *  cosh  p,        sin  p  cosh  ,6_ 

isu  u6     -cos  Bn  sinh  pn  sin  p.,  cosh  B  -cos  B,  sinh  3, 
rl      rl     '  1      '1     '1      '1 

sin  p  cosh  B  -cos  6  sinh  p_ 

idv  u£    -rsinh  p  -sin  pn  •   u  q        • — <T^ ~ 

rl     '  1       sinh  8  -sin  p_ 

sin  p  cosh  p  -cos  p  sinh  p 

its  S6  •   u  a1    ■      o  ] ^   S_ 

&  sinh  p.. -sin  p..  1 

sinh  p1 

sinh  p  -sin  B 

p  (sinh  p  -sin  B  ) 
'1      'l     '1 

p1(sinh  B  -sin  B  )   B  (sinh  B  -sin  8  ) 

-cos  p 
-cos  p  cosh  p 
1-cos  p_  cosh  p 
1-cos  p  cosh  p. 


ica 

u3 

sinh  p. 

isu 

u2 

-sin  p1 

imr 

e7 

X?! 

its 

u5 

ics 

ul 

-cos  p_ 

imr 

vh 

*cosh  P1 

iad 

lh9 

+1 

idv  u5  -r  p..  (sinh  B  -sin  p_  )    /  .   Q : — r— r 

rl      '1     rl   p  (sinh  B  -sin  B  ) 
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Order         Input  to  MRA         In  MRA         Transfer 

from  MRA 

(continued) 

1-cos  p  cosh  p 
its  g7  P^sinh  p^-sin  pi)     S2 

ica  u3  sinh  p  sinh  p_ 

isu  u2  -sin  p..        sinh  p  -sin  p_ 

its  v$  sinh  p  -sin  p    sinh  p  -sin  p. 

1  J.        J. 

ica  elO  2                2 

imr  el  *  P1               2px 

imr  u2  *  sin  p           2p  sin  p 

imr  u3  *  sinh  p.,  2p..  sin  p  sinh  p 

2p_  sin  p  sinh  p 

idv  u5>       -j-sinh  p  -sin  p_  .  ,  ,  — : — ■=. 

rl     rl    sinh  {L-sin  p 

2p  sin  p  sinh  p 

its  g8  .  ,  0 : o S 

&  sinh  p  -sin  p  3 


This  group  of  orders  will  calculate  the  end  moments  and  the 
angle  of  rotation  at  the  inboard  end  of  the  shaft  selected.  The  previous 
results  stored  in  g5,  namely  m _ ,  and  previously  calculated  values  of 
©n„,  si,  s2,  s3  are  combined  into  the  formulas  shown  in  the  flow  chart 
to  determine  the  values  at  the  end  of  the  second  section,  m-,9.   This 
process  is  repeated  three  times  to  produce  m,  ,  the  moment  we  are 
looking  for.  9_,  is  also  determined  although  it  is  not  utilized. 
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The  method  of  saving  registers  that  is  used  in  this  section  should  be 
noted.  mn1  is  in  g5.  As  each  successive  moment  is  calculated  it  too 
is  put  in  g5.  The  preceeding  moment  is  lost  by  this  process,  but  it 
is  of  no  importance  since  only  the  final  end  moment  is  required. 

The  last  six  orders  have  a  different  function.  The  first  five 
of  these  are  utilized  to  calculate  the  value  of  p  for  the  next  span 
which  is  used  in  the  calculation  of  the  moment  and  si,  s2,  s3.  The 
relation  of  0  and  £L  is  shown  in  the  flow  chart. 

The  last  order  is  used  to  cycle  the  program  back  to  h8.  ict 
causes  the  computer  to  return  to  h8  and  at  the  same  time  increases  the 
index  of  the  cycle  counter  by  one.  When  the  index  equals  the  number 
previously  set  in  the  criterion  by  an  icr  instruction  ict  will  reset 
the  index  to  +0  and  pass  control  to  the  instruction  immediately  following. 


Order 


Input  to  MRA 


In  MRA 


Transfer 
from  MRA 


ica  g5 

its  glO 

ics  g5 

imr  g6 

its  ul 

ics  f2+c 

imr  g8 

imr  g° 


m 


01 


-m 


01 


*S_ 


-C 


XS. 


*e 


01 


m 


01 


m, 


01 


-m 


01 


*moisi 


'moisi 


-c 


•c  s_ 

a  3 


•°aS3e01 


m 


01 


■moisi 
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Order 

Input  to  MRA 

In  MRA 

Transfer 
from  MRA 

(continued) 

iad  ul 

-%A 

-°aS3e01-m01Sl 

its  g5 

-CaS3e01-m01Sl 

^2 

ics  glO 

"%! 

~™01 

imr  g7 

as2 

-m01S2 

idv  f2+c 

^C 

-m01S2 

isu  g6 


-s„ 


"m01S2 


-  S. 


its  g9 


-m01S2 


-  S. 


9 


12 


ica 

f7+c 

C 

imr 

e6 

Kp 

imr 

el 

X.TT 

idv 

e2 

tL 

its 

e7 

ict 

h8 

TT 

i8o  y 

180  V 


This  group  of  orders  vrill  print  out  the  contents  of  g£,  which 
are  the  moments  at  the  inboard  end  of  the  shafto  The  first  instruction 


Ik. 


is  used  to  return  the  counting  of  cycles  back  to  the  135>»  This  is 
done  since  we  want  135  end  moments  printed,  or  one  for  each  (3  value 
used. 


Order 


isc  1 


ica  g5 


Input  to  MRA 


in 


3U 


imoa+p  123U.56  *  10~  c 


In  MRA 


m. 


% 


Transfer 

from  MRA 


This  group  of  orders  does  several  things  which  will  be  dis- 
cussed in  the  order  of  their  occurance0 

The  first  three  orders  are  used  to  increase  the  value  of  {3 
in  e6  by  2  before  the  next  cycle  starts.  This  operation  is  performed 
135  times  also,  the  [3  starts  at  90  and  goes  to  360  in  steps  of  2  , 

The  fourth  order  is  used  to  break  the  normal  sequence  and 
add  a  correction  before  returning  to  the  main  program,  a  process  known 
as  patchinge  The  correction  added  is  given  below.  The  first  instruc- 
tion is  merely  a  new  flad  used  to  get  the  number  +1.0  in  the  computer. 
The  second  instruction  increases  the  index  of  counter  1  by  one  and  per- 
forms the  same  function  as  explained  above,  namely,  permit  this  part  of 
the  program  to  cycle  135  times.  The  third  and  fourth  instructions  are 
merely  to  print  a  series  of  zeros  in  the  print  out  to  separate  the 
results  into  blocks  of  135  lines  of  frequency  vs.  moment  for  each  order 
being  investigated.  The  last  instruction  is  to  pass  control  to  the 
next  patch. 


7*. 


Order 

Input  to  MRA 

In  MRA 

Transfer 

from  MRA 

ica 

e6 

P 

P 

iad 

elO 

+2 

p+2 

its 

e6 

p+2 

P+2 

isp 

3h9 

isp 

hl2 

ict 

h6 

h9,  i  STOP 

START  AT 

hi 

lh91+1.0 

ict 

h7 

ica 

lh9 

+lo0 

+1.0 

imoa+123i|CC 

isp- 

2h9 

This  group  of  orders  will  do  several  things.  First  it  passes 
control  to  the  flad  hl2  which  introduces  the  subroutine  of  resetting 
the  value  of  $   to  90  at  the  completion  of  the  program  for  one  order. 
Then  counter  zero  is  selected  to  begin  calculations  for  the  next  order. 
The  instruction  OUT  is  used  to  enter  the  Whirlwind  mode  and  the  next 
three  instructions  are  used  to  increase  the  address  of  I4I16  by  two.  This 
is  normally  done  automatically  by  the  use  of  counters  and  +c  at  the  end 
of  each  instruction.  However,  difficulty  was  encountered  and  this 
"manual"  stepping  was  used  to  save  time.  The  last  two  instructions  are 
used  to  return  to  the  CS  mode  and  the  main  program  at  one  instruction 
before  h9. 
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Order 

hll,  +90,0 
-2h9/isp  hl2 
hl2,  ica  hll 
its  e6 
isc  0 
OUT 
calih6 
ade  11 
ade  11 
tsbh6 

IN 
isp  h?-l 


Input  to  MRA 


In  MRA 


+90.0 


285 
+1 
+1 


+90,0 
+90,0 


285 
286 
287 
287 


Transfer 
from  MRA 


+90,0 


I4I16 


The  main  program  is  reentered  as  explained  above  and  the  last 
three  instructions  listed  below  are  performed.  The  first  Increases  the 
index  of  counter  zero  by  one  and  when  the  entire  program  has  been  repeated 
three  times  passes  control  to  the  second  instruction  which  stops  the 
computer. 

The  last  instruction  is  used  to  start  the  calculations  at  the 
proper  place.  It  is  always  the  last  instruction  fed  into  the  computer 
since  the  computer  scans  the  entire  program,  locates  the  flads  and  stores 
the  information  before  it  s tarts e   If  this  instruction  were  anywhere 
else  the  computer  would  start  computing  before  all  the  information  was 
fed  into  it  and  properly  stored* 
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Order 

ict  h6 

h9,  i  STOP 
START  AT  hi 


Input  to  MRA 


In  MRA 


Transfer 
from  MRA 
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PROGRAM  OF  STRESS  CALCULATIONS  OF  NON- SYMMETRICALLY 
LOADED  TRANSVERSE  WEB  FRAME 

The  working  program  of  the  web  frame  stress  calculations  will 
consist  of  a  series  of  sub-routines  joined  together  by  a  sequence  control 
program  (see  Figure  V).  To  facilitate  programming  and  program  alteration, 
the  instructions  and  input  data  were  placed  in  four  separate  sections 
and  these  sections  were  represented  by  the  flexowriter  tape  number  that 
contained  the  instructions  on  data.   Scantling  data  and  other  information 
relating  to  a  particular  problem  were  placed  in  tape  359-217-100 «,  The 
sub-routines  were  placed  in  359-217-101,  and,  the  sequence  control  pro- 
grams were  placed  in  359-217-221  and  359-217-222.  Calculations  for  a 
different  though  similar  ship  will  require  the  complete  changing  of 
tape  359-217-100.  Tape  359-217-101  is  general  for  similar  problems.  Tapes 
359-217-221  and  222  determine  the  nature  of  the  calculations,  that  Is 
itfhether  a  check  calculation  or  a  design  calculation  will  be  calculated,, 

In  the  instructions  for  the  input  data  in  tape  359-217-100  abso- 
lute addresses  were  assigned  to  flads .   The  absolute  addresses  were  necessary 
for  block  in,  bl    out  procedures  used  in.  the  auxiliary  storage  drum.  The 
use  of  flads  is  still  retained  due  to  the  ease  of  writing  and  categorizing 
letters o  The  flad  and  absolute  address  assignments  for  storage  of 
generated  data  were  also  placed  In  tape  359-217-101. 

The  following  list  of  data  are  the  specific  input  data  for  a 
DP  type  vessel.  Due  to  the  quantity  of  numbers  in  each  ategory,  only  a 
sample  of  the  total  group  is  presented  belowc 
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32  I  xL,  +00.00' 
+3c266' 
+6.508' 
+9o917 

Total  of  25  double 

+0.000 
82  |  yl,   +00.00' 

25  numbers 

+26.33 
332|bl,   +2LI..00 

25  numbers 

+12.00 
282|b2,   +0.3125 

25  numbers 


fc  Tape  359-217-100 
(2i|,6) 

Assigns  the  first  value  of  the  x-coordinate  of 
the  shell  to  register  32  which  is  the  first 
register  in  the  case  memory.  Note  that  the 
coordinates  are  in  feet  measured  from  the  origin 
at  the  keel  centerline,, 

length  numbers  which  occupy  50  registers,-, 

Value  at  the  main  deck  centerline. 

First  value  of  y-coordinates  of  the  shell. 


+0.2,400 


1|82 


b3,   +12.0 

25  numbers 

luO 

te|blis   +0.500 

25  numbers 

+0,250 
382 | b5,  +0.7500 

25  numbers 

+0,375 


y-coordinate  of  main  deck  centerline. 

Value  of  web  depth  (d  )  at  the  origin  measured 

w 

in  inches.  First  value  stored  in  register  332, 


Value  at  main-deck  centerline. 

Thickness  of  web  (t  )  measured  in  inches. 

w 


Width  of  flange  (w^,)  measured  in  inche; 


Flange  thickness  (t„)  measured  in  inches, 


Shell  thickness  (t  )  in  inches. 
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632  e5>j  ~97h.O       Weight  and  shear  loads  at  each  station.  The 

weights  are  determined  in  averaging  the  weight  of 
50  numbers       vessel  and  shear  for  one  foot  of  length,  ^he   length 

per  foot  times  the  number  of  feet  of  girth  between 
-97luO       stations  gives  the  weight  load  per  station. 

932  e25  +25o5l7       Height  of  water  pressure  head.   In  this  problem 

the  design  water  head  was  the  height  of  the  deck 
edge. 

978  el,  +8.0         Longitudinal  distance  between  web  frames.  The 

stresses  c  omputed  were  for  the  loads  on  one  foot 
of  logitudinal  distance.  Each  web  frane  is  assumed 
to  support  the  loads  imposed  in  the  longitudinal 
distance  between  frames. 

pal  =  25  Present  parameter  which  represents  the  number  of 

stations  from  the  origin  deck  centerline. 

pa2  =  18  Present  parameter  which  represents  the  number  of 

stations  from  the  origin  to  the  deck  edge, 

pa3  =  h  Number  of  stations  from  the  origin  to  the  bottom 

of  the  starboard  stanchion. 

pal;  =  k  Number  of  stations  from  the  origin  to  the  bottom 

of  pori.  stanchion  when  measuring  counter  clockwise. 

pa£  =  18  Number  of  stations  between  bottom  and  top  of  the 

starboard  stanchion. 

pa6  =  18  Number  of  stations  between  bottom  and  top  of  the 

port  stanchion. 


The  following  instructions  are  the  sequence  control  instructions 
and  serve  to  utilize  various  sub-routines  to  perform  the  calculations  of 
stresses  for  a  web  frame  of  one  particular  set  of  scantlings  and  load 
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conditions.  To  alter  the  web  frame  and  to  produce  a  rational  design, 
another  set  of  instructions  designated  tape  number  359-217-222  will  be 
required. 

fc   359-217-221  (Name  of  tape) 
(2U,6) 

(Selects  number  arrangement  of  double  length  registers. 
This  arrangement  utilizes  2k   binary  numbers  for  the 
mantissa  and  6  binary  digits  for  the  exponent  of  the 
number). 

pi,  isp  gl  gl  forms  sin  9,  cos  9  and  As 

isp  g2-2  g2-2  forms  horizontal,  and  vertical  water  loads 

H  ,  V  . 
w'  w 

isp  glO-2  glO-2  forms  the  total  vertical  loads  by  summing 

the  values  of  the  weights  and  vertical  water 
loads. 

isp  g3  g3  forms  sin  ^  and  cos  <f>, 

isp  glj.  Forms  the  reciprocal  of  the  moment  of  inertia 

(1/Im  a  )f   of  the  section  (A  )  and  distance 
from  the  shell  to  the  neutral  axis  (d,T  .  ). 

OUT  The  following  orders  are  block  out  instructions 

for  transferring  information  from  the  case 
memory  to  the  auxiliary  drum. 

20U8  (Address  of  first  register  of  the  drum 
storage  block). 

Selects  case  output  routine  to  the  auxiliary 
drum. 

800  (Number  of  registers  necessary  for  storage). 

(Case  memory  address  of  the  first  register  being 
blocked  out). 


Offsets  the  x  and  y  coordinates  of  the  neutral 
axis  from  the  shell  coordinates. 


ca 

c8 

si 

967 

ca 

cl2 

bO 

32 

IN 

isp 

g6-2 
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isp  gl  Forms  sin  9  and.  cos  9,  and  As  of  neutral 

axis. 

Forms   sin  /,    and  cos  <f>  of  neutral  axis. 

+0.0  +0.0 

+0,0 
+0.0 

Forms  m/I  at  each  station  due  to  applied  loads. 

Forms  integrals  of  m  ,  p  ,  q  ,  R  and  S,  into 
a  matrix  (See  Figure  VI). 

Uses  0 rout's  method  to  solve  the  matrix. 

Alters  V  and  H  to  include  redundants  V-q  , 
H+p  ,  R+V  and  V+S. 


isp 

g3 

ica 

c6 

its 

til 

its 

b3 

isp 

g7 

isp 

S9 

isp 

g23-2 

isp 

g2lrl 

ica 

dl+£0 

its 

tU 

imr 

d3 

its 

b3 

isp 

g7 

OUT 

ca 

clO 

si 

963 

ca 

cl7 

bi 

332 

+mQ  +mQ 


'    0  0  0 


+mo 


mo/3:o 


Forms  moments  at  each  station.  This  moment 
includes  the  effect  of  the  redundants 


2 3ii8  ~  Drum  address  of  the  first  register  of  the 
drum  storage. 

Selects  "block-in"  procedure  for  transferring  data 
from  the  drum  to  the  case  memory. 

+1^0  -  Number  of  registers  being  "blocked  in." 

Blocks  in  d  ,  t  and  t„. 
w      s     f 

Address  of  the  first  register  in  the  core  memory 
that  will  receive  the  first  number  moved  in  from 
the  drum. 

ca  ell  +27l|8  -  Address  of  the  first  register  of  the  drum 

storage. 
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si  963  Selects  "block  in"  routine  from  the  drum 

storage. 

ca  cl8  *$0   -  Blocks  in  50  registers  of  information 

0f  dN,A. 
bi  732  Address  of  the  first  register  in  the  core  memory 


that  will  contain  d^  . 


IN 


p2,   isp  g23>-l        Forms  (T ,  (K,   (T   and  (T   for  all  stations  around 

the  girth., 

The  following  instructions,  called  tape  35>9-217-2224  perform  the 
stress  calculations  for  a  web  frame  for  a  given  set  of  scantlings  and 
compare  the  computed  stresses  with  allowable  stresses.  If  the  computed 
stresses  do  not  fall  within  the  range  of  allowable  stresses,  then  the  instruc- 
tions below  will  alter  the  scantlings  of  the  web  frame  to  produce  the 
desired  stresses.  For  this  program  the  scantlings  are  altered  by  increasing 
or  decreasing  the  width  of  the  flange.  Notice  that  tape  359-217-222  repeats 
all  the  instructions  contained  in  359-217-221  and  in  addition  contains  a 
method  of  altering  the  flange.  While  this  may  seem  unnecessarily  repetitious, 
two  complete  sequence  programs  were  deemed  more  convenient  to  handle  and  less 
susceptible  to  error  than  one  partial  tape  with  a  connective  addition. 

fc  359-217-222 
(2li,6) 

pi,  isp  gl  gl  forms  sin  ©,  cos  9  and  As. 

isp  g2-2  g2-2  forms  V  and  H  . 

w     w 

isp  glO-2  glO-2  forms  V 

isp  g3  g3  forms  sin  ^  and  cos  /. 

isp  gk  gk   forms  lAN#A#,  At,  and  d^A 
OUT 


ca 

c8 

si 

967 

ca 

cl2 

bO 

32 

IN 

isp 

g6-2 

isp 

gl 

isp 

g3 

ica 

c6 

its 

tli 

its 

b3 

isp 

g7 

isp 

S9 

isp 

g23-2 

isp 

g2li-l 

ica 

dl+50 

its 

tb 

imr 

d3 

its 

b3 

isp 

g7 

OUT 

ca 

clO 

si 

963 

ca  c!7 

8U. 


20liB  -  Drum  Address* 

Selects  drum  output  routine. 

OOO  -  Nrmber  of  registers  being  transferred. 

Core  memory,  address  of  first  register  being 
blocked  out. 

Forms  (x^,  jr^). 

sin  Q9   cos  0,  and  As  of  neutral  axis. 

sin  $  and  cos  fi   or  neutral  axis, 

+0,0  +0.0 

+0.0 

+0.0 

Forms  m  /I  at  each  station 
n' 

Forms  integrals  of  m^9   p  ,  a,  R  and  S.   (See 
Figure  VI).        v       °'  ^ 

Matrix  solution. 

Alters  loads  to  include  redundants.  q«,  pn,  R  and  Sf 

+m0  +m0 


^0  "(A) 


+mQ 


mo/]:o 


Forms  moments  at  e  ach  station* 

23I18  ~  Drum  address. 

Selects  "block  in"  from  drum  internal  routine. 

+15>0  -  Number  of  registers  being  blocked  in. 


8$. 


bi  332  Address  of  first  register  in  core  memory  to 

receive  block  in  from  drum, 

ca  ell  27U8-Drum  address  of  first  register  in 

,rblock  in"  group, 

si  963  Selects  "block  in"  internal  routine. 

ca  cl8  +5>0-Number  of  registers  being  blocked  in. 

bi  732  Address  of  first  register  in  core  memory  to 

receive  blocked  in  data. 

IN 

isp  g25>-l  Forms  er  .  <3",  ,    <^  and   <S   for  all   station. 

a       b       f  s 

OUT 

ca  cl°  +2li°8-Drum  address  of  first  register  to  be 

blocked  in. 

si  963  Selects  "block  in"  internal  routine, 

ca  cl8  +50-Number  of  registers  being  blocked  in0 

bi  I482  Core  address  of  first  register  receiving  blocked 

in  data. 

aO  976  Forms  cornter  '00   limit  number  of  alteration 

cp  p2-l±  cycles  to  3(l  +  number  in  976) . 

spO  Stops  computer  after  completion  of  3  cycles 

IN 

ica  c3  -1.0  -lo.O 

its  t3  -1.0 

Note  that  the  above  counter  was  used  for  the  test  runs  of  this 
program.  In  any  case  the  alteration  procedure  should  be  limited  to  pre- 
vent the  possibility  of  the  computer  obtaining  an  unreasonable  answer.  As 
an  example,  should  the  web  frame  be  to  large  the  web  flange  can  be  reduced 
by  predetermined  amounts.  However,  large  reductions  could  give  a  negative 
length  of  flange.  A  future  programmer  may  want  to  change  the  value  used 
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in  the  counter  and  therefore  the  number  of  cycles e  However  a  counter  must 
be  used  in  all  cases. 

Now  it  becomes  necessary  to  compare  the  values  of  the  computed 
stresses  with  allowable  values.   To  do  this  the  computed  values  must  be 
plus  and  a  small  program  must  be  written  to  insure  that  all  values  are 
plus.  Also,  the  maximum  stress  value  at  each  point  of  equal  scantlings 
must  be  determined.  Remember  that  the  web  frame  is  symmetrical  port  and 
starboard,  but  the  values  of  stress  are  not  symmetrical 6  After  determining 
the  maximum  value  of  stress  at  each  station,  the  values  are  compared  with 
allowable  values  and  the  scantlings  altered. 

ict  pal 

p2,    ica  e5+c  +  CO  ^(starboard  side) 

icp  ml  If  CT  is  minus  perform  sub-routine  ml.      If 

plus,    do  next  instruction, 

its   tl  +  <7~ 

f 

ica  e£+pb8+pb8  +(T  +CT (port   side). 

icp  ml  Checks  sign  of  CT* 

its  t2  +^ 

isu  tl  -Cm  starboard)   (/.(port)  -  (^(starboard) 

icp  m2  Checks  sign  of  (/"(port)  -  (T  (starboard)  and  if 

minus,  goes  to  m2, 

ica  t2  +  (7"  +(T 

f  f 

p3,  isu  c20  -(+20,000)       CT  -(+20,000) (Lower  value  of  stress 

limit). 

icp  m3  If  (0y-20,000)  is  minus,  go  to  m3  to  reduce 

scantling  values. 
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isii  c21  -8000  ((T-20o 000) -8000  (8000  /in2 

f 

is  the  range  of  stresses). 

icp  pU  If  value  is  minus,  perform  sub-routine  pij.  which 

moves  forward  to  the  next  station,, 

ica  b3+c  +w^  +w„ 

f  f 

iad  cl  +  .5  w  +.5   (increases  flange  width). 

its  b3+c  w  +0£ 

ics  c3  -(-1.0)  +1.0 

its  t3  +1.0 

pU,    OUT 

ca  p2+3  ica  e£+pb8+pb8       ica  e5+pb8+pb8 

su  c7  -2  ica  e£+pb8+pb8-2 

ts  p2+3  ica  e5+pb8+pb8-2 

IN 
ict  p2 

ica  t3  -1.0  -1.0 

icp  p6 

OUT 

ca  cl9  +2J4.98*Drum  address  register. 

si  967  Selects  drum  block  out  procedure. 

ca  cl8  +50-Number  of  registers  being  blocked  out, 

bo  I482  Core  memory  address  of  first  register  of  group 

being  blocked  out0 

ca  c8  +20l+8-Drum  address  register  of  first  number  in 

a  group  being  blocked  in0 

si  963  Selects  drum  block  in  procedure. 

ca  cl2  800-Number  of  registers  being  blocked  in0 
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ca  cl2  +800  +800 

ad  cl7  +1^0  +950 

ad  c2  +30  +980 

ts  t8  +980 

ca  t8  +980 

si  963  Selects  drum  block  in  procedure. 

ca  cl2  +800 

bi  980 


The  two  block  in  procedures  put  numbers  and  instructions 
back  into  the  core  memory  in  preparation  for  another  cycle  of  stress 
calculations  of  the  altered  scantlings «   The  block  out  instruction 
prior  to  the  block  in  instructions  puts  new  values  of  web  flange 
width  into  the  columns  of  numbers  on  the  drum.   The  first  block  in 
instruction  moves  numbers  into  the  core  memory c   The  second  block 
in  instruction  came  as  a  result  of  "debugging"  operations  when  it 
was  determined  necessary  to  replace  all  the  addresses  in  the  in- 
structions back  to  their  starting  value  involving  the  port  side  it 
was  necessary  for  the  programmer  to  write  instructions  that  would 
decrease  the  addresses  of  the  registers  containing  the  starboard  side 
values.   In  other  words  it  was  necessary  to  move  up  a  column  of 
numbers  in  the  core  memory.  At  the  beginning  of  the  second  cycle 
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the  instruction  addresses  were  in  error  by  a  number  equal  to  the  number 
of  cycles  in  that  particular  sub-routine.  By  blocking;  in  the  original 
instructions  which  remain  on  the  drum  into  their  assigned  positions 
all  the  addresses  were  reset  to  their  original  value. 

For  cycles  within  a  basic  cycle  of  operations  it  is  necessary 
to  write  a  reset  routine  that  returns  the  addresses  back  to  their 
original  values  at  the  end  of  each  inner  cycle. 


IN 


isp  gk 
isp  pl+11 


Starts  cycle  at  glu 

These  two  instructions  skip  the  block  out 
procedure  in  the  main  cycle,  Since  the  block 
out  has  occured  during  the  first  cycle,  the 
values  remain  in  the  drum  and  need  not  be 
blocked  out  again. 


p6,  i  STOP 


The  following  sub-routines  are  used  in  conjunction  with 
the  comparison  and  alteration  sequence. 
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ml,    ita  ml+2 
imr  c3 

isp 

m2.    ica  tl 


x-loO 


♦<5 


(-icOx-cf )  =  +  <rf 

(Makes  CJT  plus). 


♦  (/„  ( s  tarboard ) 


isp  p3 

m3.    ica  b3+c 


+w, 


+w, 


isu  cl 


-.5 


W    -o5 


its  b3+c 


w  ~05 


isp  p3+7 

cl9, 
c20, 

c219 


+2li98-Drum  address, 


+20,000.0-Lower  limit  of  allowable  stress 
value e 

+8>000o0-Range  of  allowable  stress  values „ 

Note  that  for  this  problem  a  very  wide  range  of  values  was 
used  to  prevent  the  possibility  of  the  computer  seeking 
an  impossible  solution*     If  the  stress  change  caused  by 
the  flange  alteration  were  to  be  larger  than  the  stress 
range,    then  a  solution  would  be   impossible  with  this  program. 


976     -2 


i  START  AT  pi 

In  programming  it  is  necessary  to  use  various  constants  to  change 
units,  to  form  equations  and  to  change  generated  values  by  arithmetic 
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manipulations ,  These  constants  must  be  assigned  registers  or  flads  to 
be  of  use  in  the  program,,  The  constants  listed  below  Here  assigned  Had 
cl,  c2,  etc.,  the  c  being  useful  in  remembering  the  nature  of  the  input* 


cl,  +0„0 


Tape  fc   359-217-101 
(2l|,6) 


c2,  +30.0  A  constant  used  to  form  the  width  of  shell 

plating  which  is  considered  to  be  effective  in 
forming  the  web  beam, 

c3j   -loO  A  constant  useful  in  changing  the  sign  of  a 

value  in  the  MIA. 

clj.,  +1.5  Used  in  calculation  of  the  moment  of  inertia. 

c5,  +6iio0  Weight  of  a  cubic  foot  of  sea  water  in  lbs/ft  . 

c6,  +0.0  Used  to  clear  registers  where  the  value  in  the 

register  must  be  zero. 

c7,  +2  Whirlwind  number  used  in  changing  addresses  in  the 

interpreted  mode, 

c8,  +20I48 

c9    2850 

9  .  o  A  Whirlwind  number  used  to  assign  addresses  on 

C  9      97)  n  "bhe  auxiliary  drum. 

cl2,  +800  Number  of  registers  in  a  blockout 

cl3,  +216  procedure. 

cllu  +1  Whirlwind  number  used  in  changing  addresses. 

cl5,  +12.0  Converts  feet  or  inches  to  inches  or  feet. 

C16,  +0.707         1/  2  m   sin  h$°   or  cos  k$° « 

cl7,  +1^0  Number  of  registers  in  block  out 

cl8,  +50  procedures- 

tl,  through  t8,       +0,0  -  Temporary  registers  for  calculations. 
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The  following  preset  parameters  are  part  of  the  inputs 
do  not  vary  with  vessel  type.  Note  that  each  preset  parameter  listed  below 
is  generated  from  the  present  parameters  in  the  variable  input  tape. 


pbl  »  pal-1 

pb2  =  pal-2 

pb3  =  pa2-l 

pbli  =  pa3  +na3  -2 


Number  of  segments  in  the  half  •'irth  of  a 
web  frame. 

Two  less  than  the  number  of  stations  on  a  half  girth, 
Used  in  moment  formation. 

Number  of  segments  between  the  origin  and  deck  edge. 

Number  of  registers  from  the  origin  register  to 
the  bottom  of  the  starboard  stanchion. 


Also  contained  in  Tape  359-217-101  are  the  storage  assign- 
ments for  p-enerated  numbers. 


132  |  al,  +0.0 

DITTO 

I|8   +0.0 

182  |  a2,  +0.0 
DITTO 
U8   +0.0 

232  |  b$.   +0.0 
DITTO 

lj.8   +0o0 

532  |  eh,   +0*0 
DITTO 
U8   +0.0 

582  |  e3,  +0.0 
DITTO 
1|8   +0e0 

832  I  d3,  +0.0 
DITTO 
U8   +0c0 


Sine  of  the  angle  storage P  Note  that  this 
instruction  puts  zeros  in  25  double  length 
registers. 

Cosine  storage* 


Girth  segment  lengths  storage 


Storage  for  vertical  water  loads. 


Storage  for  horizontal  water  loads, 


Storage  for  the  inverse  of  the  moment  of  inertia 

1;I 

-1     N.A, 
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882|  d2,  +0.0 
DITTO 
U8   +0.0 

7321  dl,  +0.0 
DITTO 
U8   +0.0 

782|  si,  Library  sub-routine  for  finding  square  roots. 

The  remainder  of  tape  359-217-101  consists  of  the  various 
sub-routines  that  perform  the  functions  listed  in  the  flow  diagram, 
Figure  V0  A  discussion  of  each  sub-routine  is  placed  at  the  beginning  of 
the  list  of  instructions  that  comprise  the  sub-routine e 

The  form  of  presentation  of  the  following  sub-routines  was 
found  to  be  an  aid  in  the  organization  and  understanding  of  the  problem. 
Since  continuity  of  thought  is  so  important  to  the  writing  of  a  work- 
able program,  the  increased  effort  in  writing  the  flow  of  information 
into  the  KRA,  accumulation  of  data  in  the  MRA  and  the  values  trans- 
ferred from  the  KRA  was  considered  worthwhile.   Other  methods  of 
organizing  material  could  have  been  used;  however,  the  authors  found 
that  this  method  was  both  effective  and  presented  the  ideas  in  a 
convenient  form. 


9h. 


The  following  sub-routine  forms  the  sin  8  and  cos  0  of  the 
frame  segments  by  the  use  of  the  basic  equations  listed  below. 


As  =  Ax2  +  Ay2 
n     n    J  n 


Ax  =  x  -x  _ 
n    n  n-1 


n    n  n-l 
sin  ©  =  Ay/As 

cos  ©  =  Ax/As 


y 


(Xn,   Y») 


AV=y»-y»-/ 


Instruction 


Input  to  MRA 


In  MRA 


Transfer 
from  MIA 


gl,  ita  gl+21  This  instruction  transfers  the  address  in  the 
control  element  to  the  address  of  the  instruc- 
tion in  gl+21  and  adds  one.   (See  glossary  of 
instructions ) , 


icr  pbl 


Sets  the  cycle  count  to  the  number  of  segments 
on  the  side  of  the  frame.  This  routine  will  be 
performed  pbl  times 


ica  xl+2+c 
isu  xl+c 
its  tl 
imr  tl 
its  t2 


+x 


n 


-x 


n-l 


+x 


x  -x  _ 
n  n-l 


XA: 


n 


Ax 


n 


Ax 


n 


Ax' 


n 


Instruction 


Input  to  MRA 


In  MRA. 


(continued) 


Transfer 
from  MRA 


9$. 


ica  yl+2+c 

isu  yl+c 

its  t3 

irnr  t3 

iad  t2 


+y 


n 


+y 


n 


-y 


n-l 


Ay. 


n 


♦Ax 


n 


n     n-J- 


2  2 

Ay     +  Ax 
n  n 


Ay 


n 


SP  si    si  is  a  sub-routine  that  extracts  the  square  root 
of  the  number  in  the  MRA  and  returns  the  square 
root  of  the  number  back  into  MRA. 


its     a5+c 


As 


>/Ax     ♦ 


Ay 


ica  tl 

idv  a5+c 

its  a2+2+c 

ica  t3 

idv  a5+c 

its  al+2*c 


+Ax 


n 


-f-Ax 


n 


+Ay. 


n 


^As 


n 


+Ax 


n 


Ax  /As 
n       n 


Ay 


n 


Ay  /As 
J  n!      n 


(cos  9)     =  Ax  /As 
n  n       n 


(sin  9)     =  Ay  /As 
n  n      n 


ict     gl+2        Puts  gl+2   in  the  control  element  and  starts   the 
cycle   over  again „ 


isp 


Address  put  into  this  instruction  by  the  ita 
instruction  in  flad  gl. 


96. 


The  sub-routine  that  follows  forms  the  water  loads  from  the 
origin  to  the  main  decko  In  this  sub-routine  the  water  loads  are  assumed 
to  be  uniformly  distributed  along  the  shell  segment  and  the  center  of 
pressure  is  assumed  to  act  at  the  mid-point  of  the  segment.  After  formation 
of  the  total  load  on  the  segment,  the  load  is  divided  by  2  and  half  of  the 
load  is  placed  at  each  end  of  the  segment.  The  approximation  of  a  con- 
tinuous beam  under  a  trapezoidal  shaped  pressure  load  by  a  uniform  load 
is  only  valid  when  the  shell  segments  are  small  in  size.  As  the  loads 
are  formed,  the  program  converts  the  loads  into  horizontal  and  vertical 
components. 

E 


©*»-  / 


rJ^  Cos  e 


V^iv)r7-  / 


faL 


P  =  Y# 


Instruction 


Input  to  MRA 


In  MRA 


Transfer 

from  MRA 


ita  g2+l8   (See  glossary  of  instructions). 

icr  pb3     Sets  cycle  counter  for  number  of  segments  In  the 
starboard  side  from  origin  to  deck  edge. 


g2,  ics  yl+c 


-J 


n 


-y 


n 


isu  yl+2+c   -y 


n+1 


"W^n  "  "(yn+lV 


97. 


Instruction 


imr  cl 

iad  e2 

imr  c3> 

imr  a5>+c 

imr  cl 

its  tl 

imr  a2+2+c 

its  eii+2+c 

iad  eli+c 

its  el|+c 

ica  tl 

imr  al+2+c 

its  e3+2+c 

iad  e3+c 

its  e3+c 

ict  g2 


Input  to  MRA 

In  MRA                 Transfer 

from  MRA 

(continued) 

0.5 

-(3Wyn)/2  "  ~yav 

♦H(water  head) (ft) 

H-y 
Jav 

*62|       /ft3 

6MH-yay) 

its  (ft) 

6l*a(ft-T„)  =  Fw 

xo.5 

F _/2   (assumes  simple  support) 

V2 

^ (cos  ©)     t 
n+1 

(cos  eWV2  -  (v.Ui 

(Vn*l 

V 
n-1 

V         ♦?     ■  V 
w     -      n        w 
n-1 

(v  ) 

w  n 

+(sin  9) 


n 


(H  )     _ 
w  n-1 


■(V2>n 


(vsin9)     =  (H  ) 

d  n  w  n 


(H  )     _  +  (H  ) 

w  n-1       w  n 


(H   ) 

w  n 


=  H 
w 
n 

(H  ) 
w  n 


cycle  back  to  g2, 


isp 


Address  provided  by  ita  instruction, 


98, 


Formation  of  vertical  forces  is  accomplished  by  the  programs 
listed  below*  Vertical  forces  are  the  algebraic  addition  of  buoyant  forces 
which  act  in  the  upward  direction  and  the  weight  forces  which  act  down- 
ward. Upward  forces  are  taken  as  positive.  Note  that  a  separate  program 

is  required  for  each  side.  The  V  forces  are  symmetrical  and  can  be 

w 

stored  in  50  registers.  The  weight  forces  are  non- symmetrical  and  the 

weight  for  each  of  the  stations  must  be  stored.  In  the  program  for 

forming  the  starboard  vertical  forces,  the  addresses  of  the  V  moves  down 

'  w 

a  column  of  registers  of  increasing  addresses  which  is  the  manner  of 
operation  of  the  el|+c  address.  To  move  the  other  direction  up  a  column 
of  numbers  a  small  Whirlwind  program  is  required  to  reduce  the  addresses 
by  2  for  each  cycle. 


Instruction 


Input  to  MRA 


In  MRA 


Transfer 
from  MRA 


ita    (See  glossary  of  instructions), 
icr  pa2 


glO,  ica  eli+c     +V 


+V 


w 


lad  e5+c         -W, 


its  e5+c 


w 


+V  -WA 
w  t 


=  V. 


V. 


End  of  starboard  side  program. 


ica  eU+c         +V 


its  e5+pb8+pb8 


OUT 


w 


+V 


iad  e5+pb8+pb8    -W 


w 


+v  -wx  =  v., 

w  t    t 


Goes  out  into  Whirlwind  to  change  addresses 
of  weight  loads. 


99. 


Instruction       Input  to  MRA       In  MRA       Transfer 

from  MRA 

(continued) 

ca  glO+i;         iad  e£+pb8+pb8     iad  e!?+pb8+pb8 

su  c?  -2  iad  e£+pb8+pb8-2 

ts  glO+L.  iad  e5+pb8+pb8-2 

td  glO+5  e£+pb8+pb8-2 

IN 

ict  glO+2 

isp 

After  forming  sin  0  and  cos  9  of  each  of  the  segments,  it  was 
necessary  to  find  the  sin  and  cos  of  the  angle  $   formed  by  the  tangent  to 
the  curve  at  each  station  and  the  horizontal  axis,  cos  /  and  sin  $  were 
determined  by  averaging  the  cosines  and  sines  of  the  angles  of  the  chords 
on  either  side  of  the  statione  While  it  was  considered  to  be  more  accurate 
to  find  the  tan  $   or  the  slope  of  the  curve  at  each  station,  the  possibility 
of  tan  $  -   Ay/Ax  becoming  infinity  rendered  this  impossible.  Further,  only 
sin  /  and  cos  <f>   were  necessary  for  the  purposes  of  the  program  and  tan  / 
would  have  required  additional  storage  and  programming.  This  implies  two 
conditions  that  do  not  exist.  Onep  the  ""ate  of  change  of  angles  should 
be  equal  and  secondly,  the 

•   (ei  *  V 

sin  ^ 

does  not  equal 


sin  0_  +  sin  0 


100, 


and 


cos 


(e1  +  o2) 


does  not  equal 


cos  9,  +  cos  9» 


However,  the  resulting  error  will  be  small  as  long  as  the  station  points 
are  judiciously  selected  to  maintain  small  nearly  equal  angle  changes. 
Should  more  accurate  angles  be  necessary  the  programmer  could  utilize 
trigonometric  identities  and  equal  angle  changes. 

An  additional  problem  in  this  sub-routine  was  to  insure  that 
the  centerline  values  of  angles  on  deck  and  bottom  represented  0  and  180 
respectively.  As  was  mentioned  earlier  the  sin  and  cos  of  13£  were  placed 
at  the  corner  station. 


Instruction 

g3,  ita  g3+20 

icr  pb2 

ica  al+2+c 

iad  al+li+c 

imr  cl 


Input  to  MRA 


In  MRA 


(See  glossary  of  terms). 


+sin  9 


Transfer 

f^om  MIA 


1 


♦sin  9, 


xo.5 


+sin  9., 


+sin  9..  +  sin  9? 


sin  9_  +  sin  9? 


sin  / 


its  al+2+c 
ica  a2+2+c 
iad  a2+Li+c 


+cos  9, 


+cos  9, 


sin  $ 


+cos  9, 


cos  9,  +  cos  9_ 
h  2 


101. 


Instruction 


imr  cl 


Input  to  MRA 


In  MRA 


Transfer 
from  MRA 


(continued) 


0.5 


cos  9,  +  cos  0? 


cos 


its     a2+2+c 
ict     g>2 


cos 


The  following  orders  sets  0  and  180  in  the 
center  line  deck  and  keel  and,  also  sets  135 
in  the  corner  angle. 


ica  c6  +0.0 

its  al 

its  al+pbl+pbl 

ics  c3  +1.0 

its  a2 

ica  c3  -1.0 

its  a2+pbl+pbl 

ica  cl6         +0.707 

its  al+pb3+pb3 

ics  cl6         -0.707 

its  a2+pb3+pb3 

isp 


+0.0 


+1.0 


-l.o 


+0.707 


-0.707 


+0.0 
+0.0 

+1.0 
-1.0 
+0.707 
-0.707 


Since  web  frames  are  built  up  sections  it  is  necessary  to 
compute  the  moment  of  inertia  at  each  station  around  the  girth.  The  width 
of  shell  plating  assumed  to  be  effective  in  forming  the  web  beam  is 
thirty  times  the  plate  thickness. 


102, 


Wc 


i. 


E 


tur- 


a 


w 


MA- 


ii 


7" 


3?& 


yf6f. 


The  following  sub-routine  uses  the  relationships  below  in 
forming  the  area,  moment  of  inertia  and  distance  to  the  neutral  axis 
of  the  web  from  the  shell* 


Af  =  «  t 


mf "  (d«  ♦  r)Af 


Jf  ■  (dw  ♦  f)mf 


A  =dt 

W     WW 


m  =  -~^A 
w   2   w 


d      t  d 
T     w       w  \ 

J-  =  .. — p  m  *»  — r — 
w   1.5  w    3 


A  =  30  t 
s      s 


t 
s   2  s 


t 

T        S 

1   ■  r— =r  m 

s   1.5  s 


Ax  =  A  +A  +A 
t    f  w  s 


m,  -  m„+m  +m 
t    f  w  s 


I.  = 


f  w  s 


m, 


^'\ 


NA 


Witt 


103. 


Instruction 

Input  to  MRA 

In  MRA      Transfer 

from  MRA 

gU,  ita  g^2U 

icr  pal 

ica  bl+c 

*% 

*«w 

imr  b2+c 

\ 

\,*\-K 

its  tl 

A 

V 

imr  bl+c 

*V 

w 

imr  d 

*+o.S 

0.5  *■  A,  *  d^  -  «, 

its  t2 

m 
w 

ica  bU+c 

+tf 

+tf 

imr  b3*c 

X  Vf 

*f * vt  -  *f 

its  t3 

A™ 

ica  bl|+c 

♦tf 

+tf 

imr  cl 

xo.5 

0.5  x  tf 

iad  bl+c 

♦d 
w 

«v  ♦  V2 

imr  t3 

XAf 

Af(dv*V2)-mf 

its  ih 

mf 

ica  b£*c 

♦t 
s 

♦t 
s 

imr  b5+c 

At 
s 

<V2 

imr  c2 

X30 

30(t  )2  „  A 

S       8 

its  t$ 

A 

10lu 


Instruction 


imr  b5+c 


imr  cl 


Input  to  KRA 


In  MRA 


(continued) 


** 


xo5 


A  *  t 
s    s 


A  t 
s  s 


m 


Transfer 
from  MRA 


its  t6 

ica  t2 

iad  tli 

isu  t6 

its  t? 

ica  tl 

iad  t3 

iad  t5 

its  d2+c 

ica  t7 


m 


+m 


w 


+m, 


-m 


+A 


w 


♦A. 


+A 


hM. 


m 


W 


m     +  m_ 
w         f 


m  +m  -m     «  M 
w     f     s         t 


K. 


+A 
w 

A  4A  _ 
w     f 

A  +A  +A^,  «  A^ 
s     w     f         t 


►M. 


idv     d2+c 


^A 


A . 


dNcA. 


its     dl+c 
g$-,   ica     t2 

imr     bl+c 

idv     ck 


+m 


w 


xd 


w 


1.5 


♦m 


w 


m  d     = 
w  w 

Ad2 

w  w 


Ad2 

w  w 


w 


lN.A. 


its     tl 


x-t 


105. 


Instruction 

Input  to 

MRA 

In  MRA 

Transfer 

from  MRA 

(continued) 

ica 

bii+c 

*f 

*f 

imr 

cl 

xo<>5 

V2 

iad 

bl+c 

d 
w 

W2 

imr 

tii 

X  mf 

m„(d  +.5tJ 

I       W             1 

-  V  W 

iad 

tl 

+1 

I  +!,. 
w     f 

its 

tl 

(W 

ica 

t6 

+m 
s 

+m 
s 

m        w       t 

idv 

ch 

+1,5 

s         s       s 
1.5  s3       3 

w  t3 

imr 

b5+c 

xt 

S 

s   s        T 
3            s 

iad 

tl 

I   +1- 

W      1 

w     f     s         about  shell 

its 

tl 

Ishell 

ics 

d2+c 

-At 

"At 

imr 

dl+c 

XdN,A. 

"AI  ^.A. 

imr 

dl+c 

XdN.A. 

"AI  4.A. 

iad 

tl 

about  shell 

2 
"Va,   "  Ishell"AI  <XA0 

its 

tl 

"^•A. 

ics 

c3 

-(-1.0) 

+lo0 

idv 

tl 

^H.A. 

Its 

d3+c 

1/I„  . 

Jf 


106, 


When  working  with  beams,  all  forces  and  moments  are  considered 
to  act  along  the  neutral  axis.  To  make  the  program  more  general,  the 
original  coordinate  inputs  were  the  x  and  y   coordinates  of  the  inner 
surface  of  hull  shell.  The  coordinates  of  the  neutral  axis  may  be  formed 
by  offsets  from  the  shell  coordinates 0     The  reader  will  notice  that  such 
a  program  allows  alteration  of  the  web  section  as  desirer1,  to  produce  the 
practical  results.  The  offset  distance  is  cL  . 

dv.A.  Cosf 


<1a/.a.   5'0  ^ 


*N0A.  =  Xs  "  VA.  Sin  *> 

yM.A.  "  ys  +  dN.A.  C0S  *• 

The  following  program  forms  the  coordinates  of  the  neutral  axise 
Since  this  program  is  a  continuation  of  the  above  sub-routine,  the  counter 
instructions  are  still  in  effect. 


Instruction 


g6,  ics  dl+c 


Input  to  MRA 


+d, 


%A. 
idv  cl5         -KL2 

its  tl 


In  MRA 


Transfer 
from  MRA 


"W* 


107. 


Instruction 


imr  al+c 


iad  xl+c 


Input  to  MRA       In  MRA 


Transfer 
from  MRA 


(continued) 


xsin  ^ 


shell 


+  x 


shell 


dN.Aasin  * 
12 


dNaAssin  * 
Xs  "  12     "  XN.A. 


its  xl+c 

ics  tl 

imr  a2+c 

iad  yl+c 


-Ht.a/12 

X.COS  ^ 


shell 


XA. 


+dNcA.A2 


shell 


N.A,  ^  , 


7s   +  iS"^  C°S  *   "  yN#A. 


its  yl+c  y^ 

ict  gU 

isp 

The  moments  at  each  station  are  formed  by  multiplying  the 
summation  of  the  horizontal  and  vertical  components  of  the  loads  by  the 
vertical  and  horizontal  projections  of  the  segment,, 


l/\ n  ~Xn-i) 


IXn-ij    Yn-i)         H5* 


A  Hn-t 


ZVn-, 


108, 


ZV   (x  -x  _ )  -  }    H  Ay  -y     _ )  ♦  ra 
n-1  n  n-1    /  v  n-1  ''n  ''n-l     n 


0 


The  following  sub-routine  forms  the  moments  at  each  station  due 
to  the  weight  and  water  pressure  loads  on  the  starboard  side.  Note  that 
the  first  fevr  instructions  put  +0.0  into  tl,  t2,  t3.  This  was  necessary 
due  to  their  use  in  forming  the  summations 0  In  addition  when  using  this 
sub-routine  to  form  m  ,  tl|  and  b3  must  be  set  to  zero.  When  using  this  to 
compute  the  bending  moment  after  solving  for  the  redundant,  the  value  of 
the  redundant  iru  must  be  put  into  tk   and  m  /ln  in  b3. 

Instruction       Input  to  MRA       In  MRA       Transfer 

from  MRA 

g7,  ita  g8+35 

icr  pbl 

ica  c6  +0.0  +0.0 

its  tl  +0.0 

its  t2  +0.0 

its  t3  +0.0 

ica  e£+c        +V  _  V  _ 

n-1  n-1 

n-2  n-1 

tl  >^  V 


n-1 

S7  v 


its 

tl 

ics 

xl+2+c 

-(x     ) 

na  n 

n-1 

imr 

tl 

*Z_J 

X 


(-x   ) 

na  n 


n-1 


n  v 


L09o 


Instruction 


its  t3 

ica  e3+c 

iad  t2 

its  t2 

ics  y+2+c 


ict  g7+8 


Input  to  MRA 


In  MRA 


.continue 


d) 


+H 


w 


n-2 


H 


w 


("7N.A.  II 


iad  yl+c         (^N.A^n-1 


n-1 


Transfer 
from  MRA 


(nn) 


+H 


w 


n-1 


H 


w 


n-1 


("7N.A.  N 

"I^HJlA"  (3rN.A.}n.l]  =  ZN 


n-1 

^ 7 


imr 

t2 

x      >     H 
/    v  w 

A   7 

7/ 

H 

\    w 

iad 

t3 

n  n— 1 
n-1 

K>» 

iad 

tli 

V   («  ) 
/    \      n 

(Mn) 

its 

tli 

(„n) 

imr 

d3+2+c 

X  i/i 

V1 

its 

b>2+c 

in  /I 
n 

110. 


The  following  section  of  the  sub-routine  forms  moments  port 
side.  Remember  that  the  coordinates  of  neutral  axis  and  the  horizontal 
loads  are  stored  for  one  side*  The  sub-routine  must  move  up  the 
stored  column  of  numbers  to  form  the  moments. 


Instruction       Input  to  MRA       In  MRA       Transfer 

from  MRA 

icr  pal 


g8,  ica     e!>+pal 

i-pal+c        +V 
^                    n 

♦V 
n 

iad     tl 

n-2 

X- 

n-1 
ZV+Vn 

n-1 

its     tl 

n-1 

ica  xl+pb2+pb2 

+aW 

+Xn*l 

isu  xl+pbl+pbl 

-xa 

X       --X       = 

n+1     n 

r   A 
X 

luir     tl 

n-1 

n-1 
A     >     V 

-  m 
n 

its     t3 

(Mn) 

ica     e3*pbl+pbl 

n-1 

+Hn-1 

iad     t2 

n-2 
XH-2 

n-2                            n-1 

X^-vi  -Z 

Ill, 


Instruction 

Input  to  MRA       In  MRA 
(continued) 

Transfer 
from  MRA 

its  t2 

n-1 
/  .i.  n 

ica  yl+pb2+pb2 

^+yN.A.  n+1        ^+7N9A. 

n+1 

isu  yl+pbl+pbl    (-y^  A  ) 


imr  t2 


n-1 


h 


[^N.A/n+l  "  ^N.AjnJ 


n-1 


iad  t3 


n  v 


(m  )_  +  (mT,)H  »  m 
n  v     n  n    n 


n-1 

7 


iad  ml+pbl+pbl+c 


m 


n 


n-1 


n 


m  +  m  = 
n 


m 


imr  d3+pbl+pbl     l/I 


n 


m/I 


OUT      The  program  now  goes  out  into  Whirlwind  to 
change  the  addresses  for  each  cycle „ 


ca  g8+3 
td  g8+U 


ica  xl+pb2+pb2      ica  xl+pb2+pb2 

Puts  the  address  of  g8+3  into  the  address  of 
g8+i|  now  reads  isu  xl+pb2+pb2  e 


su  c? 
ts  g8+3 
ca  g8+7 


-  (+2W.W0 )  ica  xl  >pb2+pb2-2 

ica  xl+pb2+pb2-2 
ica  e3+pbl+pbl     ica  e3+pbl+pbl 
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Instruction       Input  to  MRA       In  MRA       Transfer 

from  KRA 

(continued) 

su  c7  -(+2  W.W.)         ica  e3+pbl+pbl-2 

ts  g8+7  ica  e3+pbl+pbl-2 

ca  g8+10  ica  yl+pb2+pb2      ica  yl+pb2+pb2 

td  g8+ll     (Instruction  in  g8+ll  is  now  isu  yl+pb2+pb2). 

su  c7  -(+2  W.W.)         ica  yl+pb2+pb2-2 

ts  g8+10  ica  yl+pb2+pb2-2 

ca  g8+l5  +imr  d3+pbl+pbl     +imr  d3+pbl+pbl 

su  c7  -(+2  W.W.)         +imr  d3+pbl+pbl-2 

ts  g8+l5  inr  d3+pbl+pbl-2 

IN 
ict  g8 
isp 

Since  all  moments  due  to  applied  loads  and  forces  had  been 
formed,  it  was   possible  to  perform  the  integrations  of  moments  to  obtain 
the  slopes  and  deflections  previously  discussed  in  the  procedure.  At 
this  point  in  the  programming  the  redundants  m ,p  ,  and  q^  must  be  con- 
sidered. mn  will  appear  as  a  constant  value  of  moment  at  each  station 
throughout  the  girth.  The  moments  due  to  p  and  cu  depend  upon  the 
value  of  y  and  x  at  each  station.  Since  the  web  frame  was  assumed  to 


redundant    On    On    0 
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be  symmetrical 5  the  slope  changes  and  deflections  around  the  entire  girth 
caused  by  the  redundants  will  be  double  the  slopes  and  deflections  on  one 
side0  Fo:  convenience  in  integration,  the  redundants  will  be  integrated 
around  the  starboard  side,, 

The  method  of  integration  used  was  the  trapezoidal  rule  which 
implies  linear  slope  of  the  moment  curve  between  stations  on  the  moment 
curve.  The  area  under  the  curve 


y^fo  +(M)   "l^n 
/  \  [J'n       KI   n+ll  2 

n=l 


[& 


+  (T}2  T  +  (T}2  +  (T}3 


1^2 
3j  2 


q/j,  2   c/5,  3  c/Ss     *       &$j     *      dSj-   8  ct  S 


c/s-,8 


As  can  be  seen  from  the  figure,  the  trapezoidal  rule  will  cause  small 
error  providing  the  curve  is  smooth  and  stations  are  closely  spaced* 
Despite  the  small  error,  the  above  method  of  integration  was  chosen 
because  it  is  readily  applicable  to  computer  programming  and  because 
random  station  spacing  was  necessary  to  approximate  a  curved  figure 
by  a  series  of  straight  lines „ 


lliu 

In  the  process  of  "debugging"  the  program  several  improvements 
in  the  sequence  of  operations  for  the  formation  of  the  integrals  were 
made.  The  results  of  these  improvements  and  errors  in  addresses  of  the 
initial  program  for  the  integration  procedures  alters  the  procedure  as 
written.  The  techniques  used  and  explained  in  the  follox-iing  pages  are 
correct. 

The  major  improvement  was  the  inclusion  of  the  integration 
routines  under  one  sub-routine.  All  the  integrals  for  nu,  p_,  c^  and 
m  /I  were  performed  under  one  cycle  counter.  Each  of  the  stanchion 
redundant  forces  R  and  S  Here  under  a  separate  cycle.  The  routine 


I' 


P0  )  I  ds 


was  not  performed.   To  obtain  the  coefficient  of  p  ,  the  coefficient  of 


mo 


which  is 


P 


•J  ds 


was  used.  The  group  of  instructions  included  under  flad  g°  were  added 
to  insure  that  zeros  were  initially  placed  in  the  registers  that  were 
to  hold  the  matrix.  Also,  the  addresses  for  forming  the  port  side 
integrals  were  found  to  be  in  error  by  two  registers. 

The  instructions  listed  below  are  correct  and  reolace  the 
integration  routines  listed  under  flads  gll  through  g22  of  the  detailed 
description.  Unfortunately,  these  alterations  were  not  made  in  sufficient 
time  to  allow  the  correction  of  the  detailed  explanation.  These  instruc- 
tions should  be  read  from  left  to  right. 
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g9,  ita  g22 

ict  g2+2 

imr  a5+c 

imr  xl+c 

imr  xl+2+c 

iad  dl+2it 

its  tl 

iad  tl 

gill,  ica  tl 

imr  yl+2+c 

its  dl+12 

imr  cl 

imr  b3+c 

iad  t3 

its  dl+52 

ica  xl+2+c 

imr  cl 

iad  b3+pal+pal+c 

its  dl+50 

icayl+pb2+pb2 

imr  cl 

imrb  3+pbl+pbl+c 

iad  tl 

its  dl+5i| 

ts  gl8+2 

ts  gl9 

ts  gl9+10 

ict  gll 

icaxl+pbli+c 

imr  cl 

its  dl+56 

imra5+pbl|+c 

iad  dl+56 

its  tl 


icr  30 

icr  pbl 
iad  dl 
imr  d3+c 
imr  d3+2+c 
its  dl+2U 
ica  yl+2+c 
imr  a5+c 
imr  yl+c 
iad  tl 

gl5,  ica  b3+c 
iad  dl+^O 
its  t3 
imr  a£+c 
gl7,  ica  xl+c 
imr  b3+2+c 
iad  dl+5U 
imr  a5+pb2+pb2 


ica  c6 

gll,  ica  d3+c 
its  dl 
its  tl 
iad  tl 

gl3,  ica  yl+c 
imr  d3+2+c 
iad  dl+2 
its  tl 
imr  a5+c 
iad  b>2+c 
its  dl+50 
ica  yl+2+c 
imr  cl 
imr  b3+c 
iad  t3 
its  dl+5h 
imr  cl 


gl9, icayl+pbl+pbl imrb 3+pbl+pbl+c 
imrb3+pal+pal+c  iad  tl 


iad  dl+52 
its  tl 

imra5+pb2+pb2 
OUT 

ts  gl9+6 
su  c? 
su  c7 
ica  dl+2 

isu  xl+pbli 
its  tl 

icaxl+pblj+2+c 
imr  cl 
its  dl+56 
imryl+pbli+c 


its  dl+52 
icsxl+pb2+pb2 
imr  cl 

ca  gl8+2 
ts  gl9+l6 
ts  gl9+3 
ts  gl9+13 
its  dl+10 
its  t3 
imrb3+pbL+c 
isuxl+pblj 
its  t2 
ica.  tl 
iad  dl+l6 


its  dl+c 

iad  d3+2+c 

gl2,  icaxL+c 

ica  rl+2+c 

imr  a5+c 

imr  d3+c 

its  t2 

its  dl+2 

ica  t2 

iad  dl+12 

imr  a5>+c 

gl6,  ica  yl+c 

imr  b>2+c 

iad  dl+52 

its  t3 

imr  a5+c 

gl8, icab 3+pbl+pbl+c 

iad  dl+50 

its  tl 

imra5+pb2+pb2 

icsxl+pbl+pbl 

imrb3+pal+pal+c 

iad  dl+5a 

su     c7 

ca  gl9+3 

ca  gl9+13 

IN 

g20,icr  pa5 
irara^+pbb+c 
iad  dl+56 
its  tk 

imrb3+pbl|+2+c 
imrd3+pbli+c 
its  dl+16 
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ica  t2 
iad  dl+16 
iad  dl+26 
iad  dl+26 
iad  dl+36 
iad  dl+36 
iad  dl+6 
its  dl+30 
its  dl+3li 
its  t3 

imrb  3+pb6+pb6 
isuxL+pb5 
its  t2 
OUT 

ts  g21+l6 
its  tl 
ica  t2 
iad  dl+18 
iad  dl+28 
iad  dl+28 
iad  dl+l|8 
iad  dl+l;8 
iad  dl+8 
its  dl+l;0 
its  dl+l|i+ 


imrd3+pbli+2+c 
its  dl+16 
its  dl+26 
its  dl+26 
its  dl+36 
its  dl+36 
its  dl+6 
ica  dl+16 
g21,    icr  pa6 
imra^+pb^+c 
iad  dl+58 
its  tU 

imrb  3+pb6+pb6-2 
ca  g21+l6 
IN 

imiyl+pb5+c 
imrd3+pb5+2+c 
its  dl+18 
its  dl+28 
its  dl+28 
its  dl+US 
its  dl+JU8 
its  dl+8 
ica  dl+18 
g22,   ispO 


its  t2 
ica  tl 
ica  t2 
ica  t3 
ica  tU 
ica  tl 
ict  g20+l 
its  dl+32 
icaxl+pb5>+c 
imr  cl 
its  dl+58 
imra5+pb5+c 
iad  dl+58 
ts  g21+7 
ica  tl 
iad  dl+18 
its  t2 
ics  tl 
ics  t2 
ica  t3 
ica  tl| 
ica  tl 
ict  g21+l 
itx  dl+U2 


imiyl+pbl4+2+c 

imrxl+pbi|+c 

imrxl+pblt+2+c 

imr  tl 

imr  t2 

iad  t2 

ica  dl+6 

ica  dl+26 

isuxl+pb5> 

its  tl 

icaxl+pb£+2+c 

imr  cl 

its  dl+58 

su  c? 

imrd3+pb5+c 

its  dl+18 

imiyl+pb5+2+c 

imrxL+pb5+c 

imrxl+pb5>+2+c 

imr  tl 

imr  t2 

iad  t2 

ica  dl+8 

ica  dl+28 
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The  results  of  each  integration  are  stored  to  form  the  matrix 
as  shown  by  Fig.  VI.  The  following  sub-routine  performs  the  integration 
of  the  change  of  slope  around  the  girth  of  the  web  frame  due  to  redundant 
moment  vu0 


AG 


0  A  A 

0         0  o   r 


,    ds 
n     n+1 


~*   n    n+1 


Instruction 

gll,  ita  gl5+9 
icr  pbl 

ica  d3+c 
iad  d3+2+c 

imr  a5>+c 


Input  to  HRA 


In  MIA 


Transfer 
from  MRA 


(See  glossary  of  instructions). 

(Sets  cycle  counter  equal  to  the  number  of 
stations  on  a  side  minus  l). 


+1/1  - 

'  n-1 

+1/1 


n 


x  ds 


n 


+1/1 

n-1 

+1/1  +  1/1  _ 

1    n    '  n-1 

(1/1  +  1/1  _)ds 
'    n     n-1   n 


iad  dl 


n-1  n 

+  Z(l/In  +  lAn-l)dSn-l  X  (V*»  +  1Art)d«n 


its  dl 


n 


(1/1  +  1/1  _)ds 

'    n    '  n-1   n 


The  instructions  below  find  the  deflection  in  the  y-direction 
due  to  the  moment  caused  by  the  shear  redundant  q  . 


Ay 


0  A  n-1 

=  J  x(q()x)ds  »  2  J   x(qQx)ds  =  2qQ  \ 
0  0  1 


n-1       n 


n-1 


n 


ds 
2 
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Instruction 

Input  to  MRA 

In  MRA 

Transfer 

from  MRA 

gl2,  ica  xl+c 

+x  1 

n-1 

+x  _ 

n-1 

irar  xl+c 

*  Vi 

<v/ 

imr  d3+c 

'  ^n-l 

(Sa^ 

-1 

its  tl 

K-l^n-l 

ica  xl+2+e 

♦x 

n 

X 

n 

imr  xl+2+c 

X  X 

n 

2 

X 

n 

imr  d3+2+c 

x  1/1 

n 

(x  )2/I 
n  '  n 

iad  tl 

♦(*  -,)2A  -, 

n-1  '  n-1 

K^ ♦ 

'Vl^l 

imr  a£+c 

x  ds 

n 

^X  + 

(xn-l)2/Jn-l  ds 

n-1  n 

7 .         .0    .  .  .9    . 

ds 


iad     dl+2l[  +  N     ds     -    (x  /I  )2  +   (x     -  )2/l     i      yU  )2/l  +(x     .  )2/l     , 

/     v     n-1       n'    n  n-1     '    n-1    /   Q  n     '    n       n-1     7    n-1 

n 

its     dl+2ii  >      (x  )2/l  +(x     n)2/l     .    ds 

/   \   n  '  n   n-1  /  n-1   n 

The  deflection  in  the  x-direction  due  to  mn  was  found  by  per- 
forming the  integration 

0  A  n-1 


m         /-*   m 

Ax  -  J   y^dB-2    7^dsB2iL 


/  A ds  -  2-oZ  (^.i/w+vv]  ^ 


i 

0  0 


The  instructions  starting  with  g!3,  perform  this  integration. 
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Instruction 

Input  to  MRA 

gl3,  ica  yl+c 

%i-l 

imr  d3+c 

x  1/1  n 

'    n-1 

its  tl 

ica  yl+2+c 

+y 

imr  d3+2+c 

x  1/1 

/  n 

iad  tl 

*  VAi-l 

imr  a5+c 

x  dsn 

^L, 

iad  dl+2 

V  W^n- 

In  MRA 


+7 


n-1 


vA-] 


+7. 


n 


y  A 

Jv!    n 


Transfer 
from  MRA 


^n-l'n-] 


(y  A  +y  -.A  J 

Jn'  n  Vl'  n-1 

(y  A  +y  ,/l  n)ds 
n  n  n-1'  n-1   n 


n 


(y  A  +y  -i/y  n  )ds 

Jn'  n  •yn-l/,'n-l/  n 


its  dl+2 


n 


(y  A  +y  ,A  J 

n  n  n-1  n-1 


ds 


n 


The  following  section  of  the  sub-routine  performs  the  integration 
of  the  moment  due  to  p  o  p  v  around  the  girth. 


AG 


P(f 


ds 


y~    y  .-i  ds 

n+l  \  n 


0 


2Por^ds  =  2p0y^+^) 

-/  ^ >  n    n+l 
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Instruction 

Input  to 

MRA 

In  MRA 

gill,  ica  yl+c 

+yn-l 

+yn-l 

imr  d3+c 

X1'1*.! 

yn-l/:Ln.l 

its  tl 

ica  yl+2+c 

+yn 

+yn 

imr  d3+2+c 

xi/ln 

7   A 
'n'  n 

its  t2 

iad  tl 

*7     -,A 
J  n-1'  n- 

-]. 

y  /I  +y 

Jn'    n  Jn- 

imr  a5+c 

X  ds~ 

(y yi+y 

Transfer 
from  MRA 


WVl 


j  A 

Jn'   n 


n 


iad     dl+10 


n-1 


n 


(y     n/l     n+y     ./I     Jds     _       >      (y  /I  +y     _A     Jds 
_yn-l'    n-1  J  n-2'    n-2       n-1  ^_ kVJn'    n  'n-l'    n-1       n 


its     dl+10 


n 


(y  /I  +y    _/l     _)ds 
'n'    n  Jn-r    n-1       n 


The  x-deflection  due  to  the  axial  redundant  force  pn  was 


V 


ds  = 


2pQ  J  -^  ds  .  2P() 
0 


n=:Lry2  (7  -,)2lds 

'n  .  Jn-1     n 


I 
L_n 


n-1 


2  * 


The  instructions  given  in  gl5,  performs  the  p  moment  x-deflection.  Note 
that  the  values  of  y  ../I  _  and  y  /I  formed  in  glii,  and  stored  in  temporary 
registers  is  used  in  gl£,  thus  saving  operations. 
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Instruction 

Input  to  MRA 

gl5,      ica     tl 

+yn-l/:En-l 

imr    yl+c 

xVi 

its     tl 

ica     t2 

+y  A 

•W    n 

imr    yl+2+c 

xyn 

ica     tl 

+W\-1 

imr     a£+c 

Xdsn 

iad    dl+12 

n-1 

+N       (y2/l  +y: 
/    .     n    n  Ji 

In  MRA 


^n-A,-! 


Transfer 
from  MRA 


^-l^n-l 


/A 

Jn'  n 

(y2/l  +y2  _/l   _)ds 
Jn'  n  ^n-r  n-1   n 


n 


its  dl+12 


n 


(y2/l  +y2  _/l  _)ds 
wn  n  n-1'  n-1   n 


ict  gll+2 
isp 


(Cycles  back  to  gll+2). 

(Address  set  by  ita  instruction  at  gll). 


Due  to  symmetry  of  the  web  frame  structure,  the  following  re- 
dundant integrals  are  zero.  Changes  in  slope  due  to  cu. 


A9  = 


qo 


ds  = 


(+x)  -ydS  + 


(+x)-j-ds 


0 


A 


qQ  J  ^ds  + 
0 


u  • 


(-x). 


Is  =  0 
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Vertical  deflection  due  to  rau  and  p„. 


Q1  A       Q1  A       01 


Ay  =    -^ds  =  n^   [  -fds  +  (  ^ds  =  ra0  [  ~ids  +  )  ^^ds  -  ° 

Jo  ^L 


"4 


'o 


0'  /    x         A        o' 


Ay  =  J    -4jl-  ds  =  pQ  J  ^ds  +  J      (-x)^ds  -  pQ  0  . 
vO  Jo  Jk 


Horizontal  deflection  due  to  <l, 


V 


r  y(qQx) 

Ax  = 


-J    -r-  ds  -  %  \  -?ds  +  /     (-x)ids  -  <*o  ° 

v6  Jo  Jk 


The  sub-routines  below  perform  the  integration  of  moments  caused 
by  water  pressure  and  weight  loads.  Note  that  two  programs  are  necessary 
since  h9   stations  are  used  in  the  integration  but  only  25  x  and  y  coordinates 
had  been  stored.  The  following  program  performs  the  integration  of  the 
moments  along  the  starboard  side  to  find  the  change  in  slope  for  one-half 
the  web. 

A  ,  2lt 

a^  rij  \     (  n  .      n-l\ds 

Ae  .    I    _ds  .    \     (_+  ) 

/0  ^     n         n-1 
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Instruction 

gl6,  ica  b3+c 

iad  b3+2+c 

imr  a5>+c 

imr  cl 

iad  dl+50 


Input  to  MRA 


+m  _ /I  . 
n-17  n-1 


+m  /I 
n  n 


X   ds 
A   n 

X  1/2 
n-1 


In  MRA 


+m   /I 
n-1  n-1 


Transfer 
from  MRA 


(m  /I  +m   /I   ) 
n  n  n-1  n-1 

(m  /I  +m  _/I  n  )ds 
ir  n  n-1'  n-1   n 

(m  /I  +m  _/l  _)ds  /2 
n  n  n-1  n-1   n 


n 


+  >   (m  /I  +m  _/l  _)ds  /2   >   (m  /I  +m  , /I  -  )ds  /2 
/  .»   n  n  n-17  n-1   n7   X  ^  n7  n  n-1  n-1   n 


its  dl+50 


n 


(m  /I  +m  ./I.  n)ds  /2 
n7  n  n-17  n-1   n 


The  segment  of  the  sub- routine  that  follows  forms  the  integration 

of  y  m  /I  to  give  the  deflections  in  the  x-direction  around  the  starboard 
n 

half  of  the  web. 


Instruction 

gl7,  ica  yl+c 

imr  b3+c 

its  t3 

ica  yl+2+c 

imr  b3+2+c 

iad  t3 

imr  a5+c 


Input  to  MRA 

+yn-l 
X  VA-1 


+7 


n 


xx   n  n 

+y  n  m  _/l  _ 
J  n-1  n-17  n-1 

X  ds 
^   n 


In  -MRA 


+7 


n-1 


y    m.  „/l  t 

17  n-1  n-J.7  n-1 


+7_ 


n 


y  m  /I 
n  n  n 


Transfer 
from  MRA 


y  .,  m  _/l  _ 
* n-1  n-17  n-1 


y  m/l+7  t  m  ../I  ., 
n  n7  n  'n-l  n-17  n-1 

(y  m  /I  +y    m  ../l  ..  )ds 
J  n  n7  n  Jn-1  n-17  n-1   i 


12lu 


Instruction 

Input  to  MM           In  MRA          Transfer 

from  MRA 

(continued) 

imr  cl 

x  1/2                 (y  m /I +y  _  m      /I   )d 
'                   n  n n  •'n-l  n-1'  n-l 

n-1 

iad  dl+52 

+  >   (y  in  /I  +y  -in  _/l  _  )ds  /2 
/  x   n  n'  n  Jn-1  n-17  n-1   n 

n 

7 


(7  m  /I  +y       m       /I     -,)ds/2 
j,     n  n    n     n-1  n-1     n-1       n 


n 


its     dl+52  >     (y  m  /I  +y     _ir     n/l     _  )ds 

/  s.  n  n'  n  ''n-l  n-1'  n-1   n 


To  obtain  the  deflections  in  the  y-direction  on  the  starboard 

side  due  to  applied  loads,  the  integral  of  x  m  /I  was  formed. 

n 


Ay  , 

A         n 

m  nx 
/  n-1  n- 

Jn-1 

-1 

m  x  ds 
n  n}  n 

I  ;  2 
n 

/"*  Xlll      \   / 

/^s  -Z 

0          1 

Instruction 

Input  to  MRA 

In  MRA 

Transfer 
from  MRA 

glS,  ica  xl+c 

+x  ., 
n-1 

+x  n 
n-1 

imr  b3+c 

Y  m  ., /I  n 
*  n-1'  n-1 

x  ..m  t/I  t 
n-1  n-1'  n-1 

its  t3 

X 

n- 

t  m  n  A  1 
-1  n-1'  n-1 

ica  xl+2+c 

+x 
n 

+x 

n 

imr  b3+2+c 

x  m  /I 
n  m 

x  m  /I 
n  n  n 

iad  t3 

+  x  nm   /I 

<x. 

i™ 

/I    +  X 

m  / 

I  ) 

-25, 


12?,, 


Instruction 


imr  a5+c 
imr  cl 


Input  to  MRA 


In  MRA 


xransier 
from  MRA 


(continued) 


*dSn 
x  1/2 


(x  _m  n/I  ,,+x  m  /I  )ds 
n-1  n-1'  n-1  n  n'  n   n 

(x  -.m  _/l  .,+x  m  /I  )ds  /2 
n-1  n-1  n-1  n  n  n   n 


n-1 


iad  dl+5ii     +  7   (x  -.m  _/l  n+x  m  /I  )ds  /2 
/  v   n-1  n-1'  n-1  n  n  n   n 


its  dl+5U 


n 


/  (x  _m  ../I  t+x  m  /I  )ds  /2 
Z i  n-1  n-1  n-1  n  n  n'  n 


ict  gl6+2 


(Cycles  back  to  gl6+2). 


To  perform  the  integrations  on  the  port  side,  the  address  changes 
if  the  x  coordinates,  y  coordinates  and  girth  segments  ds  must  move  up  a 
column  of  registers „  This  is  accomplished  by  subtracting  2  from  the  address 
for  each  cycle  performed. 

Formation  of  the  change  in  slope  due  to  applied  loads  (port  side) 


A9 


0' 


n 


ds  = 


A 


25 


■y     m         m  _  ds 

*  I    I  J  2 

n    n-1 


126, 


Instruction 


Input  to  MRA 


In  MRA 


Transfer 
from  MRA 


icr  pbl       (Sets  cycle  counter  to  number  of  points  in  a  side) 


gl9,  ica  b3+pal+pal+c  +mn..]/In_i 


iad  b3+2+pal+pal+c   +m  /I 

n  n 


m  ,/l  - 
n-1'  n-1 


imr  a5+pbl+pbl 
imr  cl 


x  ds 
A   n 

Kl/2 


(m  /I  +m  ./I  .) 
n  n  n-1'  n-1 

n  n  n-l  n-_L   n 

(m  /I  +m   /I  Jds/2 
n  n  n-1  n-1   n 


25+n-l 


25+n 


iad  d3+50      +>   (m  /I  +m  _/l  Jds  /2   >  (rn  /I  +m  _/I  Jds  /2 
X  >   n  n  n-1'  n-1   n'    /  v,  n  n  n-1'  n-1   n 


1 
25+n 


its  d3+50 


">  '(m  /I  +ni  -/I  -Jds  /2 
X  x.  n  n  n-1  n-1   n 


Formation  of  horizontal  deflection  port  side 


0         25+n 

/m  y      ^- 7    m  _y  „   m  y  ds 

A  25 


Instruction 


Input  to  MRA 


g20,  ica  yl+pbl+pbl  yn_1 

imr     b3+pal+pal+c  x  m     ../I     _ 

^       ^  n-1     n-1 

its  tl 


In  MRA 


7 


n-1 


y  ^m   /I  , 
Jn-1  n-17  n-1 


Transfer 
from  MRA 


y  .m   /I 
n-1  n-1'  n-1 


ica  yl+pbl+pbl-2      +y 


n 


+y 


n 


127. 


Instruction 


Input  to  MRA 


In  MRA 


Trars "Qr 
from  MRA 


(continued) 


imr    b3+pal+pal+2+c        %  m  /I 


y  m  /I 
Jn  n7    n 


iad 

tl 

+m    ny       /I     n 
n-i  n-1     n  l 

imr 

a£+pbl+pbl 

xds 
n 

imr 

cl 

Xl/2 
25+n-l 

-^ y 

iad 

dl+52 

>       (y  m  /I  +y       m 
/     .    'n  n     n  Jn-1  n- 

y  m  /I  +y       m       /I 
n  n     n     n-1  n-1     n-1 

(y  m  /I  +7     -.m    _/I     .  )ds 
n  n     n     n**l  n-1     n-1       n 

(y  m  /I  +y     _m       /I     _  )ds  /2 
nn     n  ^n-l  n-l     n-1       n 


25+n 


(y  m  /I  +y     _m       /I       )ds  /2 
n  n     n  J n-1  n-1'    n-1       n' 


its     dl+52 


Formation  of  vertical  deflection  port  side  due  to  applied 


loads. 


0' 

/xm 

A 


25+n 

-^   m    _x    _       myds 


n-1 


I     '   2 
n 


Instruction 


ica  xl+pbl+pbl 


Input  to  MRA 


+x 


n-1 


In  MRA 


+x 


Transfer 
from  MRA 


imr  b3+pbl+pbl+c      +m  -,/I  1 
its  tl 


n-1 

x  nm  _/I  _ 
n-1  n-1'  n-1 


x  ,m   /I 
n-1  n-1'  n-1 


128* 


Instruction 


imr  b3+pal+pal+2+c 


iad  tl 


mr  cl 


iad  dl+5k 


Input  to  MRA 


In  MRA 


Transfer 
from  MIA 


(continued) 


ica  xl+pbl+pbl-2      +x. 


n 


V  Itl     /I 

A    n    n 


+x 
n 

x  m  /i 

n  n     n 


+x     nm       /I 

n-1  n-1     n-1 


imr     a^+pbl+pbl+c  *  ds 


n 
xl/2 


(y  m  /I  +y     _  m     _  /I     _  ) 
Jn  n'    n  Jn-1  n-1     n-1 

(y  m  /I  +y     _m     _/I     ,  )ds 
Jn  n'    n  Jn-1  n-1     n-1       n 

(y  m  /I  +y    ,m    _/I     _  )ds  /2 
''n  n'    n  'n-l  n-1'    n-17     n 


2g+n-l 

>      (x  m  /I  +x     _m       /I     _)ds  /2 

X     „      n  n'    n     n-1  n-l     n-1       n 


25+n 


(x  m  /I  +x     _m     _/l     _  )ds  /2 
n  n    n     n-1  n-1'    n-1       n 


its     dl-»£l| 


2^+n 


(x  ra  /I  +x     ,m    t/I     _  )ds  / 

n  n    n     n-1  n-1'    n-1       n 


n 


The  following  program  changes  the  addresses  of  the  instructions 
in  gl°  and  g20  that  utilize  x,  y  and  ds.  The  new   addresses  must  move  up 
a  column  of  registers. 


OUT    (Program  goes  into  Whirlwind  mode  to  change  addresses). 


129. 


Input  to  MRA 

imr  a5+pbl+pbl 
-2  W.W. 


-2  W.W. 


-2  W.W. 


In  MRA 

imr  a5>+pbl+pbl 
imr  a5>+pbl+pbl-2 


Transfer 
from  MRA 


ica  yl+pbl+pbl-2       ica  yl+pbl+pbl-2 


ica  yl+pbl+pbl-li 


Instruction 

ca  gl9+2 
su  c7 
ts  gl9+2 
ts  g20+6 
ts  g20+17 
ca  g20+3 
ts  g20 
su  c? 
ts  g20+3 
ca  g20+lU 
ts  g20+ll 
su  c7 
ts  g20+lli 

IN 

ict  gl9 
isp 

The  above  integrations  provide  the  three  equations  for  solving 
an  unsupported  ring  frame  of  three  degrees  of  indeterminancy.  To  solve  the 
subject  problem,  the  equations 


ica  xl+pbl+pbl-2       ica  xl+pbl+pbl-2 


ica  xl+pbl+pbl-lj 


imr  a5+pbl+pbl-2 
imr  a5+pbl+pbl~2 
imr  a5>+pbl+pbl-2 

ica  yl+pbl+pbl-2 

ica  yl+pbl+pbl-li 

ica  xl+pbl+pbl-2 

ica  xl+pbl+pbl-lj 


(Puts  computer  back  in  interpreted  mode). 


/ 


2  m(x-x  )ds 

a 


,E 


and 


h   m(x-x,  )ds 


130. 


must  be  formed  and  solved  simultaneously  with  the  three  previously  formed. 
The  designations  E  9   B   B  and  B,  are  the  end  points  of  the  two  stanchions, 

The  sub-routine  that  follows  determines  the  vertical  deflection 
between  B  and  B  of  the  starboard 


B        .      B 
2  m__(x-x_)    ^ — r   m_(x_-x_)   m  n  (x  ,-x  ) 


Ay  « 


E. 


n    a 


■ds 


a. 


a  .   n-1  n-1  a  -\ds 

I  _     }2 
n-1 


B 


2  r 


■I 


B. 


i  (x  -x  )  ,    m  _  (x  _-x  )ds 

n  n  a  ds    n-1  n-1  a   n 

I     2  I  „      2 


n 


n-1 


Instruction 

Input  to  MRA 

ita  g21+ia 

g21,  icr  pa£ 

ica  xl+pbl+c 

n-1 

isu  xl+pbli 

-x 
a 

its  t3 

imr  a5+pnl;+c 

X  ds 

imr  cl 

xl/2 

its'  tl 

imr  b3+pbii+c 

x  (m/1)^ 

n-2 

In  MRA 


iad  dl+56 


x 


n-1 


(x   -x  ) 
n-l  a 


(x  _-x  )ds 
n-1  a 


Transfer 
from  MRA 


n-l  a 


n-1 


(x  ,-x  )ds/2 

n-1  a   ' 

(x   -x  )ds/2 
n-1  a 

un_1-Xa)(m/Dn_1dsn/2  rt 


i 


B. 


Instruction 


its  dl+56 


Input  to  MRA  In  MRA 

(continued) 


131. 


Transfer 
from  MRA 


n-1 

y    (x  -,-x  )(m/l)  -.ds  /2 
Z. ».  n-1  a     n-1  n 


B. 


ica  xl+pbl|+2  +x. 


n 


+x 


n 


isu  xl+pbli 

-X 

a 

(x  -x J 
n     a 

its     tli 

(x-x  ) 
n     a 

imr  a5+pbi|+c 

x  ds 

(x  -x  )ds 
nan 

imr     cl 

Xl/2 

(x  -x  )ds/2 
n     a 

its     t2 

<W2(xn-xa) 

ijnr  b3+pbl+pbl+2 

+c 

x(m/l) 
n 

(m/l)    (x  -x  )ds/2 
n     n     a 

B1+n-l 

B  +n 

iad  dl+56 

-7 

(m/1)    (x  ■ 
_s.              n     n 

-xa)ds/2 

J 

-7- 

(m/1)    (x  -x  )ds/2 
-^            n     n     a 

B. 


B. 


its  dl+56 


B1+n 

F 


(m/1)  (x  -x  )ds/2 

n  n  a   ' 


The  following  instructions  perform  the  integration 


B 


2  y(x-x  ) 


ds  =  p 


7. 


(x  -x  )  ds    y  _  (x   -x  )  ds 

n  n  a    n  ,  J  n-1  n-1  a    n 


n 


n-1 


132. 


Instruction 

Input  to 

MRA 

In  MRA                        Transfer 

from  MRA 

ica     tl 

(xn-l"Xa 

)ds/2 

(x      -x  )ds  /2 

n-1     a       n 

imr  d3+pbli+c 

*  ^n-l 

(x       -x  /I     Jds  /2 
n-1     a     n-1       n 

its     tl 

(x     ,-x  /I     ,  )ds  /2 
n-1     a!    n-1       n' 

imr  yl+pbU+c 

xVi 

y     .(*       -x  )/l     ,ds  /2 
n-1    n-1     a       n~l     n 

B  +n- 

•2 

B1+n-l 

iad     dl+l6 

•>: 

^n-l^n-l 

-x  )/I 
a  '   n« 

ids  /2 
-1     n 

>     7     _(x     .,-x  )/l     _ds  /2 
/     ^'n-1     n-1     a  '    n-1     n' 

B. 


its     dl+16 

ica     t2 

imr  d3+pbl;+2+c 

its     t2 

imr  yl+pbli+2+c 

B^n-1 

iad     dl+16 

*> : 

(x  -x  )ds  /2 

nan 


(1/D 


n 


xy. 


n 


B1+n-l 


B. 


(x  -x  )ds  /2 
nan 

(x  -x  )/l  ds  /2 
n     an     n 


Z7     ,  (x     ,-x  )/l     ,ds  /2 
^n-l     n-1     a  '    n-1     n' 


(x  -x  )/l  ds  /2 
n     an     n 


(x  -x  )7  /i  ds  /2 

n     an     n     n 


B  +n 

T\~"7       (x  -x  )/l       ds  /2       >      y  (x  -x  )/l  ds  /2 
X    .    n-1     n     a       n-1     n         z. »/n     n     a  '    n     n 

B,  B„ 


B  +n 


its     dl+16 


">      7  (x  -x  )/I  ds  /2 

Z *.    n     n     a       n     n 

B. 


133. 

The  vertical  deflection  between  the  ends  of  the  starboard  stanchion 
due  to  the  redundant  cr  force  was  found  by 


B 


B, 


ax  .,  (x  _-x  )  ds    y  (x  -x  )  ds 

n-1  n-1  a  n    n  n  a    n 

I  -  2       1     2 

n-1  n 


/2  x(x-x  ) 


Instruction 


Input  to  MRA 


In  MRA 


Transfer 
from  MRA. 


ica  tl 

imr  xl+pbh+c 


(x    -,-x  /I     _)ds  /2  (x      -x  )/l     _ds  /2 

n-1     a7    n-1       n  n-1     a       n-1     n' 


x    x. 


n-1 


x     .,  (x       -x  )/I     _ds  /2 
n-1     n-1     a  '    n-1     n 


B  +n-2 

-r 


B  ,+n-l 


iad     dl+26  +  \   'x       (x     n-x  )/l     _ds  /2  >    x     _  (x    _-x  )/l     _ds  /2 

/    .     n-1     n-1     a  '    n-1     n'  X    v  n-1     n-1     a       n-1     n 

B,  B, 


B  +n-l 

1 


its     dl+26 


B, 


ica     t2 

imr  xl+pblj+2+c 

> 

B  +n-l 

iad     dl+26 

>: 

(x  -x  )/l  ds  /2 

n     a       n     n 


x    x 


n 


(x  -x  )/I  ds  /2 
n     an    n 

x  (x  -x  )/l  ds  /2 
n     n     an     n 


B  +n 


\   r  x     _(x       -x  )/l     _ds  /2     \  'x  (x  -x  )/I  ds  /2 

X    .     n-1     n-1     a       n-1     n'        X    ^  n     n     a"    n.   n 

B-  B„ 


B  +n 


its     dl+26 


I 

B. 


x  (x  -x  )/I  ds  /2 
n     n     an     n 


13k. 


The  vertical  deflection  between  the  ends  of  the  starboard  stanchion 
due  to  the  redundant  column  axial  force  R  was  found  by 


B_ 


X~X  f 

x-xa)(-y^)ds  «  R^ 


(x  -x  )  ds  (x  _-x  )  ds 

J3 a_  n  n-1  a   n 

I  2  1-  2 

n  n-1 


Instruction 

Input  to  MRA           In  MRA 

ica  t3 

(x   -x  )              (x   -X  ) 
n-1  a               n-1  a 

imr  tl 

X  (x  ,-x  )/I  _ds  /2    (x   -x  ) 
n-1  a  '  n-1  n       n-1  a 

B  +n- 

-2                  B  +n-l 

iad  dl+36 

(x   -x  )2/l  nds  /z        \     (x   - 

Transfer 
from  MRA 


B. 


B 


its     dl+36 


B  +n 


\     (x     _-x  )2/l     _ds  /2 

Z. ^    n-1     a     '    n-1     n 


B. 


ica 

tl* 

(x-x  ) 

n  a 

imr 

t2 

x  (x  -x  )/l  ds  /2 
n  a  '    n  n 

dl+36 

B  +n 

-1 

iad 

-N 

"(x   -x  )2/I   d 

B 


(x-x  ) 
n     a 

(x  -x  )2/I  ds  /2 
n     a     '    n     n 


B^+n-1 

>    '(x  -x  )2/I  ds  /2 
z_ *.     n     a         n     n 

Bi 


its     dl+36 


B  +n 
1 


\    '(x   -X    )2/l    ds    /c 

/   v     n     a         n    n' 


135. 

The  vertical  deflection  between  the  two  stanchion  ends  starboard 
side  due  to  mn  was  performed  below. 


52      (x-xj 

0  J  ~T~  ds  "  m0 

B. 


(x  -x  )  ds         (x       -x  )  ds 

n     a         n  .        n-1     a         n 

— ^ ~~z —  +    — ^ ~z — 


n 


n-1 


Instruction 

ica     tl 
nad     t2 


Input  to  MRA 


(x     _-x  )/I     -.ds  /2 
n-1     a       n-1     n 


4-    (x  -x  )/I  ds  /2 
n     a       n     n 


In  MRA 


Transfer 
from  MRA 


(x       -x  )/l     .ds  /2 

n-1     a  '    n-1     n 

(x       -x  )/I     _ds  /2 
n-1     a  '    n-1     n 


+(x  -x  )/l  ds  /2 
n     a       n     n 


B  +n-l 


B  +n 
1 


iad     dl+6            -f    >             (x     _-x  )/I     _  >       (x       -x  )/l       ds    , 

/    x  n-1     a       n-1  / ^      n-1     a       n-1     n/2 

B_  '  RT 

1           +(x  -x  )/l     ds  /2  -1       +(x  -x  )/l  ds  /2 

nann  nann 


its  dl+6 


ict  g21+l 


isp 


B1+n 


B. 


n-1  a   • n-1 


+  (x  -x  )/l  "I  ds  /2 
n  a   nl   n 


(Cycles  back  through  the  stanchion  integration 
program) . 


The  following  sub-routine  performs  the  integrations  around  the 
port  side  which  will  give  the  total  vertical  deflection  between  the  ends 
of  the  port  stanchion.  The  entire  sum  should  equal  zero,  since  the  stanchion 
is  assumed  incompressible. 


136. 


The  instructions  that  follow  perform  the  integration 


B 


Ay 


m  (x-x,  ) 


B 


B. 


ds  = 


Jl,-. 


T. 


m 


m 


n-1 


ds 

2 


Instruction 

Input  to  MRA 

ita  g22+l8 

g22,   icr     pa6 

ica  xl+pb5+c 

Xn-1 

isu     xl+pb£ 

-*b 

its     t3 

imr  a5+pb5+c 

X  ds 
^       n 

imr     cl 

X    1/2 

its     tl 

imr  b3+pb6+pb6 

*    (m//l)n-l 

B  +n 

-2 

iad     dl+£8 

•> 

T 

U/I)n-I(xn-1 

In  MRA 


Transfer 
from  MRA 


n-1 


its  dl+58 


vr: 

% 

(x     ,- 

n-1 

V 

<Vi 

•*b)ds* 

K-i 

-^5^/2 

(x       - 
n-1 

■xb)dan/2 

(m/D 

n-l(xn- 

■l-\)dSn/2 

B  +n-l 

>: 

B3 

m  (x 
n     n- 

•r^^n-iV2 

B  +n-l 

> 

m        (x     . 
n-n     n-_ 

■%  y\. 

ica  xl+pb^+2+c 
isu  xl+pb5> 


+x 


n 


*b 


+x 
n 


Instruction 


Input  to  MA 


In  MRA 


137. 


Transfer 
from  MRA 


(continued) 


its  tit 

imr  a5+pb5+c 

X  ds 
n 

imr  cl 

X  1/2 

its  t2 

x  - 


n'\ 


(\-\)dSn/2 


(x  -x,  )ds  /2 

n  d   n 


imr  b3+pb6+pb6-2     X  (m/l) 


n 


^VvW2 


Bj«-2 


B  +n 


iad     dl+58  (m/l)n_1(xn-xb)ds/2      ^>    U/lMx^Ms^ 


B. 


B. 


its  dl+58 

OUT 

ca 

g22+l6 

ts 

g22+7 

su 

c7 

ts 

g22+l6 

IN 

Bj!n 


^(mA)n(xn-xb)dsn/2 

B_ 


imr  b3+pb6+pb6-2  imr  b3+pb6+pb6-2 


imr  b3+pb6+pb6-li 


imr  b3+pb6+pb6-2 


imr  b3+pb6*pb6-l| 


The  following  section  of  the  sub-routine  performs  the  integral 


B, 


yU-Xj^) 


'o 


ds  =  p 


0 


Vi^n-A1  ,  r^vV 


B. 


n-1 


n 


ds 


n 


38. 


Instruction 

Input  to  MRA 

ica     tl 

(x     ,-x,  )ds  /2 
n-1     d       n 

imr  d3+pb£+c 

l(1/«ri 

its     tl 

imr  yl+pb5>+c 

x7n-l 

B  +n- 

•2 

iad     dl+18 

> 

y  (x  -x.  )/l  ds 

In  MRA 


ir&usier 
from  MRA 


(vrV/WV2 


B  +n 


>        y  (x  -x,  )/I  ds  /2  \       y    _  (x     ,-x  )/I     ,ds  /S 

X ,   'n     n     b  '    n     n  X        *  n-1     n-1     b       n-1     r/ 

B.  B„ 


its     dl+18 


B0+n-l 


/       y       (x     ,-x  )/I     ,ds  / 

X k      n-1     n-1     b       n-1     n 


'2 


ica     t2 

imr  d3+pb5+2+c 

its     t2 


X   1/1 


imr  yl+pb5+2+c  xy, 


(x  -x.  )ds  /2 

no       n 

(x  -x,  )/l  ds  /2 
n     d       n     n 


n     o  •     • 


n 


y  (x  -x.  )/I  d3  /2 
n     n    P       n     n' 


1 


+n-l 


iad     dl+18 


B£L 


y     ,(x       -x  )/l     _ds  /2 

^     n-1     n-1     b       n-1     n 


yn(x-x)/lds/2 

_^  n     u       i»     n 


its     dl+18 


>    y  (x  -x  )/l  ds  /S 

X     .  n     n    b       n     n' 


139. 


The  following  section  of  the  sub-routine  performs  the  integration 


E, 


x(.x-x,  ;  \ 


% 


ds 


%A^ 


B 


x  (x  -x,  )    ,  x     n  (x       -X,  ) 

n     ii     b     as  n-1     n-1     P 


n 


n-1 


as 


Instruction 

ica  tl 

imr  xl+pb5+c 


Input  to  MRA. 


In  YRk 


Trans 
from  MRj 


x  x 


(x   -x,  )/I  _ds  /2     (x   -x.  )/l  -ds  /2 
n-1  t>   n-1  ir        n-1  d   n-1  n 

x   (x   -x,  )/l  -ds  ,/2 

n-1  n-1  b   n-1  n 


n-1 


B  +n-2 


B  +n-l 


iad    dl+28  ^    Xn(VXb)/Inds/2     ]T  Vl(lri^>/1;riV 

B.  B„ 


B  +n-l 


its     dl+28 


x    _(x     ,-xJ/l     _ds  /2 

n-1     n-1     d  '    n-3     n 


ica     t2 

imr  xl+pb5+2+c 

iad     dl+28 


X    x 


n 


(x  -x,  )/I  ds  /2 

no       n     n 

x  (x  -xj/l  ds  /2 
n     n     b  '    n     n 


2 


+n-l 


B  +n-l 


B. 


x     -,(x  -x.  )/l     „ds  /2       y        x     ,  (x  -x,  )/I  ds  /2 
j,      n-1     n    b  '    n-J.     n         Z. ^     n~l     n     ®       n     n 


B  +n 


its     dl+28 


x  (x  -x,  )/I  ds  /2 

_^  n     n     D  '    n     n 


lho. 


The  following  section  of  the  sub-routine  forms  the  vertical 
deflection  caused  by   the  axial  stanchion  force  redundant  S. 


B, 


Ay  = 


(x  -x.  Kx-x,  ) 

\„     ^  ds  «  S 


II 


B. 


^-f  (x  -x,  )  ds     (x  t-x,  ) 


B. 


2-, 


I     2 


n 


n-1 


ds 


n 


Instruction 

ica     t3 
imr     tl 


Input  to  MRA 


W*b 


In  MRA 


K-rV 


from  MRA. 


X(x     n-xJ/l     ..ds  /2  (x    ..-xj2/!    ,ds  /2 

n-1     b  '    n-1     n'  n-1     b     '    n-1     n' 


iad     dl+1+8 


B  +n-2 

^       7- 


(x     0-xJ2/I       ds  /2 
n-2     b         n-2     n 


B  +n-l 


(x     -,-xJ2/I     -ds  /2 

n-1     d     '    n-1     n 


B. 


B. 


its     dl+li8 


B  +n-l 
3 


1?.. 


(x    ,-x.  r/I     nds  /2 

n-1     b  n-1     n' 


ica     tli 

(vV 

imr     t2 

x  K-^W2 

B  +n-l 

iad     dl+U8 

Zl  (xn-r^)2/In-ld 

(Vxb) 

(x-x)2/Ids/2 
n     d         n     n 


B  +n 

-* r 


B. 


>       (x       -x  )2/I  ds  /2 
X    >      n-1     b     '   n     n' 

B„ 


its     dl+l;8 


E i+n 


(x  -or  r/I  ds  /2 

_*   n    b         n     n 


B. 


ua. 


The  following  section  of  the  sub-routine  finds  the  vertical 
deflection  caused  by  the  redundant  moment  nu. 


Ay  =  m 


0 


B, 


(«„) 


B 


ds  =  m 


0 


k 


B. 


(V*b}     ffn       (*n-lV 
2 


I 
u       n 


n-1       J 


dr 


n 


Instruction 

ica  tl 
iad  t2 


its  dl+8 


Let  g22+l 


Input  to  MIA 


In  MRA 


Transfer 
from  MRA 


(vV/W2  (vrxb)/IB-ids/2 


B  +r 


B  *n 


iad    dl-B  J^  [(Vl^ri  zlfVr^'^n-^V2 

B3      *(V^)/In]dSn/2       B3  +fen-^)/Indsn/2] 


B  +n 

-5 


E 


I  [<vrV/Vi 

♦xn-VJn]  dV2 


(Cycles  back  to  g22+l). 


The  redundant  stanchion  forces  affect  the  deflections  for  the 
web  as  a  whole  in  much  the  same  manner  as  the  "three  basic  redundant  forces 
affect  the  deflection  around  the  two  stanchion  sections.  An  analysis  of 
Figure  VI  shows  that  the  integrals  of  p^,  q  ,  and  m  for  the  traverse 
around  the  stanchion  loop  is  the  same  as  the  integrals  of  R  and  S  about 


12|2. 


the  entire  web  girth.  Therefore,  the  coefficients  of  the  redundants  which 
are  the  integrals  may  be  used  in  several  places .  As  an  example  the  coeffi- 
cient of  the  deflection  magnitude  due  to  m_  between  B  and  B  is  the  same 
as  the  coefficient  of  the  angular  rotation  magnitude  due  to  the  redundant 
R  about  the  entire  girth.  Moment  due  to  R  is  R(x-x  )  and  the  rotation 

cL 

due  to  the  moment  is 


R  \     xa     ds  =  R  ]     ^  a  ds   +  J   —J —  ds  +  /      a  ds. 

Bx  B2 

The  deflection  between  B     and  B     d\.ie  to  m     is 


(x-Xg) 
mo  |-T^ds' 


The  following  instructions  utilize  the  computed  coefficients  to 
fill  similar  values. 

Instruction  Input  to  MRA  In  MRA  Transfer 


from  MRA 


ica 


dl+6 


B. 


r3 

(x-x  )/lds         +  1   (x-x  )/Ids 

cl  J  3. 


B3 


Its  dl+30  j      (x-x  )/Id 


Bi 


B„  B„ 

2  p2 


ica  dl+16  /   y(x-x  )/lds        +  J      y(x-x  )/lds 


Bi  Bi 


Instruction 


Input  to  MRA 


In  MRA 


(continued) 


1U3. 


Transfer 
from  MRA 


its     dl+32 


r2 

♦  y(x-x  )/lds 

i.     a 


ica 


dl+26 


B 

C2 

+    I       x(x-x  )/lds 
B_ 


its     dl+3i| 


x(x-x  )/lds 

3, 


ica     dl+8 


B, 
J       (x-x^/lds 


its     dl+UO 


(x-xb)/Ids 


B. 


ica     dl+18 


B, 

J      yCx-x^/Ids 
B_ 


its     dl+2|2 


ica     dl+28 


3, 


/      x(x-x  )/lds 
B_ 


B, 


7(x-xb)/lds 


B. 


lU3. 


Instruction  Input  to  MRA  In  MM  Transfer 

from  MRA 


ica 


dl+8 


(continued) 


ica  dl+26  +  J   x(x-xa)/Id 


% 


J      (x-x^/Ids 


its  dl+UO 

ica  dl+18 

B3 

its  dl+U2 

ica  dl+28 

B| 

J      x(x-x  )/Ids 


% 


its  dl+32  +  I   y(x-x  )/Ids 


bi 

its     dl+3ii  I      x(x-x  )/lds 

i     * 


/ 


(x-x,  )/Ids 


h 


J      y(x-xb)/Ids 


lldu 


Instruction  Input  to  HRA  In  MRA  Trans  f^ 

from  MRA 

(continued) 

B, 
its  dL+hk  |   x(x-x  )/Ids 


isp 

The  following  sub-routine  reduces  the  matrix  of  deflection 
equations  "by  Crout's  method  and  finds  values  for  all  redundants.  The  sub- 
routine, designated  by  S2,  is  a  routine  contained  in  the  library  of  the 
computer  laboratory.  The  library  routines  have  been  tested  and  tried 
and  are  known  to  be  correct.  For  this  reason  only  the  method  of  entering 
the  routine  will  be  given. 

Instruction  Input  to  HRA  In  KRA  Transfer 

from  MRA 

ita  g23+5 

OUT 
g23,  SP  S2 

dl  Address  of  first  column. 

5>  Order  of  matrix. 

1  Number  of  columns  of  constants 

IN 
isp 

The  instructions  below  alter  the  vertical  and  horizontal  forces 
to  include  the  redundant  forces  p  ,  a   ,  R  and  S.  Note  that  each  redundant 
is  applied  as  a  load  at  its  point  of  action. 


\  -  \  -  P0      \  -  \  -  %  *   E  +  S- 


U£. 


Instruction 

ita     g2lx+25 
g2lj.,   ica     e5> 

isu     dl+£U 

its     e5> 


ica  e3 
isu  dl+52 
its  e3 

ica  e^+pbl). 

isu  dl+56 
its  e5+pblj. 


Input  to  MRA 


In  MRA 


+V 


n 


+V 


n 


"o 


V  - 
n 


% 


The  following  instructions  form  H-p 


0 


+H 


VJ 


+11 


w 


-p, 


H  "Po 
w  *0 


The  following  instructions  form  V-R 


+V 


ii 


+¥ 


-R 


ica  e5+pb5+pa5+pa5>  +V 


ii 


n 


iad     dl+56  +R 

its  e5>+pbl;+pa5+pa5> 

ica  e5>+pb6+pb6  +V 

isu     dl+58  +S 

its  e5>+pb6+pb6 

ica  ob6+pb6-pa6-pa6    +V 
iad  dl+5'3  +S 

its  pb6+pb6-pa6-pa6 


n 


¥  -R 
n 


+V 
n 


V  +R 
n 


+V 
n 


V  +S 
n 


+V 
+V+S 


Transfer 
from  MRA 


V  -o. 
n  V 


H-P 


V  -R 
n 


V  +R 
n 


V  +S 
n 


V+S 


12*6. 

The  Whirlwind  instructions  that  follow  change  the  address  of 
g2U+l6  and  g2U+l8. 

Instruction  Input  to  MRA  In  MRA  Transfer 

from  MRA 

OUT 

ca  g2h+l6  ica  e£+pb6+pb6-2  ica  e3>+pbo+pb6-2 

su     c7  -2  ica  e5+pb6+pb6-l| 

ts  g2li+l6  ica  e5+pb6+pb6-I) 

td  g2l*+l8  e5+pb6+pb6-!i 

IN 

ict  g2U+l6 

The  following  sub-routine  finds  the  stresses  at  each  station  due 
to  the  axial  loads  (p),  the  shear  loads  (q)  and  the  bending  moments  (m). 
After  the  various  stresses  have  been  formed  the  direct  stresses  and  the 
maximum  compression  and  tension  bending  stresses  are  combined  to  form  the 
total  stresses  at  the  flange  and  at  the  shell  surface.  The  relations!"  '  s 
listed  below  are  utilized  in  the  sub-program. 

p  =  H  cos  /  -  V  sin  / 

q  ss  H  sin  $  -   V  cos  fS 

Direct  stress  0     =  p/A 
a 

Shear  stress  CTT  =  q/A 
Bending  Stresses  =  my/I 

m    shell  =  m/I  {±T   .   +t  ) 
m    '     I..A.  s 

fm=  flangem  =  m/I  (y^-d^) 


T-  =  rr  +^  f 


Ili7, 


The  following  instructions  use  the  above  relationships  for  the 
starboard  side  of  the  frame. 


Instruction 

ita  g27 
g25,  ica  c6 
its  t8 
its  t7 
icr  pal 

ica  t8 

iad  e£+c 

its  t8 

imr  al+c 

its  tl 
ica  t7 

iad  e3+c 
its  t7 
imr  a2+c 
isu  tl 


Input  to  MRA. 


+0.0 


In  MRA 


n-1 


+v 


X  sin  ^ 


n-1 


+H 


X  cos  <f> 
-V  sin  ^ 


Transfer 

from  MELA. 


+0.0 


n-1 


n 


n 
)     V  sin  <j> 


n-1 


+  1 


II 


n 


+    ?     H 


n 


+0.0  in  t8 
+0.0 


n 


v 


n 


in  $ 


II 


/      H  cos  <j> 
/    H  cos  <f>  -  y   V  sin  ^  =  p 


-r' 


Instruction 


idv  d2+c 

inrr  el 

its  xl+c 

ica  t8 

iinr  a2+c 

its  tl 

ics  t7 

imr  al 

isu  tl 

idv  d2+c 

iinr  el 

its  sl+c 


Innut  to  mk                         In  I-1RA 

-•    "jsfeT" 

fro: 

(continued) 

f-A                                                p/A 

X  3                                          8  p/A 

n 


+    >     V 


X  cos 


/ 


n 


X  sin  / 


n 


V  cos  <f> 


tA 


X   8 


p/A  =cr 


n 
/     V  ccs 

|t 

n 

> 

V 

cos 

^ 

n 

n 
/  i  H  sin 

^ 

n 

/ 

n 

> 

COS 

/  = 

:    O 

q/A 

8   q/A 

8  q/A  ~07 


Note  that(7T  was  formed  using  the  entire  area  as  effective  in 
shears.  This  is  an  error  which  should  be  corrected  since  only  the  area 
of  the  web  is  effective  in  shear. 


Ics     cll+c 


-d 


N.A. 


"dN.A. 


111*. 


Instruction 

Input  to  I 

■IRA 

In  MRA                        Transfer 

from  MRA 

(continued) 

isn 

b£+c 

-t 
s 

-(dH.A.tV 

imr 

cl5 

X  12 

-12<dN,A.+V 

imr 

b3+c 

X   m/I 

-12  ra/KdK<A>+ts)  =cr; 

imr 

el 

X  8 

-96  m/KdH>A-+ts)  =  -  cr 

iad 

xl+c 

+ 
a 

-fcr  ~<r)  =cr 

sm     a           s 

its 

a^+c 

err 

ica 

bl+c 

+d 
w 

+d 
w 

iad 

bli+c 

+tf 

d  +t. 
w     f 

isu 

dl+c 

~*XJL. 

(W-W 

imr 

cl5 

X  12 

^WW 

imr 

b3+c 

X  m/I 

12  m/Kd^tf-dH>A-) 

imr 

el 

X    8 

96  m/I(d  +t -d     ,    )  =  CT 
w    f    N.A.           fm 

iad 

xl+c 

+ 
a 

fm     a         f 

its 

e5+c 

(Ti 

ict 

g25+3U 

The  following  instructions  calculate  the  stresses  around  the 
port  sidee 


icr  pbl 
g26,  ica  t8 


n-1 


X> 


n-1 


+  >   V 


iad  e5+pal+pal+c      +V 


n 


Ij 


imr  al+pb2+pb2       X  sin  $ 


n 


sin  ^  /  V 


ica  t? 


n-1 


+  ?      H 


n-1 

z: 


H 


isu  e3+pb2+pb2        +H 


n 


H 


i5o. 


Instruction 


Input  to  MM 


In  MRA 


Transfer 
from  MRA 


its  t8 


n 


V 


its  tl 


n 


♦sin  /  >   V 


its  t? 


n 


7I 


imr  a2+pb2+pb2       X  cos  / 


iad  tl 


n 


+sin  <f>  /      V 


n 


cos  </>  /     H 

n  n 

cos  /  /      H+sin  /  /  V  =  p 


idv  d2+pb2+pb2 


^-A 


P/A=C^ 


its  xl+pal+pal+c 


a 


15>1 . 


Instruction  Input  to  MRA.  In  HRA 

(continued) 

n  n 

ica     t8  />     V  /. 

imr  a2+pb2+pb2  X  cos  <f> 


V 


n 


cos 


tz  J 


Transfer 
f rora  MM 


its     tl 


ica     t7 


n 

z 


H 


n 

z 


H 


cos  / 


n 

Z* 


imr     al+pb2+pb2  X  S^J]-  $ 


sin 


n 

/  Zh 


isu     tl 


n 
-cos  <f>    y 


V 


n  n 

+  (sin  jz£    />      H-cos  ^    /      V)  «  q 


idv  d2+pb2+pb2 


its  sl+pal+pal+c 


ics  dl+pb2+pb2 

"dN., 

isu  b5+pb2+pb2 

-t 

s 

imr  b3+pal+pal+c 

Xm/I 

isu  xl+pal+pal+c 

-(T 

a 

q/A  =  <7^ 


^o 


-d 


N.A. 


-(dM.A.+V 


N.A.   s 


-(  ♦    ) 
a  sm 


its  a5+pal+pal+c 


&t 


152, 


Instruction 


Input  to  MRA 


In  MRA 


Transfer 
from  MRA 


(continued) 


ics  dl+pb2+pb2 

_dN.A 

iad  bl|+pb2+pb2 

+tf 

iad  bl+pb2+pb2 

+d 
w 

imr  b3+pal+pal+c 

Xia/I 

iad  xl+pal+pal+c 

-<r 

_dM.A. 
"dN.A.+tf 

(vtrtu.) 


m. 


'/I(WdN.A.)  "<?S. 


its  e5+pal+pal+c 


f 


OUT       (The  following  instructions  change  the  addresses  of 
the  above  program). 


ea  g26+3 
su  c? 
ts  g26+3 
ts  g26+17 
ca  g26+7 
su  c7 
ts  g26+7 
ca  g26+9 
su  c7 
ts  g26+9 
ts  g26+lli 
ca  g26+ll 
su  c7 
ts  g26+ll 


imr  al+pb2+pb2 
-2 


iad  e3+pb2+pb2 
-2 

imr  a2+pb2+pb2 
-2 


idv  d2+pb2+pb2 
-2 


imr  al+pb2+pb2 
imr  al+pb2+pb2-2 


imr  al+pb2+pb2-2 


iad  e3+pb2+pb2 
iad  e>pb2+pb2-2 

imr  a2+pb2+pb2 
djnr  a2+pb2+pb2-2 


iad  e3+pb2+pb2-2 


imr  a2+pb2+pb2-2 
imr  a2+pb2+pb2-2 


idv  d2  +pb2+pb2 
idv  d2+pb2+pb2-2 


idv  d2+pb2+pb2-2 


1*3. 


Instruction 


ts  g26+l8 

ca  g26+20 

su  c7 

ts  g26+20 

ts  g26+25 

ca  g26+21 

sii  c7 

ts  g26+21 

ca  g26+26 

su  c7 

ts  g26+26 

ca  g26+27 

su  c7 

ts  g26+27 

IN 
ict  g26 
g27,  isp 


Input  to  MRA 


In  MRA 


Transfer 
from  MRA 


(continued) 


ics  dl+pb2+pb2 
-2 


isu  b5+pb2+pb2 
-2 

iad  bl;+pb2+pb2 
-2 

iad  bl+pb2+pb2 
-2 


idv  d2+pb2+pb2-2 


ics  dl+pb2+pb2 
ics  dl+pb2+pb2-2 


ics  dl+pb2+pb2-2 
ics  dl+pb2+pb2-2 


isu  b£+po2+pb2 
isu  b5+pb2+pb2-2 

iad  bii+pb2+pb2 
iad  bU+pb2+pb2-2 

iad  bl+pb2+pb2 
iad  bl+pb2+pb2-2 


isu  b£+pb2+pb2~2 


iad  b!|+pb2+pb2-2 


iad  bl+bb2+oh2-2 
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